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ABSTRACT OF THE DISSERTATION 
Electromagnetic Plasma Wave Propagation 
Along a Magnetic Field 
by 

Craig Lee Olson 

Doctor of Philosophy in Physics 
University of California, Los Angeles, 1970 
Professor Burton D. Fried, Chairman 

Part 1: Steady-state excitation of transverse plasma waves along a 

magnetic field. 

The linearized response of a Vlasov plasma to the steady-state 
excitation of transverse plasma waves along an external magnetic field 
is examined. Previous research is reviewed. Assuming a delta- function 
excitation mechanism, and performing a detailed Vlasov-Maxwell equation 
analysis using Fourier-Laplace transforms, the plasma response is found 
to consist of three terms: a branch-cut term, a free-streaming term, 

and a dielectric-pole term(s). These terms are examined analytically, 
and numerically evaluated for a case in which the driving frequency 
((u.) is slightly below the electron cyclotron frequency (as ) . In addi- 
tion to the least-damped pole term, it is shown that the free-streaming 
term is always significant and the branch-cut term is significant when 
ojj ~ !*) ce . The infinite sequence of pole terms (that result from the 
infinite number of roots of the appropriate transverse dispersion rela- 
tion) is shown to be negligible except at positions very close to the 
place of excitation. Effects of Krook model collisions are investigated 



•J 

* 


as well as effects of a finite-width excitation mechanism. 

Part II: Steady-state transverse plasma wave echoes. 

The phenomenon of plasma wave echoes, introduced by Gould, Malm- 
berg, O'Neil, and Wharton for the case of longitudinal electrostatic 
waves, is extended to the case of transverse plasma waves that propagate 
along an external magnetic field. It is shown that a transverse echo 
results in lowest order only when one excitation is transverse and the 
other is longitudinal. For this case, the second-order (nonlinear) 
plasma echo response is computed from the Vlasov-Maxwell equations up 
to an integral over the velocity variable v . Transverse echo charac- 
teristics are discussed and several experiments are suggested. The 
integral over v is evaluated by the method of steepest descent; the 
results are explained physically and also evaluated numerically for some 
specific cases. Several extensions are considered: (1) effects of 

Fokker-Planck collisions and finite-width excitation mechanisms, 

(2) effects of temperature anisotropy (T, t T ) , and (3) effects of 
propagation just off-axis (k a ^ 0). Lastly, transverse echoes for the 
case of no external magnetic field are examined. 

In addition, several appendices concerning Part I and Part II are 
presented, the most important of these being (1) a thorough investiga- 
tion of the roots of the transverse dispersion relation (k^ | | for 
real w and cony lex k, and (2) the method of steepest descent as used 
in evaluating the branch-cut and phase-mixing integrals that occur in 
Part I and Part II. 
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I. STEADY-STATE EXCITATION OF TRANSVERSE PLASMA WAVES 


ALONG AN EXTERNAL MAGNETIC FIELD 


1. Introduction 

A sizable portion of the plasma physics research done to date 
involves the theory of waves in plasmas. Of fundamental interest are 
calculations or experiments concerning the damping (or growth) of the 
basic types of wave motion a plasma can support. Two types of problems 
are readily discernible: (1) the initial-value problem in which a wave 

is imposed on the plasma at some specific time, and it is desired to 
know how the plasma response decays (or grows) in time; and (2) the 
steady-state problem in which a wave is continuously excited at some 
point in a plasma, and it is desired to know how the plasma response 
decays (or grows) in distance. We will be concerned with the steady- 
state problem since it should be more readily realizable experimentally. 

Purely longitudinal waves which experience spatial Landau damping 
have been examined in theory extensively (since Landau's original 
paper^) by Gould^ and by Johnston, ^ et al.^’ 5 ^ and in experiment 
by Wong, D'Angelo, and Motley, et al.^ However, the corresponding 
work for purely transverse waves that propagate along an external mag- 
netic field (and therefore experience cyclotron damping when the driving 
frequency is near the cyclotron frequency) has received relatively little 
attention to date. 

Thus the problem to which we address ourselves is that of determin- 
ing the response of a Vlasov plasma (a hot correlationless plasma) to a 
localized steady-state excitation of transverse electromagnetic waves 
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that propagate along an external magnetic field. We are specifically 
interested in obtaining the form of the response when the excitation 
frequency is just below the electron cyclotron frequency, i.e., when 
electron cyclotron damping is significant. One primary objective of 
this theoretical calculation is to obtain the cyclotron-damped response 
and see if it is adequately represented by consideration of just the 
"least-damped root" of an appropriate transverse dispersion relation. 

We will find in general that it is not. 

Another motivation for this work is that the response of a single 
transverse excitation should be well understood before embarking on 
studies of the nonlinear response of two spatially separated excita- 
tions, such as occur in Part II of this dissertation where we shall 
investigate transverse plasma wave echoes. The production of an echo 
(as will be explained in Part II) depends on the free-streaming term of 
the first-order distribution function produced by the first of the two 
excitations. Thus in computing the response of a single transverse 
excitation now, we will pay particular attention to the transverse 
free-streaming waves that are produced. 

We present a brief summary of previous research. First we review 
work on the related problem of the steady-state excitation of longitu- 
dinal plasma waves. Experimental verification of spatial Landau damping 
was reported by Wong, D'Angelo, and Motley^ who compared their exper- 
imental results with a simple theory that involved just the least-damped 

root of the appropriate longitudinal dispersion relation. Use of just 

(2) 

the least-damped root in this case was questioned by Gould who rigor- 
ously investigated the first-order (linear) response for the problem 
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assuming an idealized dipole-grid excitation mechanism. Gould's results 

showed that the true response was indeed more complicated than that 

obtained from considering just the least-damped root and that the 

response was approximately exponentially damped only over a restricted 

region in space. Further theoretical work including the effects of 

finite spacing of the pair of excitation grids and also a calculation 

(3) 

of the second-order (nonlinear) response was performed by Johnston. 

The longitudinal dispersion relation mentioned above has an infinite 

(7) 

number of roots. Fried and Gould have examined in detail the loci of 

the roots of this dispersion relation for real k and complex w . 

Their results are pertinent to initial value problems in which a fixed 

real k is imposed at a specific time and then the initial perturbation 

f 81 

decays (or grows) in time. Similarly, Kuehl, Stewart, and Yeh^ J and 
Derfler and Simonen v J have examined the roots of the longitudinal dis- 
persion relation for real u> and complex k . Their results are perti- 
nent to steady-state problems (as discussed above) in which a fixed real 
u> is continuously excited and the resultant spatial perturbations decay 
(or grow) in distance. We shall comment more on their results in 
Appendix D. 

We now consider the transverse excitation problem which we find 
has only been partially investigated in the past. The principal paper 
is that of Shafranov^ 0 ^ who considered a related problem, the half-space 
problem. (The full-space problem that we shall consider is different 
than the half-space problem in that, in addition to other effects, free- 
streaming waves are produced by those particles that pass through the 
excitation region. In the half-space problem particles never pass 
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through the excitation region but are reflected at the half-space 
boundary where the excitation is applied.) Shafranov obtained a branch- 
cut term and a pole term for the response (as we shall find also). 
However, we note the following errors and inadequacies in his work: 

(i) The saddle-point contours for the k (actually kc/u ) inver- 
sion integral are drawn to the wrong asymptotes. Also the primitive 
contour was closed in some unexplained manner whereas it should really 
be deformed away from its initial position. In addition the resultant 
saddle point contributions contain phase factors that do not occur in 
usual saddle point theory. 

(ii) Only one root of the appropriate transverse dispersion rela- 
tion is considered whereas in fact there are an infinite number of them. 

(iii) Useful values of the cyclotron damped root are not given. A 
cubic equation is given, which if solved, gives the value of the cyclo- 
tron-damped root at w * a) (u is the electron cyclotron frequency) . 

Ce ce 

The only other expression given for the value of the cyclotron- damped 
root is an asymptotic result which cannot be derived rigorously as we 
shall show in Appendix C. 

(iv) No estimates are given for the range of validity of the branch- 
cut term. In addition it is stated that the branch-cut term is the most 
important term at large distances from the place of excitation for 
w < oi < w D (w D is the right-hand cutoff frequency); we shall find 
that the branch-cut term is significant only very close to the place of 
excitation. 

Thus Shafranov 's work represents only a first attempt at solving 
the problem. 


4 



Other research on the half-space problem includes Sach's^^ work 
on the effects of collisions on ion cyclotron waves. In this work two 
least-damped roots were obtained by keeping two terms in the asymptotic 
expansion (of the plasma dispersion function that occurs in the disper- 
sion relation) with the result that the roots are not valid unless the 

collision frequency is appropriately high. Another related work is that 

( 12 ) 

of Platzman and Buchsbaum on the reflection and transmission of 
transverse waves into a semi -infinite plasma. In this work a Lorentzian 
zero-order distribution function was used (to simplify numerical calcula- 
tions) and attention was focused on determining the effects of collisions 

on the reflection and transmission coefficients. Lastly we note that in 

( 13 ) 

an experiment directly related to our problem, Crawford et al. have 
attempted to measure spatial electron cyclotron damping but were unsuc- 
cessful because of high collision rates. The results they did obtain 
were compared with a simple, least-damped root theory. Clearly a thor- 
ough rigorous theoretical investigation of the steady-state transverse 
excitation problem is still needed, and that is what we shall present 
shortly. 

The transverse dispersion relation mentioned above has previously 
been only partially examined. (We state explicitly that, as used hence- 
forth, the transverse dispersion relation refers to wave propagation 

( 14 - 19 ) 

exactly along the external magnetic field.) Many authors have 

derived the transverse dispersion relation but all of them considered 
real k and complex co and then, at most, obtained an asymptotic 
expression for the value of the least-damped root. Recently Kamimura 
and Hasegawa^*^ have presented numerical results concerning the value 
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of the least-damped root (for real k and complex u> ) for several 

choices of the pertinent parameters. At present, however, there exists 

no work that classifies and examines the infinite number of roots of 

the transverse dispersion relation for real k and complex w . 

Previous work on the roots of the transverse dispersion relation 

for real u> and complex k is even scarcer. Only isolated instances 

may be found such as in Stix^ 17 ^ and in Shafranov^ 10 ^ (as mentioned 

above) j in both of these the value of the cyclotron-damped root is given 

exactly at u> = w . Also, in Shafranov an asymptotic result is given 
cc 

for w far from u> , a result which, as we shall show in Appendix C, 

c© 

cannot be derived rigorously and is therefore incorrect. Crawford et 
al.^ 13 ^ have presented some numerical results for the cyclotron- damped 
root (for u) 2 /tu 2 = 40 where ui is the electron plasma frequency) 

pG C© P© 

but no detailed study of the behavior of this root was made. Thus a 
classification and thorough investigation of the infinite number of roots 
of the transverse dispersion relation for real w and complex k is 
still needed and that is what we will present in Appendix C. 

We have performed a detailed theoretical investigation of the 
steady-state transverse excitation problem assuming an external excita- 
tion of the form 


E (?,£■) = E.(z)fe x co s(cj l t) + e Y s/v(<^t)] 
— exr ' L 


( 1 . 1 ) 


E f (z) = 

Before commencing with our presentation we summarize the most important 
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results that will be obtained: 

(1) From a rigorous derivation starting with the Vlasov-Maxwell 
equations we find that the first-order (linear) current response consists 

> of three terms as given in Eq. (2.85); a branch-cut term, a free-stream- 
ing term, and a dielectric-pole term. 

(2) Numerical evaluation of these terms for a specific case in 

which u. is slightly below u ^ (with to /o> ^ = 0.4, c/a = 1120 ; 

1 ce pe ce 

these being typical values for the UCLA Q machine) shows that, in addi- 
tion to the pole term, the branch- cut and free-streaming terms have sig- 
nificant amplitudes and penetration distances (see Figs. 7a-f). Thus 
the plasma response is not correctly given by consideration of just a 
"least-damped root." 

(3) An analytical investigation of the three terms shows that the 

penetration lengths of both the branch-cut and free-streaming terms 
would approach infinity as approached u> ce were it not for the 

fact that certain dielectric functions are present (in these terms) 
which effectively shield these terms and actually cause their penetra- 
tion lengths to go to zero as cu. approaches w , provided 

X c© 

(u) 2 /ui 2 )(c/a) » 1 . If (oo 2 /w 2 ) (c/a) « 1 then no such shielding 
p© c© p© c© 

effect occurs. In addition we show that the infinite number of residue 
terms (that occur because the transverse dispersion relation has an 
infinite number of roots) is negligible except at distances extremely 
close to z - 0. 

(4) The physical nature of transverse free -streaming waves and of 
the free-streaming term in (2.85) is discussed. Explanations are given 
for (i) why transverse free -streaming waves have negative phase 
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velocities, (ii) why transverse free-streaming waves experience no cyclo- 
tron damping, (iii) why the free-streaming term in (2.85) has an 
accompanying pole contribution for u)j < 0 or > w ce but not for 

0 < to. < a) , and (iv) why the free-streaming term in (2.85) always 
1 c© 

diverges at z = 0 and how this divergence can be removed. 

(5) Effects of collisions on the three terms in (2.85) are examined 
using the Krook model. It is shown that previous work by Crawford et 

ai ( 13 ) (concerning the locus of points where cyclotron damping equals 
collisional damping on a T g vs. plot) is in error (see Fig. 

10). The correct loci are then computed (see Fig. 11). From the results 
it is shown that cyclotron damping measurements (with negligible colli- 
sion interference) are just barely feasible on the UCLA Q machine (for 

n = 10 cm ). In general, cyclotron- damping measurements would be 
6 

more feasible at higher energies (T >0.2 eV) and lower densities 

(n < 10 10 cm' 3 ). 

© 

(6) Effects of a finite-width excitation mechanism on the three 

terms in (2.85) are examined. In place of (1.2), we consider two choices 
of Ej(z) [square-shaped in (7.3) and Gaussian-shaped in (7.11)]. Both 
choices yield qualitatively the same results which are: (i) the diver- 

gence of the free-streaming term at z * 0 is removed, and (ii) that when 
the width of the excitation region is greater than the characteristic 
length of a specific term in (2.85), then that term is reduced in ampli- 
tude. Thus, under certain conditions it is possible to diminish the 
size of the branch-cut and free-streaming terms while leaving the 
dielectric-pole term essentially unchanged. 

In the course of obtaining the above results we shall refer to 
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Appendices A, B, C, E. Appendix A contains the solution of the linear- 
ized Vlasov equation and Appendix B examines symmetry properties of the 
transverse dispersion relation. A thorough investigation of the roots 
of the transverse dispersion relation for real w and complex k is 
given in Appendix C, as mentioned earlier. Lastly, Appendix E contains 
an explanation of the method of steepest descent and a detailed account 
of the paths of steepest descent used in evaluating branch-cut and 
phase-mixing integrals. Material relegated to the appendices is 
actually text material in the usual sense and it contains many results, 
both analytical and numerical, that should be of general interest. This 
material was placed in the appendices because it must be referred to 
both now in Part I (the single transverse excitation problem) and later 
in Part II (the transverse echo problem). 

We proceed to give the results of our investigation of the single 
transverse excitation problem. 
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2 . 


Derivation of the Linear Response from the Vlasov-Maxwell Equations 


2. 1 Formulation of the problem and method of solution using 
Fourier- Lap lace transforms. 

We wish to compute the response of a hot correlationless plasma in 
an external magnetic field to the spatial excitation of transverse elec- 
tromagnetic waves that propagate along the direction of the external 
magnetic field. The plasma is described by Maxwell's equations, and 
the Vlasov equation which for a plasma in external magnetic and electric 
fields is 


a/ 


.v* 


ax 


/ 


£ -+■ 
TOTfU. 



B 

~roT^t 




i£ 

av 


O 


( 2 . 1 ) 


where f = f (ac, is the distribution function (normalized to unity) 
of the species whose particles have charge q and mass m , c is the 
speed of light, and the subscripts on ,E and are used to emphasize 

that these are the total fields (plasma fields plus external fields). 

/\ 

Taking the zero-order external magnetic field to be ^ = B Q e z and 
using cylindrical coordinates for the velocity v (as shown in Fig. 1) 
we find 


a/ a/ 

)L' ~£x 
ar — 


a y~ 
as * 5 




E +£ -) -~^Vx(BJ-& )|* O (2.2) 

-p exr c - v P exc J a V 


where ft = qB Q /mc (a signed quantity) is the cyclotron frequency, and 
the subscripts (p) and (ext) refer to (plasma) and (external) , 
respectively. 


10 



We assume a perturbation expansion of f , E , and B with 

— p -p 


= £<y) + + ••• 

Mj, ( ~ E / ( r) ^ O + • • • 

(*,*) = B m/ fx,r-) + <2S,r) + * * * 


(2.3) 


Bp (x,t) = 


where 0, 1, 2, refer to terms of zero, first, second order. We take 

E and B to be first order. Then in zeroth order, (2.3) is 
-ext -ext 


- A 




= o 


(2.4) 


which means the zero-order velocity distribution is 


/<>) = £ (v^) 


(2.5) 


In first order (2.3) is 


l± + v .l£-a^ 

it 


a + f_ e -h p +■ p-y* (i? +• b )\' 

9<p on —/ c ' / —cxr ' J 


(2.6) 


To obtain the first-order (linear) response we must solve (2.6) together 
with the appropriate Maxwell equations, 


V *(£,+€"*) = - ^ J^(X+-e*r) 




(2.7a) 


(2.7b) 


where 




( 2 . 8 ) 
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and n is the density of the species under consideration. To handle 
all means of excitation we have kept , ^g Xt > ^ext ^ rst “ 

order terms) in (2.7) but will later consider their effects individually. 

We shall use Fourier-Laplace transforms to solve (2.6) - (2.8). 

For an arbitrary function F(z) the Fourier transform and its inversion 
are defined by 


F(A) 


+ OO 



— oo 


- ikz 

F{z) e 





F(F) e 


+ i 


(2.9) 


Similarly for an arbitrary function F(t) , the Laplace transform and 
its inversion are defined by 


F(tv) 

F(r) 


OO 


= / F ( f) e 


+ t cut: 


+ oo t i ^ 


== { <*& 

J Z7T 

-Oo + / & 


F(to ) e 


- lout 


JF»_ Co>o 


( 2 . 10 ) 


where in the former we have assumed Im w > 0 initially to insure 
causality (then e +lti)t 0 for t 00 ) and in the latter we require 
6 to be such that the w contour is above all singularities of F(w) . 
The combined Fourier-Laplace transform pair and its inversion may be 
written 
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OO + HO 


F(k^) = /-* J <*z k(z,t) C 


-i(k^-<cjr) 


U> „Os> o 


( 2 . 11 ) 


- Oo 


+03 

-Oo +-t<# - 03 


+<;' (kz-Ccjtr) 


(2.12) 


Using (2.11) we shall transform Eqs . (2.6), (2.7) to k, u) space, 
solve these equations in k, w space (a relatively easy task), and then 
use (2.12) to transform back to z, t space (a relatively difficult 
task) . 


2.2 Solution in transform space (k,to) 

Taking Fourier- Lap lace transforms of (2.7a,b) gives 




(2.13a) 


/ x i / v 6c/ / __ . 

-k x +- <9 = — — ( </" + v/- ) — (if + ^ ) 

— exr/ c ' / — e*r' c / ex? ' 


(2.13b) 


where we have assumed Im uj > 0 and 


(r-o) = o 

e Cr^o') = £ extf (r *0) *o 


In (2.13a,b) Jl, £;, and £ are all now functions of k, w . Combining 
(2.13a,b) we obtain 
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(s +*„ r )(*V-«JV =*[*■- 


(2.15) 


Since we are dealing with transverse sources, k*J = k*E = 0 

— — GXl ■€ XI 

which 


ich implies k_*E^ = k’Jj = 0 . Thus (2.15) is 





97ri«j (Z, +-J~exr) 


(2.16) 


Equations (2.13a) and (2.16) summarize Maxwell's equations for our 
problem. 

We now solve (2.5) for fj which will be used later in computing 
Jj . Taking Fourier- Lap lace transforms of (2.6) gives 


(£-V -«j) / - 


a 


H 

20 


= - ? (*) 


(2.17) 


/ W + C * * (f ^c-r)] • 


if 

y o 

2V 


(2.18) 


where we have used (2.13a) to eliminate + J ext • Expanding the 
triple vector product in (2.18) we obtain 


-- - 


f /_ ^zrX^, f S e „ 


**7 o ^ , 


/W 




(2.19) 


where we have explicitly used 

/ = ke z 


(2.20) 
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In Appendix A it is shown that (2.17) has the solution [taking into 


account the necessary periodicity in <j> , fj (4>+2n) = (tj>) , 

fact that k x = 0 ] 

/ 2^ 

/M = 

e - / ° 

where 

Cu - ArV^ 


To perform the integral in (2.21) we need expressions for E* 

g 

— * 9v occur in (2.19)] in terms of 4> . To this end 

useful to change to rotating coordinates. 

In terms of the unit vectors 


A 

e. 


- £($.t*8 y ) 


the velocity vector v is 


v = e + v_ +- e_ v + + 


where 


V 4-^ 


— ( V ± l Vw) 
VI * y 


= -p ( Vj_ cos <?> t i v x s//v 0 ) 


Vz 

“ H Vj - € 


±i<fi 


Also 


and the 
( 2 . 21 ) 


( 2 . 22 ) 

and 
it is 


(2.23) 

(2.24) 


(2.25) 


16 



<? 

iV 


A 


> V + 


/V i 

€ - TvT 


A 

€ 


2 ^ V 2 


where 



j-s'**- 

n >Yl. 


C 

-£ VI Vj_ 



Now writing j: * e_ a E 0 with summation over a = ± i implied, 
noting that 


A 



A 



3 O 


A 

e± 



/ 


we obtain 


f • V 



-/ crpf 

e 


^ ^ — e ~ i<r0 /- v — • 


Thus g(4») in (2.19) becomes 


~ I fz E ext<r ) e 




(/-& 

\ Cu I 


ii 


&Yl 

a, 


There is no 3f Q / 84> term because of (2.4). Performing the 
gration in (2.21) we find 



t(5r*e w J * lr * 

nr> Fz i (ctf-AVg -<r-C2.) 


\ ' />, / 3V + 


O/ / ^ 


«>v 2 


(2.26) 

(2.27) 
and 

(2.28) 

(2.30) 

N. 

(2.31) 
inte- 

(2.32) 
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We may now compute as given by (2.8), 


Noting that 



V 


( 2 . 8 ) 


3 ^ 

Y</v = Kvl v ^e 


-- Vz 


Yi_e 


+ <£ 


/s . 

e V 




(2.33) 


we see that since f ^ (<f>) ^ e - * 01 ^ the only term of v e * 0< ^ d 3 v that 
survives the <j> integration in (2.8) equals e^O^/v^) 2lT • There- 
fore writing * e -a^la we °^ ta ^ n 




[t - lA 

[( / ~'zr)^ x * ^ ^v g . 

(Cu-kw ~ crsi) 


(2.34) 


where we have explicitly indicated a sum over species (e * electron, 
i * ion) . We assume the zero-order distribution function is an 
anisotropic Maxwellian 




i 




(2.35) 


(a = [v z 2 ] 1/2 * [2KT z /m] 1/2 t b = [v^ 2 ] 1/2 = [2 kT /m] 1/2 , k = 
Boltzmann's constant). Then (2.34) reduces to 


^ = s (e+e ) 

/er ^ ^77 V '<r exco-/ 


too 


J- 

- Oo 






(kv -o<j+ crxi) 


(2.36) 
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Returning to Maxwell's equations and writing J 

Cav 

we find (2.16) is 


e J 
-a ext cr 


( £ /ir S)(k<z-<xj) ~ + ^ (2.37) 


We must now decide whether to consider (1) an external field source or 
(2) an external current source: 

(1) External field source (E _ + 0. .J „. = 0) 
—ext —ext 

Using Maxwell's equations for just the plasma fields, (2.37) is 


( ) = Vtt i co <J^ C 


/ er 


/cr 


(2.38) 


and since E total a = E la + E ext a we find using (2.36) that 


E(k,oi) 

t °7/>c cr 


E (*>^) 

exr o- 


(2.39) 


where e T (k,w,a) , the transverse dielectric function, is 


— . f 

^ri^cr) = /-2l_ C^cW)/ 

* rtr» 


"VuU'- £*<-!)] 


( kv^. ~ou + axi) 


(2.40) 


2 2 

(u)p * 4fmq /m is the plasma frequency). 

(2) External current source (J „ / 0, E = 0) 

■ ,_ CXt C Xt 

Using Maxwell's equations for the total fields, (2.37) is 

E (k z c z -<^> z> ) = ^ xr r ) (2.41) 

/cr ' 

and using (2.36) we find 
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(2.42) 


E (*,<*/) 

Total CT 


*•«>') = 


Vzri Cu 
APcP-Ca' 


J (*, to) 

e*r cr 


€ r (£,aJ,<r) 


where e T (k,u>,a) is still given by (2.40). Had we considered 

J - J, + J we would have found 

total o 10 ext o 


J ~ (*,“0 = 

total c r 


to, <r) 


(2.43) 


where again e^k.w.a) is given by (2.40). 

Noting the similarity of the results for an external field source 
or an external current source, we arbitrarily choose to consider the 
former and will therefore continue with (2.39) and (2.40). Then 
specializing to the case of isotropic temperature (T^ * T z ) and per- 
forming the integral over v in (2.40), we obtain 

m 



(2.44) 


where u , u» - IflB /me I , a all refer to the electrons; w , , u . ■ 
pe ce o r 1 

|flB o /Mc| , A all refer to the ions; and Z is the plasma dispersion 
function^*) defined by 


Too _ t 2 


2 ± (D - 

— Oo 


<&r e 

■t-r 


for Im C^0 and its analytic continuation for Im C^0 


(2.45) 
The ± signs 
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refer to different branches of the Z function; the choice of which 
branch to use will be discussed shortly. 

To complete the solution of our problem in k, w space we 
explicitly specify the external field source to be 


E 

— ex*r 


(?,t) 



COS (O'/'O 


+ e y SIN 



(2.46) 


so Uj > 0 corresponds to a right-hand excitation while < 0 corre- 
sponds to a left-hand excitation. Writing = e_ a we find 


E = f/z) e iCrUJ ' T 

exr cr 


(2.47) 


Taking Fourier- Lap lace transforms of (2.47) gives 


E = 

cr (ca-hcru^) 


(2.48) 


We have now solved our problem in k, a) space since we now have 
expressions for f^k.w) , ^ (k,<*>,a) , and E ext o (k,u>) . Equations 
(2.32), (2.39), and (2.48) combine to give 


/ t ( )L, Ca ) ~ 




~ i,r f >£ 


s>r> 


Vi 






(cu l-crc*,') (us-kVg~ trfl) € r C/f^er) 


(2.49) 


We proceed to perform the inverse Fourier-Laplace transform to obtain 
f (z.v.t) to be used in calculating JjU.t) , the final (current) 
response . 
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2.3 Inversion of solution to position space (z,t) 
The inverse transform of (2.49) is 




tr) = 



E t Ch') C 


2 (k?-CuT) 




B - 


• i*"* >A 


(2.50) 


where 6 > 0 is chosen so the primitive to contour, as shown in Fig. 

2a, is above all singularities of the integrand of (2.50). The primi- 
tive k contour lies along the real k axis, as shown in Fig. 3a. 

The proper manipulation of these contours will show us whether to use 

e T (k,oj,a) or e T (k,u,a) in (2.50). 

+ 

Branch-cut considerations 

We reconsider the Z function of (2.45). As defined, the real 

axis of the £ plane is a branch cut of Z(£) • If initially Im ^ > 0 

then Z + (5) should be used and its analytic continuation for Im £ <_ 0. 

If instead Im 5 < 0 initially, then Z_(0 should be used and its 

analytic continuation for Im 5 0 . Regions where Z + (q) are 

originally defined are shown in Fig. 4a. Note that Z + (£) must be 

analytically continued across the branch cut (the real £ axis) into 

the region where Z (£) is originally defined. [A convenient represen- 

+ 

tation of Z + (£) that illustrates these points is 
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FIGURE 2. CONTOUR IN THE u PLANE 


(a) Primitive contour (Im to = 6>0). 

(b) Contour lowered (6-+0) and deformed around pole 

at (only the pole at to=w^ is shown). 

(c) Contour lowered below the real axis (so on 
horizontal part exp(-iwt)^ exp(-|lm w|t) -*■ 0 
as t ») . 





FIGURE 3. CONTOUR IN Tllb K PLANE. 

(a) Primitive contour. 

(b) Primitive contour with specified. Also shown, for 
Re u > w ce , are the branch cuts, the pole at k Q = (u>-u> ce ) / 
v z (for v z >0), and the least-danced roots of e^. + ( k^ ,u>)=0. 

(c) Contour pulled above for z>0. The path of steepest 
descent has been chosen for the branch-cut integral. 

(d) Net result: branch-cut integral plus two residue terms. 




m//\ 

f/imf/t 


Z AC) 

k plane 


z SC) 



V/ui) f / 
/ / / / / / / 


ik plane/ / / 

ZAC) 


k 

use Z AC) 

plane 

use ZAC) 

. * - 






FIGURE 4. BRANCH CUTS IN THE 5 AND K PLANES. 

(a) Branch cut in the ? plane. 

(b) Mapping of the branch cut in the c plane to the k 

plane, k=(co-w )/c a , for Re u > 2 . Also shown 

cc c ce 

is the pole at k =(w-u> )/v for v > 0 . 

1 o ce z z 

(c) Branch cut in the k plane for Re w < u 

ce 

(d) Resultant choice of Z + (s) along the primitive contour 
in the k plane. 
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(2.51) 


z ± (r) - 


/ 


+ oo 



-zr 2 


-/ 2 

■f t'fir %(f) e 


where a + (C) = 0 and o_(£) 
1 
2 


-2 for 

Im C > 0 

-1 

= 0 

0 

< 0 


and P stands for principal part.] 

We set o = -1 for now. Then the arguments of the Z functions 
that occur in (k,io,a) in (2.44) are 



Ar cl. 



to + 


4- A 


(2.52) 


where e, i, refer to electron, ion. Consider C e • The branch cut in 
the ? e plane (the real axis) maps onto the k plane according to 




Co - Cc/ ce 

~r* 

e 


(2.53) 


Recalling that initially Im a) = 6 > 0 and considering a single point 
on the primitive go contour (i.e., an go of the form Re go + i6 ) we 
find the branch cut in the k plane for Re to > go to be as shown in 

Fig. 4b. Note that the points Z, - t ® when mapped onto the k plane 
pinch the origin. If instead we have Re go < Go cg then the branch cut 
in the k plane would be as shown in Fig. 4c. The variable ^ gives 
similar branch cuts (i.e., like Fig. 4b for Re go > -go ^ and like Fig. 
4c for Re go < -go .). In any case note that along the primitive 
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contour in the k plane (the real k axis in Fig. 3a) we now kmw 
that for both £ e and we should use Z + (and therefore ) 

for k > 0 and Z (and therefore e,j, ) for k < 0 . This result is 
shown both in Figs. 4c and 3b where for clarity we have drawn a dashed 
line for that part of the contour where e^. is to be used and left a 
solid line for that part of the contour where £,j, is to be used. 

We may now calculate the integrals over k and to that occur in 
(2.50). We shall perform the k integration and then the w integra- 
tion since we have found this ordering to present fewer complications 
than the reverse ordering. 

The k integration 

In Fig. 3b we show the primitive k contour with the choice of 

e,j, or e,j, specified according to the above discussion. For a 
+ 

typical point on the to contour (we have chosen Re to > to ce and 
Im to = 6 « Re io for this figure) the branch cut is shown as well as 
the pole at 


Zc / — < ^'c<? 



for v z > 0 and also the poles at k^, where 

e r+( = ° (2.55) 

Only the least-damped roots of (2.55) in the plane are considered 

now. [The roots of e T (k ,to) = 0 are discussed in Appendix C.] 

In Fig. 3c we have pulled the contour above for Z > 0 being 
careful to deform the e T side of the contour around the poles at k Q 
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and k. . Only the e„ side of the contour "sees" the poles where 
1 + 

(k 1 ,a)) = 0 . Note that if Re as < W £e then the e. f side of the 
contour would deform around the pole at k Q for v z > 0 . (If v z < 0, 
the pole at k Q lies in the lower-half plane and is missed altogether.) 
We have chosen to deform the e T and tj paths to meet along a 
curved path that will presently be shown to be the path of steepest 
descent (P.S.D.) for this portion of the integral. This integral will 
be called the branch-cut integral but note it is not evaluated along 
the branch cut. As we move the outer part of the contour (the semi- 
circle in Fig. 3c) out to infinity, the contribution to the integral 
from this section of the contour becomes negligibly small (since the 
integrand goes as e^ Z /k ) . 

Thus we are left with the three contributions shown in Fig. 3d; the 
branch-cut integral and the residues from the poles at k Q and kj . 

We assume the latter to be a simple first-order pole, i.e., we assume 




* O 


k-k. 


(2.56) 


Explicitly the results of the k integration are then 


+ °° i_k? 


2 7T 
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■Z77~ 
O 

P.S. D. 




e r(^j <r ) 

+ 






+ 





e 


[0/-01 i \ 

l v z ) 


v.pr-'-'') 



(2.57) 


for 


z > 0 , and where H(v ) is 

//(VJ = / for 

O for 


the Heaviside function 
V ? £ O 

<o 


and 

Ct = S g-/v ( <^k ~ 

@ = S G- A/ (" ) 


Equation (2.57) is valid for a = ± 1 (not just a = -1 as used in 
the preceding discussion). 

The branch-cut integral in (2.57) may be evaluated as follows. 
From (2.44) and (2.51) we find 
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oo 



(kT = - ma 2 kT, = i MA 2 ). We evaluate this integral by the method 
v e 2 ’ 12 

of steepest descent discussed in Appendix E. For z > 0 the electron 
term of the integrand (i.e., the first term within [ J ) has a saddle 


point at 


fr. = e 




<x - SGrH ( fia-Cu + cr<u^. e ) 


and the ion term has a saddle point at 


{(v-c rb'ci') v,} 2 ] 3 


— sg/v ( cc > — < r c # , c ;') 


(2.60) 


The correct choice of the P.S.D. and graphs of it are given in Appendix 
E. Computing the first term of the asymptotic expansion about the 
saddle point we obtain 
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oo 


/- 

o 

f. s. />. 





r / . v? i 

— / ex' 

- 3 

"2 .ex. 


2 e 2 


2 ^pe ^ 


i/J ( X^c 2 - &- 2 ) ( ^ - trjz) 



(2.61) 




+ 


I ^ 


->[■ 


(oj- O-Oj,}} 0£ 2 




ZOL 


VF 


Z "“*? 2 


^ (t4 / c-u?)(*'-criti{t s r ^~~j) e 7 (\.,^<r) € z (\,^) 


which is a good approximation for 


^(c^-uy ce ) 

(X 


£ / 


(2.62) 


(as discussed in Appendix E) . In using the method of steepest descent 
we have implicitly assumed that the integrand of (2.58) had no poles 
near the P.S.D. For w near w the poles at k Q and lie near 

the real axis. Thus, ignoring possible interference of poles at kj's 
caused by roots of e_ (k^to.a) = 0 other than the least-damped ones, 
we conclude that use of the method of steepest descent is justified. 


The a) integration 

We now perform the u> integration in (2.50), which from (2.57) 
and (2.61) is 
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Result of branch-cut 
integral, i.e. (2.61) 




In Fig. 2a the primitive w contour lies above all singularities of the 
integrand of (2.63). We let 6 -* 0 deforming the contour around any 
poles. In Fig. 2b we have lowered the contour until 6=0 and 
deformed it around the pole at u = -ow^ (this is the only pole shown 
in the figure). In Fig. 2c we have lowered the contour into the lower- 
half plane. Since for t > 0 the contribution from the horizontal part 
of the contour vanishes (e -v e I Im q as | Im u> | or t ■+• <® ) , 

we are left with just a residue term at u> = -ou^ . Before writing down 
this term we consider contributions from poles other than the one at 

0) = -OU)j . 

The remaining poles of the integrand of (2.63) occur in the u) 
plane where 

(i) (2.61) has poles 
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(ii) 

(iii) 

(iv) 



In general none of these poles will occur at w = so none of them 

will interfere with the response at frequency . However, we shall 
investigate all four of the above cases to determine if there is any 
response at frequencies other than . 


CASE (i) : Consider just the first (electron) term of (2.51) which 

has poles where any of the following are satisfied: 

= o 

CO - k st V ? +■ <ril = O (2.64) 

e r ± C^ SG ,^ = ° 

2/3 

[Note from (2.59) that k 'v (w+afi) so none of (2.64) would pro- 

s © 

duce simple poles.] By the following argument we show that roots of 
(2.64) in the i o plane cannot produce any legitimate pole terms. 

First note that the position of the P.S.D. for the branch-cut integral 
in the k plane depended on the value of u) [see (2.59)]. Now any 
a) that solves (2.64) represents a pole in the k plane exactly on the 
P.S.D. that corresponds to that u . Recall that in obtaining (2.61) 
we assumed the integrand of (2.58) contained no poles on or near the 
P.S.D. Thus if we could find an w that solved (2.64) we would have 
to evaluate the branch-cut integral (2.58) by a different method (only 
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for that u) ) say by integrating along the real axis in the k plane. 
But then the factors in (2.64) would never even occur in the calculation 
and thus no pole term would result from the w that solved (2.64). We 
conclude that any roots of (2.64) in the u plane can be ignored and 
that (2.61) is a valid representation of the branch-cut integral only 
for w evaluated at w = - 0 u^ . 


CASES (ii) AND (iii) : For definiteness we set a = -1 and 

consider electron waves so fl = -oi . Then the residue term that 

06 

results from case (ii) is 


where 



(2.65) 



(^M) 


= O 


(X = SG-/V 




a , 

Cu - C<j 

ce> 


( 2 . 66 ) 


and the residue term that results from case (iii) is 






where because of (2.55) we must require 


(2.67) 
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( 2 . 68 ) 


€ ;r (/' 3 /fv z +-^ e ) = 0 p = SGrA/{t ^a. k) 


/V 

By making a change of variable in (2.65) from to to kv + to , i.e., 

z cc 


/v 

CO 


k 


+ ^ce 


v. 


2 ^Cu 


Cu-tu 




k=k 


(2.69) 


we see that the two residue terms, (2.65) and (2.67), cancel identically 

with the following considerations. 

For reference, the roots of £. r [ (to-to /v ) ,oo] =0 in the complex 

1 ^ cc z 

to plane are shown for a typical case in Fig. C-17b in Appendix C. Note 
there are branches which have Im to > 0 (corresponding to growing 
waves) which would be included in the to integration. The correspond- 
ing roots of e T (k,kv z +w ce ) =0 in the complex k plane are shown 
in Fig. C-17a. Note there are brandies which have Im k < 0 (corre- 
sponding to spatially growing waves) but it seems these branches would 
be excluded in the k integration since the primitive k contour runs 
along the real axis. In fact, however, if we review the lowering of 
the <o contour, and follow the pole at k = to-to /v [the equivalent 

"k" in e T (k,kv +ui ) ] , we would see that since initially Im to > 0 
I z cc 

+ 

(and v > 0 always for z > 0 ) , as Im to goes to zero and then 
z 

becomes negative, the pole approaches the k contour from above and 
then indents it so it would be included in the final result even though 
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it ends up in the lower-half k plane. (If v < 0 the pole would 
start beneath the primitive k contour and would therefore never make 
a contribution.) 

We conclude that the pole terms produced by cases (ii) and (iii) 
cancel identically. 


CASE (iv) : This case is not allowed in that we assumed a simple 

first-order pole occurred at k^ in performing the k integration 
[see (2.54)]. Even if a simple pole of higher order occurred in the 
k integration, the term corresponding to case (iv) would be the first 
non- zero derivative and would therefore not produce a pole in the 
integrand of (2.63). 

Thus in (2.63) the only pole we have to consider is the one at 
a) = -Ot^ , which gives 



-•[ 


(2,61) evaluated at 


to = -oioj 


B ! ( v> 7 e 


-crf^+JQ) 


j* +i 


0~£t/ t 


, cr) 

‘M V ’ / 


(2.70) 


+ 


e 




* 




for z > 0 and where (k^ ,o) = 0 . Equation (2.70) is valid 

-aB 


for a = ± 1 (not just a = -1 as used mainly in the preceding 
discussion) . 
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Final expression for fjCz.Vjt) 


Using the symmetry properties (B.6) and (B.8) from Appendix B, 


r “] a||? ^ 

\€ r (-cr^-<r<us,cr)\ - ^ (" cr k , ~ , er) 

L ± - 1 + 


3 /(r 


1 * 


3 (-/r *) 


(2.71) 


and the symmetry property 

[£-/£)]*= ^(-X-*) (2.72) 

[which follows from the fact that E^(z) is real], and then summing 
over a = ± 1 [which has been implied since (2.28)], we find f^(z,v,t) 
equals twice the real part of either the a = +1 term or the o = -1 
term. We shall use the latter. Then specializing to electron waves 
(£2 = -od ) , and dropping the ion term in (2.61), we obtain 

Cv 








{Cu-Cc/ ce ) IX £ 


2 a. 


/ . V? 

2 -2«T 


-l^T 
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i k, Z -i<v,tr 


r 

\ »«,(*,*>) 

l v z J 

' ik 


k-k, 


for 2 > 0 and where 


->/ — 


^~s&~ 


r- + *0 

-i?pe ZT 0 


syr> V = 




2{C^-^e)^y 


/ 


? <z‘ 


3 

e 


(2.74) 


oc ~ ?ca/ ^ce) 

(*,."/) * ° 

"f 

Equation (2.73) gives the first-order perturbation to the electron's 
zero-order distribution function. Note that although f^z.v.t) in 
(2.73) is for electrons, both ions and electrons enter into the 
dielectric functions [see (2.44)]. The first term of (2.73) will be 
called the branch-cut term (since it came from the branch-cut integral), 
the second term will be called the free -streaming term (for reasons to 
be discussed shortly) , and the last term will be called the dielectric- 
pole term (since it comes from the root(s) of the transverse dielectric 
function) . 
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2.4 Calculation of J^(z,t) 


The first-order current response is 


J V 


(2.75) 


From earlier considerations [(2.33), etc.] the integration over <p and 
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(2.76) 


for z > 0 . Note that only particles with v^ > 0 contribute to the 
free-streaming term but that all particles contribute to the branch-cut 
and dielectric pole terms. 

To facilitate computation of the free-streaming term integral in 
(2.76) and because at present we have no other preferences, we assume 
an idealized delta- function source with 




(2.77) 
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so 


E, (*) ‘ §, 


(2.78) 


where is a constant with units of electric potential. Recall that 

f (v ) was assumed to be Maxwellian in (2.35). Then the integrals in 
o v z 

the branch-cut term and the dielectric-pole term in (2.76) are of the 
form 



where Z is the plasma dispersion function defined earlier in (2.45). 


The free-streaming term phase-mixing integral 

The free-streaming term integral in (2.76) is 


<^v z € «■ e 
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(2.80) 


and we shall perform this integral in the complex v z /a plane by the 
method of steepest descent described in Appendix E. The appropriate 
saddle point of the integrand of (2.80) is 


Vs 
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(2.81) 


and the first term of the asymptotic expansion about this saddle point 
is 
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(2.82) 


which is a good approximation for 

^ . !Z (2.83) 

A complete account of the paths of steepest descent associated with 
(2.80) and an error discussion concluding with (2.83) are given in 
Appendix E. 

The primitive contour of (2.80) in the complex v z /a plane is 
shown in Fig. 5a. In Figs. 5b, c, d, the primitive contour has been 
deformed to the appropriate P.S.D. depending on the value of . The 
saddle point v s / a [defined by (2.81)] is shown as well as the 
dielectric pole at 


CL 


V£_ - 

~ kcx. 


where 



(2.84) 


Note that there is a pole contribution for to. < 0 or to, > to but 

r l 1 ce 

not for 0 < to. < to . Thus the free-streaming term (2.80) has the 

■ ■ A CC 
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FIGURE 5. CONTOURS IN THE vja PLANE USED IN EVALUATING (2.80). 

The primitive contour in (a) is deformed to the appropriate 
path of steepest descent in (b),(c),(d) according to the 
value of . The saddle points (v £ /a) and poles (Vj/a) 

are defined by (2. 81) , (2. 84) respectively. 
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saddle point contribution (2.82) for all values of , and a residue 

contribution only for a), < 0 or w, > oj . These results, including 

1 1 ce 

the absence of a residue term for 0 < uj, < a) , will be discussed 

i ce 

further in a special section on transverse free-streaming waves. 

Final expression for J^(z,t) 

Combining the results of the above saddle point integration with 
(2.79) and (2.76) we arrive at our final result for the current 
JjCz.t) 

( Z,tr ) = 

— / 



for z > 0 and where 
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o 
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O < cu 

< ^ce 



In (2.85) the terms remain in their original order, i.e., branch-cut 
term, free-streaming term, and then dielectric-pole term except that 
the S(w) portion of the last term is the residue contribution from 
the free- streaming term for < 0 or > w ce . 

Before discussing the behavior of the terms in (2.85) in general, 
we shall present the results of a numerical evaluation of these terms 
for a specific case. 
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3. Numerical Evaluation of (2.85) for Response with Electron Cyclotron 
Damping 

We present the results of a numerical evaluation of (2.85) for 

values of the excitation frequency w.. just slightly below w . This 

x cc 

choice of frequencies was made in the interest of studying electron 
cyclotron damping but in addition this choice clearly illustrates the 
various sizes and types of terms that comprise (2.85). 

First we must consider the effects of all the roots of the trans- 
verse dispersion relation. We have drawn the equivalent of Fig. 3d for 
0 < 0 )^ < ii> ce in Fig. 6, in which we display the positions of all roots 
of e T (k,w ) = 0 , not just the least-damped ones. (The roots of the 
transverse dispersion relation are discussed in Appendix C and are 
summarized in Fig. C-15.) Also shown in Fig. 6 is the branch-cut P.S.D. 
contour as chosen in Appendix li (see the discussion concerning Figs. 

E-l - li-3 for details). As used earlier in Figs. 3 and 4, e T (e,j. ) 

- + 

is to be used on the dashed (solid) line contours in Fig. 6. 

We see from Fig. 6 that in addition to the least-damped root 
(marked k^ ) we must consider the effects of the infinite sequence of 
roots and also the pure imaginary root. In the next section we shall 
show that the infinite sum of residues produced by the infinite sequence 
of roots is negligible except at positions extremely close to the place 
of excitation. Accordingly we shall neglect the infinite sequence of 
roots for now. Note that the pure imaginary root in the upper-half k 
plane in Fig. 6 is missed by the deformed contour altogether so it never 
enters the calculation for 0 < < w ce . Thus the only root we have 

to consider is the least-damped (cyclotron-damped) root. 
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FIGURE 6. FINAL CONTOUR IN THE K PLANE FOR 0 < < “ ce • 

Roots of e^(k,u>j) 3 0 are shown by dots and roots of 

e,p (k,u^)=0 are shown by crosses. Also shown are the 

branch- cut P.S.I). contour (which differs from that of 

Fig. 3d because now w < u ) and the contours that 

ce 

encircle the various poles ( e T _ [e T+ ] is to be used 
on the dashed [solid] line contours) for z>0. 
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We comment on some characteristic lengths. The branch-cut and 
free-streaming terms in (2.85) have the characteristic length 


X = 


Co -co 
/ c 


(3.1) 


ce | 


while the pole term has the wavelength 
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For our example we shall find A. ~ c/to . Thus we have chosen to 

X cc 

normalize all distances to the basic length c/u> , i.e., we define 

cc 
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co c e 


For reference, if f = to /2 tt = 2.25 GHz (a typical value), then 

C6 CG 

c/to ~ 2 cm. 
ce 

In evaluating (2.85) we selected values of the plasma parameters 
attainable on the UCLA Q machine, i.e., c/a = 1120 and to /to = 0.4. 

pc cc 

Then writing 


= 

— / 






(3.4) 
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with given by (2.85) we numerically calculated J x vs. z for 

several values of to./w in the range 0.985 < w./w < 0.999 . The 

1 ce — 1 ce- 

res u Its are presented in Figs. 7a-f. The actual values of w./u and 

the corresponding values of the parameter T defined by 

f (3.5) 

are given in the caption of Fig. 7. All of Figs. 7a-f are drawn to the 

same vertical and horizontal scales to aid in comparing the wavelengths 

and amplitudes of the various terms. And in each of these figures the 

length \ q is indicated explicitly while the length A^ (which is 

simply one wavelength of the pole term) is easily discernible. 

We briefly note the following features in Figs. 7a- f. 

Figure 7a: The branch-cut term is negligible, the free-streaming 

term is significant for 0 z ~ 3 A q and the pole term is significant 

from z = 0 or to large values of z since it is only weakly damped. 

Note that A « A. . All of these features hold for smaller values of 
o 1 

W./U) 

1 ce 

Figures 7a-f: As o),/aj._ increases, the branch-cut term grows in 

' 1 C© 

amplitude and penetrates further into the plasma, reaching peak penetra- 
tion at about w./w = 0.996 (T = 4.5) . The dashed branch-cut curves 
computed by setting e^, (^ se> u)j) equal to one in (2.85) show that if 

the dielectric function e~ (k , uj ) were not present, the branch-cut 

i ^ se x 

term's penetration length would continue to grow with A Q approaching 
infinity as w. approached w . The free-streaming term always 

X CG 


r = 


^ce 


-/ 
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FIGURES 7a- f 


In Figures 7a-f we present plots of [defined in (3.4)] vs. 

Z [=Z/(c/o) ce )] for the branch-cut, free-streaming, and cyclotron- 
damped pole terms of (2.85). ( oj /u> = .4 , c/a = 1120 ) All 

figures are drawn to the same vertical and horizontal scales, which 
are fully labelled only in Figure 7a. In Figures 7b-f, 
was set equal to one in (2.85) to produce the dashed branch-cut 
curves, and e T (u. -w /v ,w.) was set equal to one in (2.85) to 
produce the dashed free -streaming curves. Values of u>j/u) ce > 

r*|(a) 1 /u ce )-l|c/a , VV (c/w ce )sl/r » 31111 k l c/u) ce are as follows 
(values of kjC/u> ce were obtained from Figure C-5) : 


Figure 

0J. /w 
1 ce 

r 

X 

0 

k. c/o> 

1 ce 

7a 

.985 

16.8 

.06 

3.64 + i0.02 

7b 

.990 

11.2 

.09 

4.28 + iO. 145 

7c 

.994 

6.7 

.15 

5.39 + i0.91 

7d 

.996 

4.5 

.22 

5.76 + il .68 

7e 

.9975 

2.8 

.36 

5.89 + i2 . 32 

7f 

.999 

1.1 

.89 

5.95 + i2.96 
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FIGURE 7d: u. /« = .996 

1 ce 
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diverges at z = 0 . As increases, the free-streaming term 

penetrates more as its amplitude decreases; roughly, maximum penetra- 
tion occurs for w./w ~ 0.994 (T ~ 6) . The dashed free-streaming 

x C6 

curves, computed by setting [(u)j-to ce /v s ) ,w^] equal to one in 
(2.85), show that if the dielectric function e T [(“^-“.g/Vg) .(Jj] were 
not present, the free-streaming term's penetration length would also 
continue to grow with X q , approaching infinity as approached 

a) . The cyclotron-damped pole term shows increased damping and a 
smaller wavelength as increases. 

Figure 7f: The branch-cut and free-streaming terms are significant 

only for 0 < z « X q = 1 while the pole term is significant for 
0 < z ~ 2 the later limit being due to the heavy cyclotron damping. 
Note that X q = X J . For still larger value of Wj/Wgg (i.e., 0.999 < 
U] / w < 1 , or equivalently 1 > T > 0 ) , the branch-cut and free- 
streaming terms' penetration lengths go to zero while X q -► 00 ; at the 
same time Xj decreases (but not to zero) and the pole terms' penetra- 
tion length decreases (but not to zero). 

We conclude that the branch-cut and free-streaming terms attain 
sizable magnitudes and penetration lengths for within a certain 

range of values. More specifically, for the present parameter values, 
we shall find in a following section that collisional damping dominates 
cyclotron damping for w ]/ w ce ~ 0-99 • Thus, for the present parameter 
values, cyclotron damping measurements would be feasible only in the 
range 0.99 5 £ 1 . And we see from Figs. 7a-f that within this 

range, both the branch-cut and free-streaming terms are quite signifi- 
cant. In fact, for 0.996 5 <V w ce 5 °* 998 approximately the first 
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1/3 or 1/4 of the distance over which the pole term (and therefore 
cyclotron damping) could be measured has dominant interference from 
the branch-cut and free-streaming terms. 

With the above specific results at our disposal, we return to 
(2.85) to discuss the general behavior of the three terms that comprise 
that equation. 
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4. Analytical Investigation of (2.85) 

In this section we shall examine the general behavior of the three 
terms in (2.85) for all values of . We shall find that the branch- 
cut and free-streaming terms have definite characteristics in each of 
the regions T » 1 , r = 1 , r « 1 , where, as given earlier, 

T = | (oj, / u> )-l| — . [For reference, note that (1) T » 1 for all 
values of except those very near u> ce , and (2) the T regions 

are symmetrically located about = w ce .] The shielding effect that 
occurs in the branch-cut and free-streaming terms, as noted in Fig. 7, 
will be discussed further. Lastly, in considering the pole term, we 
shall investigate the infinite number of residues that occur (because 
there are an infinite number of roots to the transverse dispersion rela- 
tion) and show that their effects are negligible in general, except at 
positions very near z * 0 (the place of excitation) . 

The branch -cut term in (2, 85) . 

This term has the characteristic length X Q = • 

Observe that 


/ 



Thus at z = 0 (and for any Wj except Wj = u> ce ) , the factor 

(k 2 c 2 - oj, 2 ) 'v k 2 e 2 -*■ 00 and therefore the branch-cut terms' amplitude 
v se 1 se 

-*■ 0 . For z ^ 0 the behavior of this term depends predominantly on 
the value of the parameter T : 

r » 1. In this region X Q is a very small distance and the 
branch-cut term effectively damps out in a distance of a few X Q 
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2/3 

because of the exponential factor which goes as exp[-(z/A o ) ] . 

Even for 0 < z < X the branch-cut term is essentially negligible 
— o 

because, although e, f O^e* 0 ^) ~ 1 > we find that 




cuJ; 
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> r 

o __ 


« / 


(4.2) 


2 2 

(for a) /u> not too large) . 
v pe ce 

is negligible in general. 


Thus in this region the branch- cut term 


T z l. As approaches u) ce the branch-cut terms' amplitude 

increases (because of the l/jlo.-a) } factor) and the term penetrates 
farther into the plasma (because X Q increases). If there were no 
dielectric functions in the denominator of the branch-cut term, then 


both the amplitude and penetration length of this term would go to 

infinity at u. = w . However, the dielectric functions are present 
1 ce 

and their shielding effect is strong enough (for r ~ 1 ) to damp out 
the branch-cut term a significant distance before the exponential factor 
does, i.e., before z = A q • (This effect was noted earlier in Fig. 7.) 
A result of this effect is that the branch -cut term attains a maximum 


penetration distance for some value of T larger than one. [In Fig. 7 
2 2 

(for w e / w ce = °*4, c/a = 112 °) maximum penetration occurred at 
T = 4.5 .] 

T << 1. In this region the shielding effect described above domi- 
nates, which may be seen as follows. At z = A Q we find | k se l ~ 1/A 0 
and 
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(4.3) 


€ r ± ^se^ 
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Thus at z = X Q , the branch-cut term goes as 
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(4.4) 


"? 2 

which, for w /w not too small, is « 1 and therefore negligible, 
pe ce 

In this case, the branch-cut term attains a sizable magnitude only for 
1 z I << X . However, for 1 z 1 « X the branch-cut term is not truly 
valid since then the saddle point asymptotic expansion breaks down (as 
shown in Appendix E) . Thus, although the peak behavior near z = 0 of 
the branch-cut term in Figs. 7e,f can be explained analytically, we 
shall refrain from doing so because this term is not truly valid there 
anyway. 

The above results hold for (u) 2 /w 2 )(c/a) >> 1 . If on the other 

pe ce 

hadn (oj 2 /w 2 ) (c/a) « 1 , then (4.4) indicates that the branch-cut 
pe ce 

term would have a sizable magnitude at z = X Q and therefore a sizable 
penetration length (since -*■ 00 for T << 1 ) . In this case the 
branch-cut term's penetration length would simply approach infinity as 
w. approached w . Thus effectively, in this case of extremely low 
densities and/or extremely high magnetic fields, the plasma does not 
exhibit the collective shielding effect that occurred above. 
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Note that 


The free -streaming term in (2. 85). 

This term also has a characteristic length X q . 

Vs 

CL 

A significant feature of the free -streaming term is that it always 

diverges at z = 0 due to the factor l/v s • Physical reasons for the 

presence of this divergence and conditions under which it is removed 

will be discussed in a following section on transverse free -streaming 

waves. For z t 0 , the behavior of the free-streaming term depends 

predominantly on the value of T : 

r » 1. In this region X Q is a very small distance and the free- 

streaming term effectively damps out in a distance of a few . The 

dielectric function e, r [ (a)j-u ce /v s ) , 0 )^] ~ 1 so the free-streaming 

a 1 c '" 

term is approximately 



which depends only on the ratio z/X q . It follows that the amplitudes 
of this term's peaks remain approximately constant but the peaks shift 
to larger values of z as approaches w cg (because then X q 

increases). Thus we find in this region that the free-streaming term 
is always significant for 0 z ~ 3 X q say. 

T = 1. As oo, approaches <o the free-streaming term penetrates 
further into the plasma but at the same time, its amplitude decreases 
as the dielectric function begins to exhibit a shielding effect. This 



61 



effect dominates as T ■+■ 0 and the effective penetration length 
recedes back to zero. From Fig. 7 we note that an "effective penetra- 
tion length" is not a clearly definable quantity but that for some 
value of T larger than one (e.g., T = 6 in Fig. 7), the free- 
streaming term does in some sense exhibit maximum penetration. 

T « 1. In this region the shielding effect described above domi- 
nates, which may be seen as follows. For |z| ~ X q , l v s / a l ~ |z/X Q | 

5 1 and we find 



Thus for | z | ~ X q the free-streaming term goes as 


/ 

V 2 3 



f(M, - OJ ce \ 

~~ £. ( 4 * * - 7 ) 
l X. 0 / Ui£ 

ra 

v s M) 


2 2 

which at z = X is « 1 provided w /w is not too small. In this 
o r pe ce 

case, the free-streaming term, like the branch-cut term, attains a 


sizable magnitude only for |z| « X q . [Equations (4.6), (4.7) for 
z = X q may be compared to the corresponding branch-cut term Eqs. 
(4.3), (4.4).] 


If (<£/%) (c/a) « 1 


then according to (4.7) the free-streaming 


term, like the branch-cut term, would have a sizable magnitude at 


z = X Q (thereby exhibiting no shielding effect) and accordingly its 
penetration length would approach infinity as approached u> ce . 
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The dielectric-pole term in (2.85). 

The pole term(s) represent the collective modes of the plasma. As 
noted earlier, a pole term has wavelength = 2ir/Re k^ = 

[2ir/(Re k^c/cOj^g) ] ( c /to ce ) which, for the least-damped pole, is of order 
c/a) e for near u) ce . The amplitude of the pole term depends 

mainly on the factor Z (E ) where E = (to 1 -oj )/k.a . The quantity 
|z + (E e )| is largest (of order one) for near u) ce ; elsewhere 

|Z + (C e )| 'v 1/E e % a/c which is typically quite small. 

The waveforms produced by the pole terms are governed by the zeros 
of the transverse dielectric function, i.e., by the roots of e T (k,u>) = 
0 for real w and complex k . These roots are discussed in detail 
in Appendix C which should be read at this time. A summary of the 
types of roots for the various frequency ranges is given just preceding 
Fig. C-14. With the results of Appendix C, the pole terms may be com- 
puted for any desired value of . 

The least-damped pole and any pure imaginary poles may be readily 
computed. However, the infinite sequence of poles that always occurs 
is more difficult to evaluate. In the following we obtain an upper- 
bound estimate for the magnitude of the infinite sum of pole terms. 

The infin it e sum of pole terms . 

For 0 < w, <03 the branch- cut P.S.D. contour in the k plane 
1 ce 

(as shown in Fig. 6) is such that only a portion of the infinite 
sequence of roots [of e T (k,oj) = 0 ] is enclosed. This situation 
occurs whenever w. < u> as may be seen from Fig. C-13. For 

X CC 

m > o) the branch-cut P.S.D. contour (as shown in Fig. 3) is again 
1 ce 

such that only a portion of the infinite sequence of roots [now of 
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e^, (k, 0 )j) = 0 ] is enclosed. In all cases the saddle point k ge , as 
given in (2.85), moves as or z is varied, so the number of roots 

enclosed also varies. Thus we have chosen to evaluate the full infinite 
sum of residues, knowing that in any actual case the net effect will be 
smaller than the estimate we will obtain. Without loss of generality 
we consider oj > u> , and therefore the roots of e_ (k,u)..) = 0 . 

1 C6 * ^ A 

Asymptotic expressions for the values of the roots in the infinite 
sequence are obtained in Appendix C. As given in (C.35) such roots of 
e T (k.ojj) = 0 in the plane (C g = [uij-o^l/ka) are given by 
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and at such roots we have from (C.32) 
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We proceed to evaluate the residue terms. 
For | £ e j t 3 we know 
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so from (4.9), at a root, we have 
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For now we assume the second term in (4.11) is negligible and take 
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To simplify this, we note that for 
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(since U J ~ 3 ) , we have 
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(4.15) 


'ZD^Te 


using (4.12). 

We may now calculate the infinite sum of residue terms, which we 
shall call £ . Using (4.8), (4.12), and (4.15), we find 
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where c = (u . -oj )/2/rra . Since 
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for cz i 0 , we conclude that 
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If cz = 0 , (4.16) remains finite since 
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We see that if inequality (4.14) holds (i.e., if T ~ 1 ) than 
(4.12) and the final result (4.18) hold. In the opposite limit 




^ »3 

C2L 


(which corresponds to T » 1 ) we find 
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(4.20) 


(4.21) 


(4.22) 


For cyclotron damping measurements we would usually have u. ~ u 

i. C6 

2 1/2 

= a) . Then at 2 = a/ to ~ A , (X, = [4-nne / kT ] is the Debye 

pe pe d d e 

length), we see from (4.18) that |£| <1 for F 5 1 and from (4.22) 

2 2 “ 

that |£| 5 c /a for T » 1 . Thus very close to the place of exci- 

OO 

tation (z = 0) the infinite sum of pole terms may have a significant 

effect. On the other hand, at z = c/u both (4.18) and (4.22) show 

ce 

that 1^1 < a/c which is typically very small. Thus the effects of 

OO 

the infinite sum of pole terms is confined to a small region near the 
place of excitation. [To substantiate the above conclusions for the 
case of Fig. 7, we summed the first several thousand pole terms numeri- 
cally using (4.11), (4.13) and indeed it was found that the sum was 
significant only extremely close to z = 0 .] 
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5 . Physical Picture of Transverse Free-Streaming Waves and the Free- 


Streaming Term in (2.85) 

We wish to investigate the physical significance of transverse 
free-streaming waves and explain why these waves have negative phase 
velocities (for z > 0 ) when 0 < id. < w , and why these waves 
experience no cyclotron damping. Then we want to investigate the free- 
streaming term in (2.85) and explain why this term has no accompanying 
pole contribution for 0 < w < u> , and why this term diverges at 
z = 0 . In addition to present considerations , free-streaming waves 
merit special attention since they are the basis of echo phenomenon 
(such as we will consider in Part II). 


5 . 1 Physical picture of transverse free-streaming waves. 

We consider the situation shown in Fig. 8a. A continuous uniform- 

density beam of electrons flows in the + z direction along an externally 

applied magnetic field ^ = e z B Q • All particles in the beam have the 

same velocity v in the + z direction and the same magnitude of 
z 

velocity |v x | in the perpendicular direction. The directions of the 

perpendicular velocities are uniformly distributed; tnis is represented 

in Fig. 8b by the uniform ring of dots plotted in perpendicular velocity 

space (v , v ) . Each dot represents a particle and the particles are 
x y 

uniformly spaced on a circle of radius Iv^J . Thus if the perpendicu- 
lar current 



(5.1) 


were calculated, the result would be = 0 as illustrated in the 
lower portion of Fig. 8b. 
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FIGURE 8. CREATION OF A TRANSVERSE FREE-STREAMING WAVE. 

The beam and excitation at z=0 are shown in (a). 
The velocity distribution and current J A are shown 
for z<0, z=0, z>0 in (b),(c),(d) respectively. 
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Returning to Fig. 8a, we have placed an idealized excitation mech- 
anism at z = 0 which lets particles pass through freely but subjects 
them to the external field 

e (^,r) = + e^sv/v (<^0] <^{z) ( 5>2 ) 


[as used earlier in (2.46), (2.77)]. The effect of this external field 
is that as a particle passes z = 0 it receives a net impulsive change 
in its momentum, which we shall call mAv . If a particle passes 
z = 0 at time t we find 


-t + 


it 


/7T?AV — J 
it 




exr 


• J- f. 

\ J 

~&2 


^2 


cxf 


(5.3) 


jf. 


V 2 


e K cos (<^r) +- £ y s-//v (cu f t) 


where fit (6z) is the time (distance) a particle spends subjected to 
the force qE g . (Actually fit, fiz -*■ 0 while Av remains finite.) 
The important point is that all particles that pass z = 0 at the same 
time receive the same Av regardless of the direction of their perpen- 
dicular velocity vectors. 

In Fig. 8c we show the effect on the velocity distribution when 
Av points upward. Each particle receives a Av upward so the circle 
of Fig. 8b is just displaced upward a distance Av . If we calculate 
^ from (5.1) we shall find a net current in the -Y direction (since 
q = -e for electrons) and this is indicated in the lower portion of 
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Fig. 8c. 

Now imagine what happens to the beam segment with the velocity 

distribution of Fig. 8c as this segment moves downstream (z > 0) with 

velocity v . Since every particle (electron) rotates in the right- 
z 

hand direction at the cyclotron frequency u) ce , it follows that the 
velocity vector of every particle maintains a fixed magnitude but 
rotates at frequency to . Thus, as shown in Fig. 8d, the velocity 
distribution is a displaced circle that rotates at frequency w ce . 

The net current J. rotates at cj also, as is shown in the lower 
portion of Fig. 8d. 

In summary we note that a beam segment with \ ^ > 0 has Jj. = 0 
for z < 0 . As the segment passes z = 0 it obtains a net Jx , and 
as it moves downstream with velocity v the current J x rotates in 
the right-hand direction at frequency w . 

Cc 

The shape of the wave produced downstream depends on the relative 
size of uj. and w . If u), = w each new segment that passes 

z = 0 obtains a J x exactly in phase with the 's of all the other 
segments that passed earlier. The established waveform, as shown in 
Fig. 9b, moves downstream with velocity v z and rotates in the right- 
hand direction at frequency u) cg . If u>j t u> ce then the excitation 
and established J x 's rotate at different frequencies, and the resultant 
waveform is a helix. For < w ce the waveform is as shown in Fig. 

9a. (Note that this waveform occurs even if ooj = 0 .) For > w ce 
the waveform curls in the opposite direction as shown in Fig. 9c. We 
emphasize that all waveforms shown in Fig. 9 move downstream with veloc- 
ity v and rotate in the right-hand direction at frequency w 

z cc 
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FIGURE 9. TRANSVERSE FREE-STREAMING WAVEFORMS ACCORDING TO 
VALUE OF . 

All established waveforms propagate with velocity v 
and simultaneously rotate at frequency w ce . 
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We compute the phase velocities of these waves as follows. If 

=0 it is clear that one wavelength equals v z times the time for 

one complete cycle of excitation at frequency tOj , i.e., X = v z (2 tt/ . 

If w ^ 0 we find 
ce 



^ 


(5.4) 


which can be seen by noting that an observer in a frame rotating with 
the electrons (at frequency u> ) sees the excitation as having a 

C6 

frequency w. - u> ce ) . Thus 


X- = 


<LO, 


'ce 


V, 


(5.5) 


If we consider one of the waves in Fig. 9 at some fixed z (> 0) we find 
the motion of the excitation at frequency Wj is exactly duplicated by 
the J, of the wave (but delayed by a time z/v, and rotated through 
an angle of zw /v radians). Thus the phase velocity of one of 

C6 Z 

these waves would be 


^ Phase 





(5.6) 


Note that v , <0 for 0 < u. < w which means the wave looks like 

phase 1 ce 

it is propagating toward z = 0 whereas in fact all the particles com- 
prising the wave are moving away from z = 0 (since v z > 0 ) . This 
apparent paradox may be explained as follows. 

Explan a tion of why transverse fre e -streaming waves h ave negative 
phase velocities (for z > 0) when 0 < a>j < w ce . 

First consider the useful analogy of a rotating barber-shop pole. 
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The pole does not translate but the stripes appear to be moving upward. 

If the pole actually translated slowly downward (while still rotating) 

then the stripes would still appear to be moving upward (but slower) to 

a fixed observer. Thus the stripes would appear to be moving in a 

direction opposite to the direction the pole was moving. This situation 

is an exact analogy to that of transverse free-streaming waves for 

0 < to, < <o 
1 ce 

More quantitatively, consider any waveform in Fig. 9 and choose a 
point on it. In time At the point moves a distance in the +z 

direction and it also rotates through to At radians in the x-y plane. 

Cc 

Therefore to maintain constant phase (i.e., to keep up with a vector 

J x that points in the same direction as a vector drawn to the original 

point) , besides moving ahead v z At we must move an additional Az due 

to the to rotation where 
ce 


AZ _ 
X 


2 rr 


or 


AZ 



^ce \ 


At 


(5.7) 


using (5.4). Thus the phase velocity is 


V 

phase 


V z at +- A? 

Azr 


'2 


T 


A g 
Ar 


from translation from rotation 


(5.8) 


= V. 




/ ce / 


[using (5.7)], which agrees with (5.6). But this calculation explicitly 
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shows that the w rotation can dominate the v translation causing 
ce z 

the phase velocity to be negative for 0 < w. < w . Hence it is 
important to consider transverse free-streaming waves not as fixed 
waveforms that just translate, but as waveforms that translate and 
simultaneously rotate. 

E xplanation of why transverse free- streaming waves experience no 
cyclotron damping. 

Another unique consequence of transverse free-streaming waves is 
that they propagate freely, experiencing no cyclotron damping. This may 
be seen as follows. Consider a beam of particles with velocity v z 
that have established a free-streaming wave with k = u>j-w /v z as in 
(5.5). Then a test particle with velocity V in the z direction sees 
a Doppler-shifted frequency u given by 

Co — U/ f — ~ OJ, ~ 

(5.9) 



We see from this result that u = u only if V = v , which means 

ce z 

that the only particles that could cyclotron damp the wave are the very 
same particles that comprise the wave. Thus there is no cyclotron 
damping effect. It follows that when many free-streaming waves are 
created, each with a different v z , they will each propagate freely, 
not being damped via absorption energy by particles in the other waves. 

5 . 2 Physical picture of the free-streaming term i n (2 . 85) . 

Referring back to (2.76) we see that the saddle point free-streaming 
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contribution to _J^(z,t) in (2.85) is just the result of a superposition 
of many free-streaming waves , each represented by 

exp(i [(w^-o) ce /v z )z-w 1 t]) and weighted by the factor ^ 0 f v 2 )^ v 2 • The 
apparent strong damping of the saddle-point contribution is due to phase 
mixing; all of the free-streaming waves created actually continue to 
propagate downstream (unless destroyed by collisions) but the sum of 
their J A 's (due to the different wavelengths and phases involved) 
quickly adds to zero. 

Explanation of why t h e free-streaming term in (2 . 85) h as an accom- 
panying pole contribution for ojj < 0 or > to ce but not for 0 < < u) cg . 

In (2.76) we note the presence of the dielectric function 
e [(co.-w /v ) ,u) 1 ] which represents the collective effects of the 

1 1 C6 Z I 

a 

plasma. We find that if some free-streaming wave has a v z (real) such 
that (jo.-w /v (real) is nearly equal to the k. (complex) of an 

X CC Z 1 

allowed collective mode [i.e., e^. (k^.w^) = 0 ] then that collective 

a 

mode may be excited according to the following considerations. We con- 
sider z > 0 as above, for which all free-streaming waves have v z > 0. 
Then in evaluating (2.76) by the saddle-point method, we note that the 
v z contour is deformed to that half of the v z plane wherein 
exp [i (w. -u /v )z] represents a damped wave. Accordingly a collective 
mode will be excited, and a pole will result [from the root of 
e T ((od,-w /v ),w,) =0 in the v plane] only if the k.. of the 

A X 06 Z X Z X 

a 

collective mode represents a damped wave also. 

For Wj > w ce all free-streaming waves have positive phase veloci- 
ties. In addition, the roots of (k^,u^) = 0 that have Re k^/wj > 0 
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also have Im > 0 (corresponding to damped waves for z > 0 ). Thus 
it is not surprising that when the contour in the v plane is 
deformed we pick up a residue contribution (see Fig. 5d) . 

For < 0 a similar result occurs (see Fig. 5b). 

For 0 < u < a) all free-streaming waves have negative phase 
velocities. But the roots of e.j. (k^,^) = 0 that have Re k^/w^ < 0 
also have Im k^ <0 (corresponding to growing waves for z > 0 ) , so 
in effect there is no collective mode that has the characteristics of 
a damped free-streaming wave. Thus it is not surprising that when the 
contour in the v^ plane is deformed (for this case) , the pole is on 
the "wrong" side of the contour and hence there is no residue contribu- 
tion (see Fig. 5c). 

lixplana tion of w hy the free-stre aming term i n (2.8 5) diver ges at 
z = 0 and how t his diver ge nce can be rem oved. 

Lastly, in (2.76) we note the presence of the factor 1/v^ . This 
factor is proportional to the amount of time a particle spends traversing 
the excitation region, as may be noted from (5.3). It is this factor 
which causes the divergence of the free-streaming term at z = 0 . 
Physically, this divergence may now be explained as follows. 

Consider an excitation like (5.2) but which, in place of the z 
dependence E^z) = 4>^6(z) has Ejtz) = ct^/A for -A/2 _< z < A/2 
and is zero elsewhere (then the 6-function character is retrieved in the 
limit A -*• 0 ) . Now if the time a particle (with velocity v z ) spends 
in the excitation region (-A/2 £ z £ A/2) is much smaller than one 
period of the excitation, i.e., if 


77 



(5.10) 


% 


then the excitation field will not have enough time to change direction 
as the particle passes. Consequently the particle will receive a net 
impulse in one (transverse) direction as a result of passing through the 
excitation region. As long as (5.10) holds, the smaller v becomes, 
the larger the inpulse the particle receives. In the limit of an infini- 
tesimal but non- zero A we find the v = 0 particles receive a huge 

Zi 

impulse as they pass z = 0 , and accordingly the total response [as 
given by the free-streaming term in (2.85)] diverges at z = 0 . 

On the other hand, if A is really finite, then from (5.10) we 

conclude that the v ~ 0 particles see a perturbing force which is 

z 

spread out over the excitation region and which rotates at frequency 
Wj . Accordingly the v z ~ 0 particles do not acquire a huge impulsive 
change in their 's with the result that at z = 0 the total response 
no longer diverges. Thus a finite-width excitation (A ^ 0) removes the 
divergence. (Two examples of finite-width excitations will be con- 
sidered shortly.) 

It should be clear from this discussion that if collisions signifi- 
cantly affect the v z = 0 particles (as Fokker-Planck collisions do, 
but simple Krook model collisions do not) then the divergence of the 
free-streaming term in (2.85) at z = 0 would also be removed (even 
for A -*■ 0 ) . 
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6. Effects of Collisions on (2.85) 


The Krook Model. 

In the above we have neglected collisions. Now we wish to deter- 
mine for what values of the temperature T and density n collisional 
effects are significant. Specifically we want to know when collisional 
damping dominates and when cyclotron damping dominates. To this end we 
may use the Krook, Bhatnager, Gross model for the collision term in the 
Landau- Bo ltzman equation. Then the Vlasov equation (2.1) is to be 
replaced by 

. - (A£. 

( 6 . 1 ) 

(H\ 

where v is an effective collision frequency. An appropriate choice 
for v is 



). 




V 


X + 






( 6 . 2 ) 


where v , v . , v represent effective collision frequencies for 
66 6i en 

electron-electron, electron-ion and electron-neutral collisions, 
respectively. 

First consider v . An extensive accumulation of experimental 
en 

( 22 ) 

data concerning v en is presented by S. Brown. From these data we 

conclude that for our choice of parameters (T = 0.2 eV, neutral gas 

✓ — 4 

mainly N 2 , neutral gas pressure p Q ~ 10 mm: these being typical 
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values for the UCLA Q machine), V gn may be neglected in the following. 
Now consider v and v . . In the literature there exists some 

66 6 X 

confusion as to what v , v should represent and whether v or 

66 61 vC 

v . or both should be used in (6.2). Often v in (6.1) is called an 
ei 

effective collision frequency for momentum transfer without any explicit 
indication of whether the momentum transfer occurs in electron- electron 
or electron-ion collisions. In other cases it is explicitly stated 
that v represents a self-collision time (electron-electron collisions 
only).^ 12 ^ In previous work on the collisional damping of transverse 

r i 3i . . . 

plasma waves, Crawford et al. claim the dominant collision process 

is electron-ion collisions and hence use v = V gi . Sachs presents 

a collision term (developed by Liboff) which for T g = T\ and no net 

drift velocity between species reduces to v * where V ge is 

Spitzer's^ 23 ^ self-collision frequency. To further confound the situa- 

tion we note that V ~ V . according to Kaufman et al. Perhaps 

66 61 

the easiest (but least satisfying) solution is that of Platzman and 

Buchsbaun/ 12 -* who avoided the issue altogether by simply choosing a 

number for v . Tanenbaum concludes that the Krook collision model is 

"really applicable only to cases where self-collisions are of predomi- 
(25) 

nant importance." 

Since our interest lies in determining only when collisions are 

significant we have, without resolving the above conflict, chosen to 

use v = v (self-collision frequency) with the understanding that 
66 

this v may differ from the true effective collision frequency by at 

/ i nr \ 

most a factor of 2. Thus we shall use^ ’ 
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_3 

where n is the electron density in cm and T is the electron 
e e 

temperature in eV. 

Then reviewing the derivation of (2.85) when collisions are 
included by (6.1) we find that (2.85) should be replaced by 



for z > 0 and where v is given by (6.3). Now 


81 



^ 5 e = 


f z{( ^s+iV-^e)*} 2 
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rr 
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Z <*? 


Ji 
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C^+iV -CU cq]^ Z 
"2 £2- 




(6.5) 


and 


<E r (& } aj) - /- 


, a 

/°C 


Uv- 


&?)*. 



Co V~Co t 


A:c 


-) 


( 6 . 6 ) 


whereas a and S^) are as given earlier. 

The main effect of collisions in the transverse excitation response 
given by (6.4) will be in the pole term; the cyclotron-damped root of 
e^, (k^) = 0 will be altered due to the presence of collisions. The 
branch-cut and free-streaming terms in (6.4) will also be affected but 
we defer discussion of this until later. For now we consider the effects 
of collisions on the cyclotron- damped root of e^, (k,a>) = 0 . 

Cyclotron damping vs, collisional damping. 

The transverse dispersion relation derived earlier for the colli- 
sionless case, as given in (2.44), is (neglecting ions) 


€ (*r,Ca) = / ~ 

'■+ 


/><? 


7L' (**-*- 


■Cu*)kt 




(6.7) 


When collisions are included using (6.1), the transverse dispersion 
relation is 


b 
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(k,oj) - / - 


CaCc/* 

i°e 






O 


( 6 . 8 ) 


where v is given by (6.3). For jc v l = | (o)+iv-w ce )/ka | ~ 3 , we may 
use the first term of the asymptotic expansion of the Z function 
[Z(C) = -1/5] in (6.8) to obtain 

e (k,co) = / + — - O 

+ (**c*-^ z X <5s ' + * y - c ^*) (6.9) 

which is simply the cold plasma dispersion relation with collisions 
included. 

To estimate where cyclotron damping is comparable to collisional 
damping, Crawford et al/ 13 "* computed the locus of points where Im k 
of the root of the hot plasma dispersion relation without collisions 
(6.7) equals Im k of the root of the cold plasma dispersion relation 
with collisions (6.9). We have numerically performed the same calcula- 
tion for two sets of parameters with the results shown in Fig. 10; 
curve 1 corresponds to Crawford's parameters (and his curves), whereas 
curve 2 was computed for our parameters . These loci show that for a 
given a) , as T g increases, collisional damping decreases and cyclo- 
tron damping increases. However, these loci contain a major flaw and 
that is the following. Use of the cold plasma dispersion relation (6.7) 
as an approximation for the hot plasma dispersion relation (6.9) is 
justified only for |c v l ~ 3 . But along the hatched portions of the 
curves in Fig. 10, ~ 2 so (6.9) should not be used. Indeed, as 
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FIGURE 10. LOCUS OF POINTS WHERE HOT PLASMA CYCLOTRON DAMPING 
EQUALS COLD PLASMA COLLISIONAL DAMPING. 

13 

Locus of points, following Crawford et. al. , where Im k of 

the least-damped root of (6.7) equals Im k of the root of 

(6.9). CURVE 1: f =2.25 GHz, w /u =4 (n =10 12 cm" 3 ). 

ce pc ce c io 3 

CURVE 2: f =2.25 GHz, u /w =.4 (n =10 cm” ). 
ce ’ pe ce e 

, (Curve 1 corresponds to Crawford's parameters.) 

These loci are meaningless however for those portions 

shown hatched (see text). 
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we shall now show, the true loci never turn upward (as w/^ ce -*■ 1 ) as 
the loci in Fig. 10 do. 

To correctly compare cyclotron damping with collisional damping, it 
is apparent we should compare the roots of (6.7) with the roots of (6.8). 
Let 

L = Im k of the root of (6.7) 

M = Im k of the root of (6.8) 

then L represents the effects of cyclotron damping whereas M repre- 
sents the effects of cyclotron damping plus collisional damping. Thus 
roughly, L/M represents the fraction of the total damping due to 
cyclotron damping. Accordingly we have numerically calculated the locus 
of points where cyclotron damping accounts for 80%, 50%, 20% of the 
total damping (L/M = 0.8, 0.5, 0.2) with the results displayed in Fig. 

11 as curves A, B, C, respectively. The curves of Fig. 11 were computed 
for the same parameters used to compute curve 2 in Fig. 10; in particu- 
lar note that curve B in Fig. 11 does not turn upward as curve 2 in 

Fig. 10 does. Also note that to measure cyclotron damping we must be 

above the shaded region in Fig. 11 for collisional effects to be negli- 
gible. For T g = 0.2 eV this means we should consider only frequencies 
in the range 0.988 ~ co/u> < 1 . 

In the shaded region for T g ~ 0.1 eV in Fig. 11, the total 
damping is very small. Thus one should not interpret curve A of Fig. 11 
to mean that collisional damping is becoming as large as cyclotron 
damping, but rather that cyclotron damping is becoming as small as 
collisional damping. This conclusion is apparent when one considers 
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FIGURE 11. CYCLOTRON DAMPING AND COLLISIONAL DAMPING COMPARED 
CORRECTLY. 

A. 20% collisional damping and 80% cyclotron damping. 

B. 50% collisional damping and 50% cyclotron damping. 

C. 80% collisional damping and 20% cyclotron damping. 

10 _ "K 

( f =2.25 GHz, u /w =.4 [n =10 cm ] ) 

ce * pe ce 1 e 1 
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actual values of Im k of the roots of (6.7), (6.8) such as those 
displayed in Fig. 12 for T = 0.204 eV. Note that in Fig. 12 colli- 
sional damping is significant only when cyclotron damping is very small 
Also note that a typical cyclotron damping length (for negligible colli 

sional damping) is (for w i/ w ce = 0.99) 

<c 

/ 2 Cnrn 

— _ » = /O C/w. 

Jfrn K Jk* k c ■ Z ' 

and damping lengths of this order should be easily observable experi- 

10 -3 

mentally. Thus for the parameters of Fig. 11 (n g =10 cm , 

“pe^ce = °' 4 ^’ experimentally, measurement of electron cyclotron damp- 
ing rates would be feasible only for T g ~ 0.1 eV. 

As the density n increases, collisional effects rapdily become 
© 

more significant since v ^ n g . For example consider values of Im k 
12 -3 

for n g = 10 cm as given in Fig. 13. We see that for 

0.988 < oj./w < 1 , 

— 1 ce — 


Im k 
Im k 


no collisions 


> 0.8 


with collisions 


but that at w./w = 0.998 , Im kc/u ~ 8 which implies a damping 
X c© c© 

length of 1/Im k = 0.27 cm, which is too small for accurate experi- 

10 -3 

mental measurement. Thus for "high" densities (n g >10 cm ) we 
must go to "high" temperatures (T » 0.2 eV) if we wish to obtain 
a range of frequencies over which the cyclotron damping length is 
unhampered by collisional damping and also large enough to be easily 
measurable experimentally. 
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FIGURE 13. 


EFFECTS OF COLLISIONS ON THE CYCLOTRON- DAMPED ROOT. 



Effects of collisions on the branch-cut and free-streaming terms. 


In the regions where cyclotron damping rate measurements are 
feasible, the collisional damping rate is small; i.e., from the above 
considerations we find 


V 

< -°ce 



( 6 . 10 ) 


Accordingly, in these regions, the free-streaming perturbations will 
eventually damp away by collisional damping (these perturbations are not 
affected by cyclotron damping) but phase mixing of these perturbations 
[as given by the free-streaming term in (2.85)] will occur long before 
the individual perturbations themselves are destroyed by collisions. 
Explicitly, when (6.10) holds, the free-streaming term in (6.4) differs 
only slightly from the corresponding term in (2.85); the analysis of 
Section 4 holds and the free-streaming term still diverges at z = 0 . 
(As mentioned earlier in Section 5, Krook model collisions do not 
remove this divergence, whereas Fokker-Planck collisions do, as will be 
shown in the echo case of Part II.) 

If (6.10) does not hold, i.e., if 


V 


then the analysis of Section 4 
which is now 


-§T » / C6.ll) 

does not hold because the parameter T 


r = 


^ + ir. / 

OCf US 

ce ce 


c 

CL 


( 6 . 12 ) 


is never ~ 1 (even for 


w. 


Cl) 


ce 


)• 


Thus the shielding effect described 
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earlier for T ~ 1 never occurs in this case. On the other hand, the 
characteristic length X q which is now 

<z 

° j CtJ+iy ~^ e | (6.13) 

no longer approaches infinity as oj. approaches tu , and the exponen- 

X cc 

tial factor in the free-streaming term still damps this term away in a 
distance of a few X 

o 

Similar remarks apply to the branch-cut term. If (6.10) holds 
then the branch-cut term in (6.4) differs only slightly from the cor- 
responding term in (2.85) and the analysis of Section 4 remains valid. 

If (6.11) nolds then the branch-cut term (like the free-streaming term) 
exhibits no shielding effect as 0 )^ approaches io ce but it still damps 
away in a distance of a few X q where X q is given by (6.13). 
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7. 


Effects of a Finite-Width Excitation Mechanism on (2.85) 


Up to now we have used the shape factor 

£,( 2 ) = (7.1) 

in the transverse external exciting field 

£ exr ( H,*r; = £)( 2 ){e x cos + e^s//v o/r] (7.2) 

(This form of excitation for transverse waves is the mathematical equiv- 
alent of the dipole-grid excitation mechanism for longitudinal waves.) 
The choice of a delta function in (7.1) is useful for computational 
purposes but unrealistic in that it is probably experimentally unattain- 
able. Indeed it is difficult to conceive of a mechanism (or even a 
limiting case thereof) that would produce a field like (7.1), (7.2). On 
the other hand, these are conceivable mechanisms that would produce a 
field like (7.2) provided E^(z) had a finite z dependence. In this 
section we shall examine the effects of a square-shaped E^(z) and a 
Gaussian-shaped E^z). 

First note that in Section 2 we have already carried the general 
calculation for arbitrary E^(z) up to (2.76). Thus we only have to 
alter the derivation from (2.76) to the final result (2.85). Just 
glancing at (2.76), (2.78), (2.85) indicates that inclusion of a new 
E^(z) will simply produce new amplitude factors for the branch-cut and 
dielectric-pole terms while the free-streaming term will be affected in 
a more significant way. 
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Square-shaped E^z). 

Consider the square-shaped E^(z) 


d> 

EX*) ~ ~~ for 

' A 


- A «c 2 « 

Z * 2 


O for 


?<-4 , 


(A > 0) for which 


£(*) = 

' ArA 


(7.3) 


(7.4) 


Using (7.4), the pole term in (2.85) acquires an additional ampli- 
tude factor of sin(k 1 A)/k 1 A . If A « X^, (Xj = 2ir/k ^) then 
EjCkj) ~ 4^ and the pole term is essentially unaltered. However, as 
A increases, the factor sinfk^/kjA oscillates and its magnitude 
decreases. If kjA = niv (n = 1,2,...) , i.e., if A = n(Xj/2) then 
Hi Ckj) = 0 identically. Therefore if we want to avoid these zeros and 
maximize E^(k^) we should have 

A« — (7.5) 

Z 


Since Xj is typically a few centimeters, condition (7.5) should be 
easily attainable. 

Both the branch-cut and free -streaming terms in (2.85) have 

(roughly) the characteristic wavenumber k = 2 tt/X (X = a/lu, -oj I) , 

o o o 1 1 ce ' 

so from the above considerations we conclude that both of these terms 
will acquire additional amplitude factors of sin(k Q A)/k o A . Thus if 


A « ° 


2 


(7.6) 
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we shall have E^(k Q ) ~ ^ and the amplitudes of both the branch-cut 
and free-streaming terms will be maximized. 

In general, condition (7.6) is harder to attain than (7.5) because 
in general \ Q « X 1 . However, in studying cyclotron damping our 
interest lies in the pole term and in fact it would be to our advantage 
to have the branch-cut and free-streaming terms smaller than as given 
by (2.85). In that case it would be useful to have A such that 

o < ~ ( 7 -7) 


because this would help diminish the size of the "unwanted" branch-cut 
and free-streaming terms. 

If we are explicitly concerned with the detailed behavior of the 
free-streaming term we must evaluate the free-streaming term integral in 
(2.76) more carefully, taking into full account the explicit v z 
dependence of E 1 [(u).-ui )/v ] . Physically the integral in (2.76) is 

more complicated now because the particle trajectories are altered con- 
tinuously over the region - j < z < j rather than abruptly at z = 0 . 
Using (7.4) for E^k) the free-streaming term integral in (2.76) 
becomes 


/ 

iffa. 


CO 



o 





l* ? -e * I 


2 J 
* \ ? 


(7.8) 


Evaluating (7.8) by the method of steepest descent we find the free- 
streaming term contained within the large brackets in (2.85) should 
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be replaced by 


-(Cc/r 


2 


(OJ f -OJ ce )A 


(c^~C^. e )oc(ZfA) 


■Z a. 




y . <?' 
2 +<:o< 2 






\c^ r os ce ) A) 

z /s 

t , . 1§' 
— -+yo< — 

•2 <3. 


2 Z 




(7.9) 


where 


V< 


S-H 


(U J/ -Uy <e )u(Z±^> 


"2. at. 


rr 

-I 0 C 3 


(7.10) 


When A -*■ 0 (7.9) reduces to the free-streaming term in (2.85). 

But for A ^ 0 note that (7.7) does not diverge at z = 0 (since when 

v -*■ 0 we have e„, [(w.-w /v ),u 1 l ~ 1 . The 1/v factor in 
s^ 1 1 1 ce s . 1 z 

± a ± 

(2.76) which caused the z = 0 divergence earlier, was effectively 
removed in (7.8). Thus, as stated earlier, a finite-width excitation 
mechanism removes the z = 0 divergence of the free-streaming term in 


(2.85) 


Gaussian-shaped E^fz). 

In place of (7.3) we now consider a Gaussian-shaped E^(z) 




2_ 

A 2 


i7r a 


(7.11) 
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(A > 0) for which 


2 

At a 

E,(At) = % e V (7.12) 

The field described by (7.11) is more realistic than that of (7.3) in 

that any excitation mechanism will be accompanied by a shielding sheath 

with an exponential dependence like that of (7.11). The consequences 

of (7.11) are very similar to the consequences of (7.3). 

Thus when E^z) is given by (7.11) we find the pole term in 

2 2 

(2.85) acquires an additional amplitude factor of exp[-k 1 A /4] which 
is £ 1 for k^A^ « 4 , i.e., for 

A«-r~ ^ A, (7.13) 

Similarly the branch-cut and free-streaming terms in (2.85) each acquire 

2 2 

(roughly) an additional amplitude factor of exp[-k Q A /4] which is 
~ 1 for 

Kl * ** 

(7.14) 

Thus as earlier, if 

O c ~\ 0 « A « X / (7.15) 

the sizes of the branch-cut and free-streaming terms will be diminished 
from those in (2.85) whereas the pole term in (2.85) will remain essen- 
tially unaltered. 

Evaluating the free-streaming term more carefully, we find that 
the free-streaming term integral in (2.76) should be replaced by 
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A 2 vj_ 

2 e € 

[Note that this integral, like (7.8), does not diverge at z = 0 .] We 
shall evaluate (7.16) by the method of steepest descent for the two 
limits | z | » A and |zj « A (the latter limit being of interest 
since from it we can obtain the size of the free-streaming term at 
z = 0 ) . 

2 

(i) | z | >> A. In this region we may consider the A exponential 

factor in (7.16) as part of the slowly-varying portion of the integrand. 

Then evaluating (7.16) shows that the free-streaming term in (2.85) 

simply acquires an additional amplitude factor of 
2 2 2 

exp (-[A (w 1 -a) ce ) ]/[4v s ]) where v g remains as given in (2.76). 

(ii) | z | << A. In this region we may consider the factor 
exp[i(w.-a) /v )z] in (7.16) as part of the slowly- varying portion of 
the integrand. Then the free-streaming term contained with the large 
brackets { 1 in (2.85) should be replaced by 




CL 

e 



- Us _ 
/ ce 





Uy /~ Uj cC 
Mr 3 



(7.17) 
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where now 



CL 



( (Us, -co ct )c< 
2 <z. 


(7.18) 


Expression (7.17) gives the amplitude of the free -streaming term at 

2 = 0 . 
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8 . 


Conclusions 


We have examined the linearized response of a Vlasov plasma to the 
steady-state excitation of transverse EM waves by an idealized delta- 
function excitation mechanism, (1.1) - (1.2). We summarize the main 
conclusions : 

(1) The plasma response is given by (2.85) which is valid for any 
value of . (For use in the pole term, the roots of the transverse 
dispersion relation for real w and complex k are discussed in 
Appendix C.) 


(2) For cyclotron-damping measurements, the following points must 
be considered: 

(a) The branch-cut and free-streaming terms in (2.85), in addition 

to the pole term, have significant amplitudes and penetration 

lengths (see Fig. 7). 

2 2 

(b) If (a) /a» ce ) (c/a) » 1 , then the branch-cut and free- 

streaming terms are shielded for T = |(o) 1 /co ) - l| — ~ 1 

1 C6 cl 

and their penetration lengths go to zero as T goes to zero. 

2 2 

If (to /co ) (c/a) << 1 , then no shielding effect occurs and 
pc cc 

both the branch-cut and free -streaming terms' penetration 
lengths go as X Q - a/|wj-u> | which approaches infinity as 
T goes to zero. 

(c) Use of only the "least-damped root" in the pole term was 
justified by showing that the remaining infinite number of 
roots of the transverse dispersion relation has a negligible 
effect except at distances very close to the place of excita- 
tion (z = 0) . 
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(d) Collisional effects dominate the damping except when u)j is 

near w (see Fig. 11). 

ce e 

(e) For u). > to , the cyclotron-damped root is so heavily 

X 06 

damped (see Appendix C) that it is probably not possible to 
measure it experimentally. 

(f) A finite-width excitation mechanism removes the divergence at 
z = 0 of the free-streaming term, and selectively reduces 
the amplitude of the three terms in (2.85). Under special 
conditions [(7.7) or (7.15)] the branch-cut and free-streaming 
terms are reduced significantly while the pole term is left 
essentially unchanged. 

From (d) and (e) we conclude that cyclotron-damping measurements 
should be performed in a narrow range of frequencies just below to 

C6 

(e.g., 0.998 ~ ~ 0*999 for the parameters of Fig. 7). This 

range of frequencies can be deduced from graphs like Fig. 11. It was 

concluded that cyclotron damping measurements are barely feasible for 

10 -3 10 -3 

T = 0.2 eV and n = 10 cm For higher densities (n >10 cm ) 
G 6 6 

one must go to much higher temperatures (T g » 0.2 eV) because the 

3/2 

Coulomb collision frequency, (6.3), goes as v 'v n /T . 

6 6 

(3) A special section on transverse free-streaming waves showed 
how these waves are created and explained many of their interesting 
properties (e.g., negative phase velocities for 0 < < to ce , absence 

of cyclotron damping, etc.). The free-streaming term in (2.85) was 
shown to be the result of the superposition of many such free-streaming 
waves, and explanations were given for some unusual properties of this 
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term (why there is no accompanying pole contribution for 0 < u)j < w ce , 
why this term always diverges at z = 0 , etc.)* Since free-streaming 
waves are responsible for the production of echoes, all of the unique 
features of transverse free-streaming waves discussed above will be 
especially significant when considering transverse plasma wave echoes 
in Part II . 
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II. STEADY-STATE TRANSVERSE PLASMA WAVE ECHOES 

1. Introduction and Physical Picture 

The phenomenon of longitudinal plasma wave echoes, introduced by 

f261 127-291 

Gould and O'Neil 1 } and subsequently observed experimentally, 1 

may be described as follows. If in a hot collisionless plasma a wave 
with wavenumber k 1 is excited at some time which subsequently damps 
away; and if another wave with wavenumber k 2 is excited at a later 
time and it also damps away; then at a still later time an echo with 
wavenumber j k 2 ~k 1 1 will appear and subsequently damp away: this is 

called a temporal echo process. The corresponding process for spatial 
echoes is also readily visualized: if an excitation at frequency 

is made continuously at some point in a hot collisionless plasma and 
the response damps away in distance; and if further away another excita- 
tion at frequency u> 2 is made continuously whose response also damps 
away in distance; then still farther away an echo at frequency l W 2 ‘ W l I 
will appear and subsequently damp away in distance. In this work we 
focus our attention on spatial echoes since they are more readily exper- 
imentally attainable. 

An echo results when the free-streaming perturbations produced by 
a first excitation are appropriately reordered by a second excitation. 
But since these perturbations are relatively sensitive to small angle 
Coulomb collisions, such collisions may destroy the formation of an 
echo. Alternatively though, the collision process can be investigated, 
as well as micro- turbulence, and also details of the particles' zero- 
order velocity distribution function by the appropriate use of echoes. 
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Thus echoes may serve as an important diagnostic tool. 

In this work the echo process described above for longitudinal 
plasma waves is generalized to the case of transverse plasma waves. 

This extension should greatly enlarge the potential importance of echoes 
as a diagnostic tool. Also we shall find transverse echoes offer a 
much richer variety of phenomenon to investigate than purely longitudi- 
nal echoes offer. In addition, a new use for echoes will be suggested 
and that is to excite instabilities in a well-controlled and localized 
manner. 

We first present a physical picture of the echo process. We choose 
to have a transverse excitation with frequency 10 ^ at z = 0 and a 
longitudinal excitation with frequency at z = L as shown in Fig. 

1. Both excitations produce waves with vectors along an applied extern- 
al magnetic field » e^ . Considering just those particles with 
cyclotron frequency S2 and velocity v z , we see that the first excita- 
tion produces a free-streaming wave with wavehumber [(Wj+fl)/v z ] , as 
discussed in Section 5 of Part I. For z ^ L the second excitation 
modulates the first's free-streaming wave to produce a wave with wave- 
number [(o) 2 +u 1 +fi)/v z ] . The net phase $ of the wave for z L 
includes the phase for 0 _< z <_ L and the phase for z ^ L , i.e.. 



f OU z -\-U/, +X2\ ( 

( ; ) L + ( 

, * n 




= [(^ 2 -F^ + n)z 


( 1 . 1 ) 


When we add the effects of many beams with different v 's (to obtain 

z 

the net current or electric field) we find that all their phases are 
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Z^O 

v J 

It echoes 


' s — v 

//echoes 


FIGURE 1. PRODUCTION OF TRANSVURSE PLASMA WAVE ECHOES FROM 
A TRANSVERSE EXCITATION AND A SIMULTANEOUS 
LONGITUDINAL {£) EXCITATION. 


different and therefore the net response phase mixes to zero everywhere 
except where 





o 


i . e . , at 

^y z L 

2 = = 

uj z +oy / + fi ( 1 . 2 ) 

which defines the location of the echo peak. 

The above physical picture is based on free-streaming particles 
that are non-interacting. If interactions are taken into account, then 
the possibility of exciting collective modes exists. The location of 
the echo peak remains unchanged (since it still depends only on the 
free-streaming terms) but the resultant echo shape may be appropriately 
enhanced. We shall study these effects via a full Maxwell-Vlasov 
equation analysis . 

In the following we consider wave propagation exactly along the 
direction of an external magnetic field and we present a concise mathe- 
matical explanation of what type echo is produced, given ,the types of 
each of the two excitations (longitudinal or transverse). We find that 
a transverse echo is produced in lowest order only when one excitation 
is longitudinal and the other is transverse. For such excitations, or 
more specifically for 


= ^(z)|jyros U'T + e y +f 2 (z-0[e ? cos 0-3) 

we then derive the second-order (nonlinear) echo response from the 
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Vlasov-Maxwell equations, up to an integral over the velocity variable. 


We then discuss transverse echo characteristics, summarizing the 
results in Fig. 2. The final integral over v z is evaluated by the 
method of steepest descent and numerical evaluation of the results for 
three specific cases are presented in Figs. 3-5. 

Several important extensions are then considered: (1) effects of 

Fokker-Planck collisions and finite-width excitation mechanisms, 

(2) effects of temperature anisotropy (T x t T z ) , and (3) effects of 
propagation just off-axis (k ± ^ 0) . 

Lastly, transverse echoes for the case of no external magnetic 
field (B^ = 0) are examined and their characteristics are summarized in 
Fig. 6. 

In performing the above calculations, reference will be made to all 
five appendices. As stated in the Introduction of Part I, the appendi- 
ces contain many useful analytical and numerical results. Of special 
significance for the present echo calculation are Appendix C (on the 
roots of the transverse dispersion relation for real u) and complex k) , 
Appendix D (on the roots of the longitudinal dispersion relation for 
real u> and complex k ) , and Appendix E (on the method of steepest 
descent as used in evaluating phase-mixing integrals). 
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2 . 


Derivation of the Second-Order (Nonlinear) Echo Response from the 


Vlasov-Maxwell Equations 


2.1 Formulation of the problem. 

The Vlasov equation for a plasma in an external magnetic field has 
already been discussed in (2.1) - (2.7) of Part I for the perturbation 
expansion 

= ;£(Y) +- (*,¥.,?) +••• 

+ £,(*,*) +••• (2-U 

+ 2? 2 (x,*) +*•• 

In the notation used there, the Vlasov equation is in zero order. 


in first order. 


-Xl 


<3 9$ 


= o 


( 2 . 2 ) 


+ V- — + = ° <2 ' 3) 


2-f 

JZ + -‘ OX ^ ^ 


and in second order, 


a 

st 


— +y- -o 3<8 i * 2 i- c - - 2 j 3 v_ 




3V 


(2.4) 


The nonlinear term (the right-hand side) of the second-order equation, 
(2.4), is the term that will produce an echo. Given two excitations, 
we will solve (2.4) to obtain the echo response. 
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Note that the right-hand side of (2.4) serves as an effective 
driving term (effective external source) for the full equation. It fol- 
lows, as we shall find shortly, that in Fourier- Lap lace transform space. 


« [(£ + ^W) + *(*+£«,**-■)] * (2.5) 

Let _E represent [(Ej + E^^) + x (B^ + ) ] for the transverse 

case or [E^ + E ] for the longitudinal case. Then in general, if 
two spatial excitations are made as shown in Fig. 1, we will have 


E = E (cj^ +■ E (UJ 


( 2 . 6 ) 


and 

^ + i) C<^ z ) 


(2.7) 


so the product E_ • f ^ contains four terms. Each term produces a 
response at the sum or difference of the frequencies involved, as 
follows : 



O, 2^, 

(2.8a) 



(2.8b) 

[s(<^ 5y ] 

^ ± 

(2.8c) 

[fCog-lj-) y>0,) 


(2. 8d) 


The first two terms are second-order corrections to the individual exci- 
tations (at zero frequency and the second harmonic) and these would 
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characteristically damp away in a few Landau damping lengths from the 
position of the individual excitations. The last two terms represent 
the nonlinear interaction between the two excitations and these are what 
drive the echo. Typically the fields (E^'s), and the branch-cut and 
dielectric-pole portions of f^ damp away in a distance 6 « L where- 
as the free-streaming part of fj^ persists indefinitely (if there are 
no collisions). Thus if we are interested in echoes only for z > L , 
we know £(u>p will be essentially zero at z = L , and therefore the 
only term we need consider in calculating the echo is (2.8d). Similarly 
for echoes at z < 0 , we need consider only (2.8c). 


2.2 Possible echo types. 

We now give a concise explanation as to what type of echo will be 
produced, given the types of excitations (transverse or longitudinal) 
that are made. Consider just the <j> dependence of f 2 in (2.5) (the 
angle <}> is defined in Fig. 1 of Part I). For a transverse excitation, 
we know that £ • ^ ^ and fj ^ e -1< ^ [see (2.30) and (2.73) of 

Part I]. Similarly we shall find that for a longitudinal excitation 
both E • and f. have no <}> dependence. The current response is 

— oV 1 


—z 


= 


'f. 

'Zrr 

ft 


2 v cC v 


Vi ~i <f> 


e +- n e 


Vj_ 

VF 




-Z v * 


J ^ 0 


(2.9) 


using the circular polarization representation of Part I [see (2.33) of 
Part I], and where we have indicated only the <(> . integration. If 
has no <t> dependence then the only term of (2.9) that survives the <}) 
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integration is the e term so the resultant will be a longitudinal 

z ^ 

wave. 


!. If e + ^ then the only term of (2.9) that survives the <{> 


integration is the e + term so the resultant J_ 2 will be a transverse 

3 


wave. Recall that ^ a - 


3v 


f^ and that for an echo calculation 


L- 3v 


comes from the second 


f comes from the first excitation and 
excitation. Let t~ represent a transverse excitation and & represent 
a longitudinal excitation. Then for excitations in the sequence listed, 
we will obtain the following echo types: 

//— / 
f£ — > t 

( 2 . 10 ) 

£t — > t 
tt —+ Jt 

These results are for the lowest-order echoes that will occur when all 
waves considered propagate along the direction of the external magnetic 
field . Later we shall consider the effects of having propagation 
just off-axis. 

The strictly longitudinal case (i 7 /-*/) is the original echo case 
investigated by Gould and O'Neil, et al., and subsequently verified 
experimentally . (27-29) Th e production of a longitudinal echo by two 
transverse excitations ( tf-+j f ) has been investigated in theory 
only. We shall be concerned only with the combinations that produce 

transverse echoes, i.e., T and 7*\ These cases have not 

been investigated previously. 


2. 3 Transverse plasma wave echoes . 

We choose to make a /" excitation with frequency o)^ at z = 0 
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and an excitation with frequency at z = L as shown in Fig. 1. 

Explicitly we choose 


= ^(2)[e x co S (^r)-f- e^s/A/C^r;] t £ 2 ( ? -L)[€ f cos(^t)] ( 2 .n) 


where the E^(z) term is the same transverse excitation used in Part I, 
and the E 2 (z-L) term is the longitudinal excitation at z = L . By 
considering echoes at both z > L and z < 0 we can investigate both 
the 'fj? case (for which the echo is at z > L ) and the j£~f~ case (for 
which the echo is at z < 0 ). For the transverse excitation we already 
have an expression for fj(z,v z> t) [i.e., (2.73) of Part I] and we 
can easily obtain an expression for £(z,t) [by taking the inverse 
Fourier-Laplace transform of (2.48) of Part I], We shall need analogous 
expressions for the longitudinal excitation and we now give a brief 
summary of them. 

In first order the Vlasov equation (2.3) is linear in IE, B_, and f^ 
and therefore the effects of the longitudinal excitation may be con- 
sidered independently of the effects of the transverse excitation. We 
assume that the longitudinal wave propagates exactly along the direction 
of the external magnetic field so k, E_, and B^ all lie along 

the same direction and f^ has no <J> dependence. There is no induced 
magnetic field so V * E_ = 0 (this is the well known electrostatic 
case) . The appropriate form of the Vlasov equation is then 





>4 


= o 


( 2 . 12 ) 


which combines with Poisson's equation (the appropriate Maxwell equation) 


111 



(2.13) 


to give 


( 7 ) 


EL ( £■> c<j) — 

TOTAL 


B 

— Ore 


Eg ( E, <o) 


where 




- ' k *2 z ' Z ' 


< 4 . 


'(£) -,„ 




(2.14) 


for which we assumed f (v ) is Maxwellian. Symmetry properties of the 
longitudinal dielectric function e^(k,uj) are discussed in Appendix B 
and a brief summary of what is known about the roots of e^Ck.w) = 0 
for real a) and complex k is given in Appendix D. 

The longitudinal part of the external excitation is from (2.11), 
where for definiteness we require iu 2 ^ 0 , 


E (Z,*) - E z (?.-Q £ z cosi^f) 

e*c 


s-t/ 

which has the Fourier-Laplace transform 


isc^ 3 t 

e 


(2.15) 


E ( E E) 


_ - _L 
- e ? 2. + 2 


£ S UU lkL 

(co + s-0/ 2 ) 


(2.16) 


Henceforth we shall drop the summation symbol, although summation over 
s = ± 1 will still be implied. Then from (2.13) and (2.16) we have 


E ( E<E) - e ? 

ror»^ 


_£ 

2 


e lkL 

(oj + e^EjCu) 


(2.17) 
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whose inverse transform is 


E (z,t) = 


TOO 

i fMr f*u 


-iCoT 


2 / -zrr 


J.Jr (^^)e(^w) 

- CO +<• 6 -r 


Noting that for real k and complex w the function e^(k,w) has no 

[7) 

zeros in the upper-half w plane, we lower the w contour (as 


shown in Fig. 2 of Part I) and obtain 


E <*,<■) € 


i g_ 


i &■(*-£) 


(2.19) 


We proceed to calculate the echo response for the ~f~£ case and 
then the £ T case. 

2. 3. i Tj£ Case 

We may rewrite the second-order Vlasov equation (2.4) as 


i y . ~ 9 (&) - 5^ (&) 

Jr * ~ 5 * 7/ ' 7z J 


( 2 . 20 ) 


where 


%( 0) 


=X[( C+ ^ r )^v,(5 i2 


(2.21) 


Equation (2.20) is of the same form as the first-order equation and thus 
has the solution (as shown in Appendix A) 
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X. 

jQ. 


. -27T£ V 


27T 

- + ^(0 + 0')} e/ 


(2.22) 


The Maxwell curl equations combine to give [in analogy to (2.16) of 
Part I] 

E Z(r C Ct ^) = (2.23) 


where 


v/ 
— 2 




5 

^ V 


(2.24) 


The gj term in (2.22) contributes to J_ 2 exactly as the g term 
contributed to in the single transverse excitation case examined 

I 

in Part I. The g 2 term in (2.22) produces what may be considered an 
external current and which we shall designate as where 


J 

— 2 exr 


= sr> 


1 


The part of f 2 that 
depends on g 2 


V 


(2.25) 


Then in analogy to (2.43) of Part I we find 


2 TOT/9L <J~ 


-r 

exr o~ 


cr ) 


We proceed to calculate J- ^ . 
r —2 ext 

For this, the T£ case, we have 


(2.26) 


= 


-/ a Cr.r) • ^ 


TOTSU. 


3V 


r _ » 

from X excitation 

at z = L 


from ^excitation (2.27) 
at z = 0 
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Combining (2.19) with the free -streaming part of (2.73) of Part I, and 
then Fourier-Laplace transforming, we obtain 



-ik? + LUJ 1 C 

e -Z- 




- oo 


27T ^ S UJ Z ) 


-/ o-<fi 


3l 




n viz (-V^) 




EM) e 


k, B + c crus, T 


e r^J~ <ru '/) 


(2.28) 


where 


k = 


- <r (Ct> -MI) 


V, 


(2.29) 


Moving the 9/3v z operator to the far left and interchanging orders of 
integration gives 





aZ-t e 


+ 0 O _ • / 

{(cu-f sco 2 -i-(<ru/) t f </k , E z (k / ') g 


too 


’ Oo 


Zrr e/t',soj^ 


'M e 


i z (~k i-k'+k) 




- oo 


(2.30) 


The integral over z equals 2tt 6 (k ’ -k+kp . 

two integrals , we obtain 

r 






Performing the remaining 




(2.31) 
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where 


c<j^ - - s <^ 2 - a- UJ f 


(2.32) 


Since g 2 (40 % e - * 0 ^ , the part of f 2 associated with g 2 is, after 
performing the <J>' integration in (2.22), 



/ 

i (co- kv^- <ri2) 




(2.33) 


where g 2 (<j>) is still given by (2.31). 

The "external current" that drives the echo is given by (2.25). 
Performing the integrations over 4> and v^ , we obtain 


</■ 

—2 exf 


e 


CO 


c 

c r 


O 




(ui- ArV a ~<rJl) 



(2.34) 


where 


C 










~i(k-k^) L 


l V g (<*-*$ ej^r***) sc*) 


Integrating (2.34) by parts then gives 


J „ 

— 2 £xf 



k <^Vg 

(co ~ cr -f2) 2 


Using (2.35) in (2.26) we find 


(2.35) 
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(2.36) 


cr (*><^) = 

— 2 TOT^i 



The actual response is the Fourier-Laplace transform of (2.36), 


cT (Z a t) 

2 T07/)c 




c^<Uj 


(?r)‘ 


J- (k,Ou) 

— 2 TOT/9 JL 


s£z-iCo>t 

e 


(2.37) 


The primitive contours in the u) and k planes are the same as given 

earlier (in Figs. 2a and 3a of Part I). We perform the k integral 

first. Branch-cut integrals arise from proper consideration of 

^(k-k^jSf^) and e t (k,oj) in the denominator of (2.37) but these give 

a damped response at z = L whereas the echo response occurs typically 

a good distance away from z = L . The dielectric poles in the k 

plane of e^k-k^.su^) and e t (k,w) also produce damped responses 

near z = L . Thus to obtain the echo response we need consider only 

the double pole at k = (u)- 0 ft)/v • Since the integrand of (2.37) con- 

ik fz-L) 

tains the factor e , we pull the primitive k contour above 

for z > L and below for z < L . Then since we have Im u > 0 , we 
obtain a contribution only for z > L , i.e.. 


jr (z,*) 

Z TOT/9C 


CO 

e * c / <yv 

- cr J £ 
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C*>) 
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k£ z {k -k^)Z 


-iik-k^L+ikz^ 


■Z7T (co-u)^) ^{k-k, i<5Cd^)e. r (£ J cJ) 


k = 


Co- 0X1 

V, 


(2.38) 


Lowering the u) contour we obtain a residue contribution at w = » 

J~ (?,*) = 

— 2 TOT~/>£ 


£ c 


^(v g ) ° iUJit (2 * 39) 

'’Z ~ <r<5 ‘j) ^ ^ ^2> ^ < " c k) ^7- C k 3 > c ‘J^) 


where 




f^C^-zO e 


ik(z -0 


[ e r (k, cu) J 


Ar = k- 






Xr 2 = 
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v * 




> (2.40) 


V. 


L = ^^2 £ 

A -o-^+Sl) 


118 



Recall that we are still supposed to sum over s = ± 1 and 
o = ± 1 . When this is done we find [using symmetry properties (B.6) 
and (B.8) from Appendix B] that the s = +1, a = +1 term is the 


complex conjugate of the s = -1, 0 = -1 term so the response is 
twice the real part of either term. Similarly the s = -1, a = + 
term is the complex conjugate of the s = +1, a = -1 term. Thus 
a new dummy index s , we have (arbitrarily using a = -1 ), 


CT * Z 2 ^ 

— 2 TOT4L. s _ + 


e + c <?( i-L) 




( k 3 i^-Z 0 )-icu 3 T 


V 3 6^ ( k^co,) €/ ( k z } s-^ z ) €_ ( ,oj s ) 

ot ^ Jr 


where 


C * — — cv? 




CO — CU, + ? CO 2 


(rE^k-k^e. tk ^~ L) 
€/ C k-k„^i) 
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^■zL 



Co, +JZ 


r Sco z 

k *~ ~^T 
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<cq 3 i-n 


OC ~ SG-Af (Co,i- n) 


p = 5-G/V 


> 


just 

with 


(2.41) 


(2.42) 
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In (2.41) - (2.42) we have explicitly specified which branches of the 
dielectric functions are to be used. 

Equation (2.41) gives the transverse echo response for TJl excita- 
tions. Note that since k^ depends on v z there will be phase mixing 
everywhere except near z = z q , the position of the echo peak. We 
defer calculation of the integral over v in (2.41) and further dis- 

it 

cussion of this, the 7 j case, until after we have investigated the 
case, which we now do. [Numerical results concerning (2.41) for 
some specific cases will be given in Figs. 3-5.] 

2.3.2 J£j~ Case 

The derivation proceeds exactly as in (2.20) - (2.26) but then 
instead of (2.27) we have 


fzW =£.\s (*,*) + 

** ' /yy) 7 ot/9i -1 — 

» L j K 


from ^excitation from ^excitation 
at z = 0 at z = L 


The expression J [*Lj 0TAL (k ,w) + ^ v * B T0TAL (k,u)] • for a trans- 
verse excitation was calculated earlier to be [see (2.31) and (2.39) of 
Part I] 
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Using (2.48) of Part I for E^ xt (k,w) , the inverse transform of (2.44) 
is 


9_ 

/r* 


£ (z,r) 

TOTAL 


h-^-v x b_ < h, 


Tot 



l 

jv 


= I— 
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e e 
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JV Z 
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(2.45) 


The expression for f^k.v.oj) for a longitudinal excitation is 


r . *\ izl 

x. . . 7 ex? > 

y / ' ’ ' ^ e £ {k } o)) lico-kv^ 


(2.46) 


which follows directly from (2.12) and (2.13). Using (2.16) for 
E (k,to) and then taking the inverse transform of (2.46), we find 

~~G Xt 

the free-streaming part of f^ZjV.t) is, for z < 0 , 


/fev) = 


f&(-vj£ 2 (fy ^ e 


/w 2 V-, 


-/ 


( 4 r 5 ^) 


(2.47) 


where k„ = -su) n /v . Equations (2.43), (2.45), (2.47) combine to give 
2 2 z 


+ °o ■£/ 
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This expression is more complicated than the corresponding expression 
for the tJL case, (2.28), mainly because of the factor 



x-Vx ±_ 1 

cr<*/ j 


which in the £/ case was just 3/3v z . 

We proceed to calculate J_ 2 T QTAL^ Z,t ^ in the same manner as used 
for the T/ case in (2.28) - (2.42). This computation is more 
involved than that of the case but the final result is similar to 

the final result of the 7^ case. Omitting the details of this calcu- 
lation, we find: 


r 


= Z 2^- 

— 2 ToT# l S = ±/ 








'c/v z £(^ E/(ty E* (*?) X> A z (^ e 


V E € r ( */, «?) €* (A, e T C 

*~S -/S’ 


(2.49) 


where 


A z ( V = 
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' x f. » g >* 
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£ ( X--/Q e, ( k-k-^) e 


(2.50) 




<vl> = £(vj 
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and the k's, C, oo^, and z q are as defined in (2.42). Equation 
(2.49) gives the transverse echo response for JlT excitations. Note 
that only particles with v^ < 0 contribute to the echo [i.e., note 
/ dv in (2.49)] and that there is a response only for z < 0 

_oo 2 

[i.e., note 0(-v z ) factor in (2.49)]. Note, however, that particles 

with all values of v contribute to the e's in (2.49) and thus 

z 

particles with > 0 may affect the echo shape. 

Having obtained the transverse echo response for both the 7 ^ and 
jtt cases we continue by investigating the consequences of these 
results, (2.41) and (2.49), and then actually performing the integrals 
over v^ to obtain echo shapes for some specific cases. 
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3. 


Transverse Echo Characteristics 


Equations (2.41) and (2.49) represent the transverse echo response 
for the species of particles that has cyclotron frequency Si and zero- 
order distribution function f (v ) . Thus there is an electron echo 
and an ion echo and since z q depends on ft , i.e., from (2.42) 

L 

° + (3.1) 

it follows that in general the electron and ion echo peaks will occur 
at different locations. Also note that since 

o/ 3 = u>, + sea z (3.2) 

it is- possible to have echoes at the sum (s = +1) or difference (s = -1) 
of the excitation frequencies. 

Recall that w 2 is always positive and that u)j may have either 
sign, < 0 corresponding to a left-hand (LH) excitation and > 0 
corresponding to a right-hand (RH) excitation. We adopt the convention 
that RH polarization of a wave means that the J_ (or ) associated 
with the wave rotates in the same direction as electrons gyrate in the 
zero-order magnetic field . (N.B. This differs from the usual 

convention in optics where polarization is referred to the direction of 
propagation of the wave.) 

Now consider just iZ£ echoes, which exist only for z q > L . 

From (3.1) this places restrictions on the allowed values of and 

w_ . Thus to obtain a echo for electrons (Si = -w ) at the differ- 

2 ce 

ence frequency (s = -1) we must have 
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<^ 2 L 


> L 

uy z -(^-oj ce ) 


or 

^2 > (^-^ C e) >0 


(3.3) 


Furthermore, since from (2.41) ^ Re {e + expfittoj-u^t]} , the 

region of u^, co ^ space defined by (3.3) will be divided into two 
regions; for the echo will have LH polarization, for 

the echo will have RH polarization. These two regions are shown in 
Fig. 2 as regions IV and V respectively. In Fig. 2, each region shown 
has distinct echo characteristics; i.e., type ( Tj£ or ), polariza- 
tion, frequency, and direction of propagation (-*■ or •*-) of the waves 
comprising the echo are all indicated. The boundaries of each region 
are obtained in the manner used to define regions IV and V above. As 
shown. Fig. 2 applies to electron echoes but it can be used for ion 
echoes by interchanging RH and LH everywhere, including along the 
axis, where then > 0 would mean a LH excitation. 

The transverse echoes we have been considering offer a much more 
abundant variety of phenomena than those which occur for purely longi- 
tudinal echoes ( We summarize the basic differences and new 

features : 

1. ) Echoes can occur at the sum or difference frequency whereas in 
the longitudinal case the echo occurs only at the difference frequency. 

2. ) The echo location (z ) and frequency (w,) depend on the 
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excitation 


; for ion echoes interchange RH and LH everywhere 




species, whereas in the longitudinal case both the electron and ion 
echoes occur at the same position and at the same frequency. By saying 
depends on species we mean that it is possible to have the electron 
echo at the sum frequency while the ion echo is at the difference fre- 
quency. We emphasize though that given u)j, there is only one 
electron echo (with z q > L or z q < 0 ) and similarly only one ion 
echo. 

3. ) The polarization of the echo may be different than that of 
the initial to 1 excitation. For example, in region IV, the excitation 
has RH polarization while the echo has LH polarization. 

4. ) The propagation direction of the waves comprising the echo 
may be away from or toward the excitation region whereas in the longitu- 
dinal case, the propagation direction is always away from the excitation 
region. The cases where the waves propagate toward the excitation 
region are those which have 0 < u> 3 < w ce and the phenomenon involved 
is the same as that of the negative phase velocity free-streaming waves 
that occur for a single transverse excitation if 0 < w < w , as was 

X Cc 

discussed in Part I. 

5. ) There is a unique echo for every value of whereas in 

the longitudinal case -u^ , oj 2 represented the same echo phenomenon as 

> the only difference being that the echo appeared on the 
opposite side of the excitation region. 

We now comment on a few of the many interesting experiments sug- 
gested by Fig. 2. 

1.) Consider echoes in region IV. With an initial high frequency 
(u)j > w ce ) RH polarized excitation and with to 2 ~ we may generate 
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an echo that is a low frequency LH polarized Alfven wave. 

2. ) We may examine the consequences of echoes comprised of waves 

which have negative phase velocities by considering echoes in regions 
II and V, since in both of these regions we have 0 < < w cg . 

3. ) Note that an echo occurs even if = 0 . This means that 
if the first excitation is simply a fixed transverse electric field, 
i.e., like 


E Ch,*) = £ x £-(z) 

— «xr 


(3.4) 


and if a longitudinal excitation is made at frequency w 2 , then a 
transverse echo with RH polarization will occur at the same frequency 
u) 2 , at the location 

L 

(3.5) 

Thus a longitudinal wave may be "converted" into a transverse wave (of 
the same frequency) which is localzied about a specific position (z q ) 
in the plasma. This novel feature occurs in both regions II and III. 
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4. 


Evaluation of Transverse Echo Shapes 


4. 1 Use of the method of steepest descent. 

We proceed to evaluate the integral over v z that occurs in (2.41) 
for the t£ echo case. [The integral in (2.49) for the echo case 
can be evaluated in the same manner.] Using a Maxwellian distribution. 



(4.1) 


assuming delta function sources [as in (2. 77)- (2 .78) of Part I] so 


(4.2) 


and considering just electron echoes (ft * ) , we find (2.41) becomes 

Cv 


W 0,0 

— 2 total 


Z. 

5 -=+/ 




I 


- / Cj ~t 

e 3 


£0-0 


(4.3) 


where the integral we must evaluate is 



(4.4) 


in which 


= / -f -U 3 (z-L)- 
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€/( tr z ,saj^) 
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A= /f-. 




(4.5) 
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In (4.4) we have explicitly indicated the v dependence of the k 
arguments of the e's. 

We evaluate (4.4) by the method of steepest descent as discussed in 
Appendix E. The appropriate saddle point of the integrand of (4.4) is 

Vo r 


s )(z-* 


2 « 


-] 


. j 7T 

e 


(4.6) 


where 6 = ± 1 for [(Wj-u^g) ( z-z q ) ] < 0 • The primitive contour 
(0 < v /a < °°) in the complex v /a plane is deformed through the 
appropriate saddle point along the path of steepest descent, care being 
taken to indent the contour around any poles caused by the e's in (4.4). 
The three possible pole terms arising from the dielectric functions 
involving o^, u> 2 , (o 3 in (4.4) will be referred to as the pole 1, pole 
2, pole 3 terms, respectively. Thus we find 

3 


I = 1 


PO/fi/T- 


+ 21 I 


< - / 


Pote i 


(4.7) 


The saddle point contribution is 




s/)t>oce 

pa/Asr 


r £ (ou 3 -u, ce )(z-? 0 ) 


A,(v s ) e 




•2 





(4.8) 


[which is valid for | (z-z ) (u)_-w ) /a | ~ 0.12 as discussed in Appendix 

O *5 C6 

E], Assuming single first-order poles in each case, the three pole 
terms are: 
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J = y a 
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A ^ ce 


where v^ , v 2> v^ are defined by 
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and 


a,m = ] / +ii>i£z£2 ' 


4 €s r (4- k n «■£>* = *’ ^ - k* 


Z* s ( *2,^) 




v=v / 
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r ,, ,, ,t,(i 
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>4 
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(4.11) 


*,(v 3 ) = i -3 + i* 3 (z 0 -0 -2 ( 7 ) - 


/f^e^ (*/*») 


* €' (V^) 
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Ar- A-- 
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A prime (*) denotes differentiation with respect to k , and 

UJj-cf ' ce 
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~ ^ce 

V- 
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S<^ 2 

T" 
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(4.12) 


where i = 1, 2, 3. The factor equals ± 1 depending on which 
direction the contour encircles the pole in question, and if there is 
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no pole = 0 . 

In (4.11), note the abundance of terms in A^^) and A^ (v^) as 
compared to Aj(v^) . These terms occur because when the third and 
fourth terms of A^(v z ) , as given in (4.5), are evaluated at pole 2 and 
pole 3 respectively, second-order poles are produced. The new terms 
in Aj(v 2 ) and A^Vj) come from the rather complicated residue 
expressions that result from these second-order poles. 

Concerning the sizes of the various terms in the A^(v^) factors 
in (4.11), we find that all terms that contain an £, £', or £" are 
typically of order one. Noting that 


* 9 (v4| 




(4.13) 


and that the pole terms are important only if Iv^^/al ~ 3 [because of 

2 2 

the exponential factors exp(-v i /a ) in (4.9)], we find | (z q -L) | 
» 1 provided 


L 

Ct. 

t ce 



/ 


(4.14) 


Since typically Lw /c is >> 1, we find (4.14) holds if oj. is not 

cc X 

extremely close to w ( j co, /co - 1 1 ^ » 1) in which case the 

underlined term in each A(v^) in (4.11) is the dominant term. If 
is extremely close to w ( |w,/u) - 1 1 « 1) then the A.(v.)'s 

must be evaluated in full. In any case, the A^v^'s are only amplitude 
factors for the various pole terms. More specifically, the A^(v^)'s 
have no z dependence [except for the slight z dependence in the 
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underlined terms in (4.11)] and, therefore, they can affect the echo 
shape only by affecting the size, but not the shape, of the pole terms. 

4.2 Numerical examples of transverse echo shapes. 

First we explain why we should not evaluate (4.8), (4.9) analyti- 
cally by using asymptotic exmapsions for the Z functions involved. In 
the saddle-point term (4.8), the arguments of the Z functions contained 
in the three e's are all identically equal to v g /a , which is defined 
by (4.6). Typically we shall want to compute the echo shape for 
0 < | (z-z Q ) (a) 3 -(i) ce )/a| ~ 10 say, over which range we have 0 < |v s |/a 
5 1.8 , for which Z(v g /a) is clearly not in its asymptotic limit. 

Similarly for the i th pole term in (4.9), the arguments of the Z func- 
tions contained in the three e's are all identically equal to v^/ a , 
which is defined by (4.10). As mentioned earlier, the i th pole term 
is significant only if |v^|/a ~ 3 [because of the exponential factor 
exp(-v^ 2 /a 2 ) in (4.9)], in which case Z(v^/a) is clearly not in its 
asymptotic limit. Thus, when the saddle point and pole terms are sig- 
nificant, the Z functions involved are clearly not in their asymptotic 
limits. Hence it is rather pointless to rewrite (4.8) and (4.9) using 
asymptotic expansions for the Z functions involved because the result- 
ant expressions would not be valid at the very places we would like to 
use them. (These remarks apply also to the purely longitudinal echo 
case. In particular, the final expression for the spatial longitudinal 
echo case of Gould and O'Neil^ 2 ^ was computed vising the asymptotic 
expansion of the Z function and accordingly their result is essentially 
of little value since it does not apply to any cases wherein the saddle 
point and pole terms are significant.) 
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Thus we have numerically evaluated (4.8), (4.9) for three cases of 

the echoes with the results presented in Figs. 3, 4, 5. (As in 

Part I, we chose “pg/^gg *0.4 and c/a = 1120 , these being typical 

values for the UCLA Q machine.) Choices of for each case are 

given in Fig. 3, which may be compared with Fig. 2 to obtain the echo 

characteristics. (The hatched regions in Fig. 3 will be shown to be the 

regions wherein the pole 1 or pole 3 terms are significant.) Values of 

the poles v^/a were obtained from Appendices C and D. Note that cases 

1 and 2 are echoes at the sum frequency (w^ + u^) whereas case 3 is an 

echo at the difference frequency (oi^ - o^) • 

In Fig. 4 we have drawn the appropriate saddle point contour for 

each case. Recall from Part I that for 0 < qj < w the free-streaming 

1 c© 

pole appeared on the "wrong" side of the contour and therefore never 
contributed to the final response (see Fig. 5 of Part I). For the echo 
calculation, poles 1 and 3 also appear on the "wrong" side of the con- 
tour for 0 < u. < w and 0 < w, < w , respectively; the "right" 
side of the contour being that for which the pole 1 term contributes 

for z < z and the pole 3 term contributes for z > z . However, in 
o o 

the echo case the primitive v z contour is deformed upward, say, for 
z < z q , and downward for z > z q , so a pole will always be enclosed 
somewhere, even if it is on the "wrong" side. In those cases where 
pole 1 or pole 3 is past the end of the contour in Fig. 4 we mean that 
the pole occurs at some v^ where |v^| ~ c and therefore the corre- 
sponding pole term is insignificant because the exponential factor 
2 2 

expf-v^ /a ) [in (4.9)] is zero for all practical purposes. Note that 
any pure imaginary poles (as in case 1) never come near the saddle 
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FIGURE 3. PARAMETER VALUES FOR THREE CASES OF tJL ECHOES. 

c U /to = .4, c/a = 1120 ) 
v pe' ce ’ 
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V 2 /a 
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CASE 3 

1.474 
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2 . 59+i . 142 
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Case I 


Case 2 
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FIGURE 4. SADDLE-POINT CONTOURS FOR THE THREE ECHO CASES OF 
FIGURE 3. 

In each case the primitive contour in the complex v /a 

2 

plane (the real axis for 0 < Re v z < c) has been deformed 
to the appropriate saddle-point contour for z <_ z q (solid 
line) and z > z q (dashed line). The saddle points are 
indicated by X's and the three poles by dots. 
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FIGURE 5. TRANSVERSE ECHO SHAPES FOR THE THREE CASES OBTAINED BY 
EVALUATING (4. 8), (4. 9) NUMERICALLY. (Lw /c = 10) 

C6 
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point contours. 

In Fig. 5 we display numerical results of echo shapes for the 
three cases. From (4.3) we see that 'v- Re I and in Fig. 5 we have 
plotted 
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'S/)£>OL£ 

fiO/A/T 
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I 




(4.15) 


against zc/to using (4.8), (4.9). For these calculations we chose 
cc 

Lc/to = 10 . (For the typical value to /2 tt = 2.25 GHz, c/to ~ 2 cm, 

CG C6 C6 

and L ~ 20 cm.) We note the following results: 

CASH 1. There are no pole contributions so the echo shape is 
determined entirely by the saddle point contribution. The echo width 
is large (compared to a/to ) because the saddle point term's charac- 
teristic length |a/to_-to | is relatively large (because io_ ~ to ). 

O CG J CG 

Also the amplitude is relatively large because the dominant term in 
A x (v s ) , i.e. , 



L c 



= ^ce 



VV 

\^ce ) 


CL 


[see (4.13), (4.14)] is ~ c/a (because / to cg ). 

CASE 2. There are contributions from the saddle-point, pole 1, 

and pole 2 terms. The saddle-point term exhibits a very narrow peak 

because |a/to_-to | ~ a/to (because to_ / to ). Also the amplitude 

of the saddle-point term is relatively small because the dominant term 

in A. (v ) , i.e., (4.16), is relatively small (because to. = to ). 
is i cc 

/ Both the pole 1 and pole 2 terms occur on the z > z q side of the echo 
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but since the pole 2 term is much smaller than the pole 1 term (primari- 

2 2 2 2 

ly because exp[-v 2 /a ] « exp[-v 1 /a ] ), only the pole 1 term makes 
a visible contribution in Fig. 2. The amplitudes of both pole terms are 
also governed by the dominant term of A^(v^) , i.e., (4.16), which for 
this case is relatively small (because ~ ). 

CASE 3. There are contributions from the saddle-point, pole 2, 

and pole 3 terms. The echo frequency ^ ust t> elow the cyclotron 

frequency and the resultant echo has a large amplitude and an extensive 

width. Both the amplitude and width of the saddle-point term are large 

because both |a/a) -w | and |ik (z-L)| [i.e., (4.16)] are relatively 

large. Similarly the pole 3 term has a large amplitude and penetration 

length. The pole 2 term occurs on the z > z q side of the echo, but 

it is not visible in Fig. 5 because it is much smaller than the pole 3 

2 2 2 2 

term (since exp[-v 2 /a ] « exp[-v 3 /a ] ). 

4.3 Physical meaning of the echo terms and where they are 
important . 

The saddle-point term (4.8) represents phase mixing of the free- 
streaming waves that comprise the echo. This term is always significant 
in the echo peak region, which is centered at z = z Q and which has a 
characteristic width or order | a/co 3 -oo ce | . 

The pole 1, 2, 3 terms in (4.9) represent the collective effects 
of the plasma in response to the free-streaming terms that result 
respectively from the u)j excitation, the w 2 excitation, and the 
effective excitation. We emphasize that only the free-streaming 

terms enter the echo calculation; i.e., velocity-dependent free- 
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streaming terms [in f^Cz.v.t) and E^Zjt)] of the form 
exp[i(w/v z )z] contribute to the echo whereas collective oscillations 
of the form expfikj^z] (where is a fixed complex number) do not. 

An important consequence of this is that the collective oscillations of 
f for the first excitation at Wj [as given in (2.73) of Part I] may 
persist past the echo position since they do not interfere with the 
established echo at or produce any new echoes. 

We now determine for what values of Wj, the corresponding 

pole terms are significant in comparison to the saddle-point term. The 
size of a pole term is determined predominantly by the size of the 
factor exp[-v i 2 /a 2 ] which (for |Re » |lm vj) is very small 
fcfr IvJ/a ~ 3 say, and therefore the pole terms are significant 
only for IvJ/a £ 3. For the pole 2 term we consider the least-damped 
root of e^fk.u) = 0 as discussed in Appendix D. For |<*> 2 I < the 

root is pure imaginary and no pole term results. For |w 2 | > w pe the 
root is complex, but only for |u> 2 | ~ 1.08 is |w/"k"a| = |v 2 /a 

$ 3 (see Fig. D-l). Thus only for |w_| 


pe 

1.08 u is the pole 2 term 
pe 1 

significant . Similarly for the pole 1 and pole 3 terms we consider 

the least-damped roots of e T (k,u>) = 0 as discussed in Appendix C. 

From Fig. C-2, we note that |o)-w /"k"a| = |v./a| ~ 3 only if 

U) 2 U > . Thus the pole 1 and pole 3 terms are significant only if Wj 

CC 

or u is very close to w (as shown by the hatched regions in 
3 ce 


Fig. 3). 


3 

Finally, in evaluating (4.4) we consider the 1/v factor which 
causes the echo amplitude to diverge at z = z q . This behavior also 
appears in the saddle point term (4.5) because, as given by (4.6), 
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3 

v -*■ 0 as z -*■ z . The 1/v factor in (4.4) is presumably due to 
s o z 

three factors of l/v z : one for the excitation at [see the free- 

streaming part of (2.73) in Part I], one for the excitation at w 2 
[see (2.47)], and one for the echo response at . Recall that for 
the single transverse excitation in Part I, the free -streaming term 
diverged at z = 0 (as explained in Section 5.2 of Part I). For simi- 
lar reasons the echo diverges at z = z q . In the next section we 
shall show that Fokker-Planck collisions and/or finite-width excitation 
mechanisms remove this divergence at z = z q . 
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5. 


Effects of Fokker-Planck Collisions and Finite-Width Excitation 


Mechanisms . 

We first discuss the effects of collisions and finite widths in 
the excitation mechanisms from a physical point of view. Recall from 
Part I that for a single transverse excitation, the inclusion of a 
finite z dependence in (z,t) reduced the amplitude of the free- 

streaming term and also removed its divergence at z = 0 . Physical 
reasons for the removal of this divergence were given in Section 5.2 of 
Part I. For similar reasons we shall find in the echo case that inclu- 
sion of finite z dependence for the i~ excitation and/or the £ exci- 
tation effectively removes the divergence of the echo response at its 
peak z = z q . 

When velocity-dependent collisions are taken into account, we can 
make the following physical picture. We consider delta- function excita- 
tion mechanisms. Then we note the echo problem as derived above 
reduced essentially to the consideration of free-streaming waves that 
propagate in one direction (the z direction) . Since = B Q e z , the 
particles that comprise these waves have helical orbits that are con- 
fined to lie along the z direction. Each free-streaming wave has a 
definite wavelength and is characterized by the velocity v z of the 
particles that comprise the wave. The effect of long-range collisions 
(such as Coulomb collisions) is principally to "slide" the particles 
along the field lines. Thus collisions can effectively spatially 
reorient the particles in a particular free-streaming wave. And those 
particles that suffer collisions may no longer be in phase with the 
wave. In the limit of many collisions there results a mixture of free- 

\ 
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streaming particles with uncorrelated phases so no net wave is 
discernible. 

The Fokker-Planck collision term, when Coulomb collisions are con- 
sidered, is particularly effective in destroying those free-streaming 
waves whose particles have velocities near v z = 0 . Since these are 
the particles that cause the divergence of the echo at its peak, we 
find that the aforementioned collisions are effective in eliminating 
this divergence. 

We now compute these collisional effects using the Fokker-Planck 
equation 





fv— - n ^ 



£ + 
~ / T*T/)c 


/ 

C 


VX^ T 

- — ✓ TOT/91 


2 7 o" 
:>v 






(5.1) 


where Dj and are the frictional force and diffusion coefficients 

respectively. Considering a single transverse excitation, we recall 
that the free-streaming part of fj goes like 


- e 






(5.2) 


[see (2.73) of Part I], When the Fokker-Planck operator acts on this 
fj we find «. 
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Since generally 



< 




< 

A 

5 V 


V 


•3 V 3 

V 

V 2 


(i a 1,2) we conclude that the underlined term in (5.3) is the 
dominant term in F.P. {fj} provided 


Thus we may write 




(5.5) 


(5.6) 


provided (5.5) holds. 

Now assume that T q TAL a ™i T0TAL damp away in a small 
distance d . Then for z > d and z » | a/o> 1 -co ce | [so (5.5) holds] 
we may rewrite the Fokker-Planck equation (5.1) as 
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(5.7) 


-3? + v 2 Ti - 12 7 ^ = - V/ £ 


Recalling that fj ^ e i< ^ > and Laplace transforming (5.7), we find 




.32 


= / 


^-^ce 




(5.8) 


where we have assumed f^t^) = 0 . Equation (5.8) has the solution 


/(-?) - e 




v e. 




(5.9) 


Thus we have derived an equation which gives f^z) in terms of the 
initial spatial value f 1 (z=0) . This derivation [(5.7) - (5.9)] 
parallels that of O'Neil for the temporal longitudinal case. 

Thus in the echo calculation (specifically for echoes), the 

free-streaming term f^ acquires the factor 

Z 


D z C 6 */ ~ ^ce) 


5 


z / 
Z' eCZ 


over the range 0< z < L . For z > L we have, in place of (5.2), 


/ 




and therefore, in place of (5.6), 


V/ ‘ * 


^/-^ce)Z + r^-^cg) 

r> v ' 




(5.11) 
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Thus for Y~£ echoes, the integrand of (4.4) should include the factor 


-/ 


ef 




L 




V 




(5.12) 


If D 2 does not depend on z we may perform the integrals in (5.12) 
to obtain the factor 


(e<s r - co> <e ) scu z 




L 


(5.13) 


in the neighborhood of the echo peak. In obtaining (5.13) we used 

l. (co, -cc/ ce ) 

( z-D » (**- q * - — ^ L 

The factor (5.13) represents the effects of Fokker-Planck collisions. 
(Numerical examples will be given later.) 

To investigate the effects of finite-width excitation mechanisms 
we assume the excitation fields are given by (2.11) with 


£<Z) -J? 


~ a- 2 

e 


2 Vrr A, 


(5.14) 


(i = 1,2) so 


/ "2 2 


E- W) =& e 


(5.15) 


Thus the integrand of (4.4) should include the factor 
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\ 




(5.16) 


Combining the above results we find that when both Fokker-Planck 
collisions and finite-width excitation mechanisms are taken into account, 
(4.4) should be replaced by 

.2 




I = 




*■ 




(5.17) 


where 


A = 


(cc/~Co ce ) Ay 


3^1 




V 


r = 


C^r^c f sc “z L 


3 


3 C^ 3 -^ ce ) 


(5.18) 


For Coulomb collisions the Fokker-Planck diffusion coefficient as 

(31) 

derived by Thompson may be put in the form 


= 








V, 


w (£) 


(5.19) 


Using the series representation 


(32) 


erf(z) s* 


vw 


?- J- + 


z «/ 


and the asynqptotic series 


(32) 
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+ * • * 


2 » / 


& /~f- ( 2 ") ~ / 





we find 




9. 7T /» f ^ &> -A. 9 

/>n 2 cz. . ? 40r 


M <<: / ( 5 * 20a ) 

CL 


*>*<?$ 


9/rm 

„ sryt 2 CL 




» / 


(5.20b) 


The integrand of (5.17) is now so complicated that in general it is 
only feasible to evaluate the integral by numerically integrating along 
the real axis in the complex v /a plane. However, we can make the 
following useful comments. 

Evaluation of (5.17). 

If (5.17) is evaluated by the saddle-point method, then the pole 
terms that occur are no more difficult to evaluate than earlier. Since 
each pole occurs at a single fixed point in the complex v z /a plane, 
we note that the i th pole term in (4.9) simply acquires an additional 
factor of 


V.z v 5 

e ' ^ 

To obtain a usable result, we assume |v z | » a 
Then the factor (5.21) is 


(5.21) 

so we may use (5.20b). 


e 



F 

(v /c2 ) 8 


(5.22) 
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where 


F = 



((U/ - c^ ce ) L 

/7n ^ d ^ 



(5 . 23) 


We shall evaluate F for a numerical example, but first we consider 

d) f z) . 

the saddle points of e where 

• (^3 _ _v£ 

= <■ v ? v* v/ (5.24) 

We note immediately that even when the v z dependence of D 2 is 

£ 

neglected, the equation that defines the saddle points, g— { 4> (v z ) } = 0, 

z 

is a seventh degree polynomial. It is apparent that it is worthwhile 
to use the method of steepest descent now only in special limiting 
cases. For example, it would be useful to know the height of the echo 
peak. Thus setting z = z q we shall investigate the limits when 
(i) collisions effects are dominant, and (ii) when finite-width effects 
are dominant. 

(i) Collisional effects dominant. 

If we have 


A 2 


v; 


« 


V/ 



s 




where v g is a saddle point to be chosen, then 

0(V) ^ OH _ 

v 

To proceed we assume |v z | » a , then using (5.20b), we 
appropriate saddle point to be 


(5.26) 

find the 
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CL 




(5.27) 


where F is defined by (5.23). Accordingly we obtain 


/. G 5 F 5 




A f (V s > £ 




(5.28) 


where we have evaluated the constant 5/4 ' 2 1.65 . If instead 

j v | << a , then using (5.20b) we find the appropriate saddle point 
to be 

(5.29) 



and 


1 0 = 20 ~ 


2 

-/. 8 2 7 

/)/ ^ 




•5' 




^ce 






7) 


(5.30) 


7 2/7 

where we have evaluated the constant j (10/3 /tt) 2 1.82 . 

Summarizing these results, we find from (5.22), (5.28), and (5.30) 
that Coulomb collisions decrease the size of the pole terms and the 
saddle point term (at the echo peak) according to the factors, roughly 

F 

POLE TERMS ~ (v,/«) 8 (5.31a) 
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ECHO PEAK 


(5.31b) 


(fr) /a £2 

(W /r ) 7<< 7 (5.31c) 

For the parameter values used earlier ( i*> /u> = 0.4 , c/a = 1120 , 

pe C6 

in - 3 

n = 10 cm ), we calculate F from (5.23) to be 




Then using values of Wj, w 2> w^ as given in Fig. 3, we find 


(5.32) 


F = 


00 (A^£f 


CASE 1 


CASE 2 


(/. J»/0 S ) ( 


c 


-) 


3 


CASE 3 


Thus if (Lw /c) = 10 the echo will be essentially destroyed in all 
06 

three cases. However, if (Lw cg /c) = 1 , the case 2 echo would be 
clearly visible. Thus for the parameters under consideration, trans- 
verse echoes could be seen, but only for small values of (Lw /c) . 

CG 

3 

Of course, since F 'v n/T , collisional effects could be reduced 
’ e * 

substantially by going to higher temperatures or lower densities. 
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(ii) Finite-width effects dominant. 


We return to our consideration of (5.24). If 


XD^Vs) 


2 

Vs* 


cz. 2 


(5.33) 


where v g is a new saddle point to be determined, then 


/(v ) 


V 


(5.34) 


and the appropriate saddle point is 


T s 

£ 2 . 




(5.35) 


and 




*/ ( v s) e 


* 


* 






(5.36) 


which may be compared with (7.17) of Part I for the single transverse 
excitation. Note that even if A^ or A^ is zero, the result (5.36) 
still holds. 
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6. Effects of Temperature Anisotropy (T^ f T ) 


We return to the collisionless transverse echo case and investi- 
gate the consequences of having . Assuming an anisotropic 

Maxwellian distribution 






( 6 . 1 ) 


(where a * /2 kT /m , b = /2 kT x /m ) , and carefully reviewing the deri- 
z 

vation for ~tJ? and £t~ echoes, we find that only the following modifi- 
cations are required: 

(i) t/ CASE: replace f (vj in (2.41) by 

vs Z 



(ii) //"CASE: no change. 

(iii) e^(k,ai): no change. 

(iv) ^(kjO)): has already been given for T x ^ in (2.40) of 

Part I. The roots of (k.w) = 0 for real to , 
complex k , and T a ^ T are discussed in 
Appendix C in (C.45) - (C.54). 

Thus for electron echoes in the "fj? case, there is an additional ampli- 
tude factor of 



( 6 . 2 ) 


and for the /t case, there is no change. 

The only significant difference from the previous echo results 
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comes from the possibility of having different dispersion character- 
istics for the transverse waves. In Appendix C, we conclude that for 
w i 0 there is no appreciable change in the positions of the roots of 
e T (k,to) = 0 from their positions in the isotropic case. Thus for 
w t 0 , there is essentially no change in the echo results given 
earlier. On the other hand, if to » 0 and there is a temperature 
anisotropy, then the fire-hose instability may occur. [We obtain the 
usual criterion for the onset of this instability in (C.53).] 

Thus a novel use of echoes is suggested and that is to choose toj 
to be in a stable region and to^ to be in an unstable region so only 
the echo would excite the instability. Use of an echo to examine 
instabilities has the advantage that at the position at which the 
instability is excited, i.e., the echo position, there is no external 
excitation mechanism (grids or current loops, etc.) to be considered. 
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7. Effects of Propagation Just Off-Axis (k^ i 0) 


Up to now, we have considered all waves to be propagating exactly 
along the direction of . In an actual experiment, this may be 
difficult to attain and undoubtedly the excitations will produce waves 
that propagate slightly off-axis so k x would be small but non-zero. 

We now investigate the consequences of having such slightly-off-axis 
wave excitations. 

Consider the longitudinal excitation first. If the excitation is 
made slightly off-axis, then the field of the wave it produces 

will have mainly a z component, but also small x and y components. Thus 


E = (7 * X) 

~ total 

where 6 represents the small transverse electric field ( 1 6 j « |e|). 
As shown in Appendix A, the corresponding driving term for the first 
order Vlasov equation is [for isotropic f Q (v)] 

9 = ^ ^ cos 

5//V (0-1-0') 

+ e | 

* V £ J 

which is to be used in (A. 5), i.e.. 


= — ~ — — 

(e' 27r>r -/) 


{ X SVa/ 




(7.3) 


where 
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Co “ /c j ^ 

nZ 


Thus 






{^S/N<p 


( IZttV \ 

(e -/) 


j ^ ^ _ 

x JVl > 

+ -fe z 


where 


1/ 4$ r - i X -rxv (0 + 0'; 

I ~ J c*V' £ 


(32) 

Using the Bessel function identity 


•f Oo 


-iAs/A'0 ' -cm<f> 

€ = ^ e JfA 

. /?? 

(7) - ~ CO 


together with 


^7 (x) - (-/^ct; (x) 


n integer 


we find 


7=re'^Va)- £ " 7 


-f co 


/7j c- OD 


/ ( V -/>,) 


Then using the series expansion 


(7.4) 


(7.5) 


(7.6) 


(7.7) 


(7.8) 


(7.9) 
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jr(x) 


X' 


2 /»' 


/- 


V ('»+/) 


n integer 
x real 


( 7 . 10 ) 


we find that for |X| « 1 (i.e., | ) much less than the reciprocal 

of the gyroradius , | v ± /Q | ) 

it 

fw = f, (v ? ,v L , <Z )) 
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(7.11) 


2 

which is correct through terms of order X 

To zero order in X we retain only the three underlined terms. 
If in addition 5=0, then only the double-underlined term remains 


//*) 


r £ 


av* 




(7.12) 


in agreement with (2.46) obtained earlier. For 6^0 we include the 
single underlined terms 



= L _/ 

on ■Zi 


-i<£ 


V £ -Jl 


£ 


+ (T 0 


C*j- * £ V +J2. 




(7.13) 


To first and second order in X , we obtain the additional myriad 
of terms in (7.11). Considering the e ±in< ^ factors (where n is an 
integer), we note that as n increases, the terms become higher order 
in X and therefore have smaller magnitudes (since |X| « 1 ). 

Knowing the <j> dependence of fj , we proceed to find the <J> depend- 
ence of e T otaL " We nee< * to ^ now $ dependence of both f ^ and 
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E nTA . to establish what types of echoes and other miscellaneous 

1 U 1 AL 

responses are generated by the new terms created by the off-axis 
excitation. 

If we compute 


= 



(7.14) 


we find that Jj and £ T0TAL no longer obey a simple scalar relation 
but instead they obey 


J~ — cr • E 
— ~~ — r*o7“/9Z. 


(7.15) 


where the elements of the conductivity tensor a could be extracted 
from (7.11) but this information is not needed for our present purposes. 
Noting that Maxwell's equations for the plasma fields give 


£ (4 6 '*-X'V a ) = -97riusiT 

PcasM/9 ' 

[as in (2.38) of Part I] and that £ TOm = *ip LASMA + ^ XT 


(7.16) 
we find 


6 • E 


e xt 


(7.17) 


where the dielectric tensor e is 

9/rictJ 

€= = / + cr 

~ (7.18) 


Thus 



tot /11 ~ ~exr 


(7.19) 


and it follows that in general £jQTAL an< * are not necessar Hy 
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parallel. Nonetheless, if the excitation is made just slightly off- 
axis, both E total and IfiXT wil1 be directed almost along the z 
direction. 

Since [from (2.33) of Part I] 


V e^V = 


a X i 

e — e 

+ Vz 


y, +L<f> y, 

+ p -=■ p +■ e v 

f e - Vi ^ € 




(7.20) 


the only terms of f ^ (<)>) that survive the 4> integration in (7.14) are 

id) 

those with a <j) dependence of e“ or those with no <{> dependence. 
Accordingly Jj has terms of the form 


A 



terms of f ^ ( d>) 


i ±i<j> 

that went as e , 


such as (7.13) 


and 


e 

z 


terms of f^(<j>) 


that had no 4) dependence, 
such as (7.12) 


This spatial behavior of eventually 

of E tqtal [via (7.16) - (7.19)] and we 
(2.27) of Part I] 

terms of 

-T«r„ c >)L L 

+i«sr 

+ q- terms of 


becomes the spatial behavior 
finally obtain [using (2.26) - 




that had no 
that went as 



±i<f> 

e 


(7.21) 


In the right half of Table I, we have summarized the above results 
concerning the 4> dependence of ^(40 and E T0TAL • for a longitu- 
dinal excitation made slightly off-axis. If we repeat the derivation 
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t EXCITATION 


/ EXCITATION 


1 •( 6+EX+ . . . ) 


1 •( E+6X+... ) 


e ^ • ( E+6X+ ... ) 


e ±X **( 6+EX+. . . ) 




e ±l2<,) *( EX+6X 2 +... ) 


±i3d> , c . 2 , 

i ( EX + . . . ) 


e ±l2 *‘( 6X+EX 2 +... ) 


±i3<J> , . 

i * • ( 6X + . . . ) 


E , •— — % i 

-^total 3v 


1 •( 6+EX+... ) 
e~^ • ( E+6X+ . . . ) 


1 •( E+ 6X+ . . . ) 
e ±1<J, -( 6+EX+ . . . ) 


TABLE 1. THE <j> DEPENDENCE OF f.(<J>) AND E. . FOR SLIGHTLY 

T l vy -total 3v 

OFF-AXIS f AND J? EXCITATIONS . 

For the t excitation E = I Q^otal^l md 6 = I ^total 5 z' * 
For the J. excitation 6 = I (§. total )J and E = Uitotai^' 
In all cases X = k,v x /fl . If the excitations are made 
just slightly off-axis, then | X | << 1 and I 5 1 << ll. tota i I and 
the dominant terms are the underlined terms. 
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in (7.1) - (7.21) for a transverse excitation made slightly off-axis, 
we note we would begin by assuming 

£ = € x f x + e y E y + (7.22) 

TOT0L 

where E = (E TO TAL^ 311(1 1^1 <<: I*Lt<3TAiJ ' But s ^ nce ( 7 * 22 ) t * le 

same as (7.1) with 6 and E interchanged, it follows from (7.1) - 
(7.21) that we can obtain the results for a ~f~ excitation simply by 
interchanging 6 and E in the results of the ^ excitation. We 
thus obtain the left-hand side of Table I. The dominant terms for both 
types of excitations are shown underlined. We have not shown the 
resonant denominators that are associated with most terms but these can 
be readily deduced from (7.11). 

We may now easily deduce the new types of echoes that are produced 
by the off-axis terms in ^(40 and ^. TO xal * 9t 7 ‘ We know ^ rom the 

±i (f) 

argument used to obtain (2.10) that if g 2 (4 1 ) % e where 

PzW ~ [ -tota'l ] 

' one ' ^the other' (7.23) 

excitation excitation 

there will be a transverse echo and if g 2 C 4>) has no 4> dependence, 
there will be a longitudinal echo. From Table I, we can easily obtain 
the <j> dependence of g 2 (4>) for all possible combinations of excitations. 
We can then construct Table II which gives the echo type and amplitude 

for the various possible excitations. We note that the dominant echo 

£ 

(which is underlined in Table II) has corrections of order X and 
that there is always an echo of the opposite type of order -g . 

To demonstrate the wide variety of echo locations involved, we 
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££ 


/* 

J?{ E 2 + E6X + ...) 
/•( E6 + E 2 X + . ..) 


TV 


/( E6 + E 2 X + ...) 
-f- ( E 2 + E6X + . . . ) 


If 


/( E6 + E 2 x + ...) 
E 2 + E6X + . . .) 


rr 


' o 

j( E + E6X + ...) 

■4 

f( E6 + E 2 X + ...) 
- 


TABLE 2. ECHO TYPES FOR SLIGHTLY OFF-AXIS EXCITATIONS. 

Shown are the echo types and their amplitudes. 
Here E and 6 refer to either the ^ or the /" 
excitations. The dominant echo types are 
underlined. 
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note that for T£ excitations we may have terms like 


i/m$> |Wwl-/| I A I 

and — EX (7.24) 

v-/7> v- r 

where m, n, r, t are all integers. An echo will occur only if (m+r) 
equals ± 1 (which gives a ~t echo) or zero (which gives an £ echo). 
We may easily deduce that the echo will have frequency 

cu 3 =Co / + Sco z S = ±/ (7.25) 

occur at location 

(asz + tn.) L 

2 

(us 2 + zn.) + s(«' / +*>s£) (7.26) 


with an amplitude factor of 


IH-/| +M (7.27) 

In summary, we note that if the ~t~ and £ excitations are made 
slightly off-axis, the lowest-order echoes for on-axis excitations as 

given in (2.10) will be the dominant echoes, but there will be correc- 

6 

tion terms of order — X also present. In addition, there will be an 

echo of the opposite type ( £ if the dominant echo is /" , or T if 

the dominant echo is £ ) of order jj , as well as many other echoes 

£ 

that are higher order in g- and A , and that involve interactions 
between harmonics of the cyclotron frequency [as in (7.24) - (7.27)]. 
However, typically we may have 



Ar ± cz. 
11 


« / 


(7.28) 
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and then all terras involving X may be neglected. In this case we 
find that only the dominant echo and the echo of the opposite type 
(which is of order ^ ) occur. And for any of the four possible types 
of echo excitation in Table II, we note that the echo of the opposite 
type always occurs at a location different from that of the dominant 
echo. Thus the dominant echo should always be clearly distinguishable. 
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8. Transverse Echoes for B = 0 


For the limit of no external magnetic field, we set w , 
equal to zero wherever they occur in the above calculations. Thus the 
final results for the ~f~i f case and the £ i~ case derived above for an 
external magnetic field, i.e., (2.41) and (2.49), carry over directly 
by simply setting w ce > =0 in them. Thus echoes with frequency 

+- ( 8 . 1 ) 

occur at location 



The appropriate echo type diagram in u)j , space is shown in Fig. 6 

(compare to Fig. 2). Note that many of the interesting features of 

transverse echoes in an external magnetic field have disappeared: 

(1) there is no echo at the sum frequency (1^1 + | o> 2 1 ) ; (2) there 

are no negative phase velocity waves (as occurred earlier for 0 < < 

oj or 0 < to < a) ): and (3) the ion and electron echoes occur 
ce 3 ce 

with the same frequency at the same location. Furthermore, echoes in 

regions 1, 2 represent the same phenomena as echoes in regions 3, 4. 

The only interesting feature that remains is the polarization change 

that occurs for |w | > (wj , i.e., for T£ echoes in regions 2, 3. 

The transverse dielectric function e T (k,co) for = 0 is 

+ 

examined in Appendix C, in Eqs . (C.55)- (C.58) . The basic dispersion 
relation is shown in Fig. C-14a where we note there is an evanescent 
region for |w| < w and a propagation region for |to| > . In 
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LH excitation 


FIGURE 6. TRANSVERSE E 


Fig. C-14b, we have plotted | to/ka | vs. w/w . Considering the 

pe 

roots of E [(w/v ) ,u>] = 0 , we note that w/"k"a = v /a and since 

* ^ Z Z 

| to/ka | ^ c/a (from Fig. C-14b) , any echo pole terms [as in (4.9)] will 

2 2 

be insignificant since they contain the factor exp[-v^ /a ] ~ 

2 2 

exp[-c / a ] which is essentially zero. (Rigorously the distribution 
function should be zero for |v | > c so there should be no pole con- 
tribution whatsoever for | co/ka | c/a .) Thus the pole 1 and pole 3 
terms will never contribute to the echo shape. The only possibility of 


echo shape enhancement is a pole 2 term on the side of the echo 
furthest away from the excitation region (i.e., z > z q for 7 jf 
echoes and z < z q for echoes) . 
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9. Summary 

We have examined the second-order (nonlinear) echo response that 
follows from the Vlasov-Maxwell equations when two spatially separated 
steady-state excitations are made. We found that given the types of 
excitations, the lowest-order echo is given by (2.10). Assuming the 
excitation (2.11), we then derived the transverse echo response up to 
an integral over v for both the t£. case [i.e., (2.41)] and the 
ft case [i.e., (2.49)]. The various characteristics of transverse 
echoes were discussed and summarized in Fig. 2, after which several 
possible experiments were suggested. 

Evaluating the integral over v z by the method of steepest 
descent, we obtained the saddle-point contribution (4.8) and the three 
pole contributions (4.9) for the t£ echo case. Results oftumerical 
evaluation of these terms for three specific cases were given in Figs. 
3-5. The physical meanings of these terms were discussed, and it was 
determined that (i) the saddle-point term is always significant, 

(ii) the pole 1 (pole 3) term is significant only if w 1 = to (w, ~ 
u) ce ) , and (iii) the pole 2 term is significant only if jo^l ~ 1.08 u^. 

We continued by examining several additional effects, some of 

which occur only for transverse echoes in an external magnetic field 

(as we have been considering). First, in examining the effects of 

Fokker-Planck collisions and finite-width excitation mechanisms, we 

found that the resultant integral over v , (5.17) was very compli- 

2 

cated but that it could be evaluated by the method of steepest descent 
in special limiting cases; we obtained approximate solutions in the 
limits when collisional effects dominate and when finite-width effects 
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dominate. We concluded that for the parameter values of Figs. 3-5 

10 -3 

(i.e., w p e / u ce = 0.4, c/a = 1120, n g = 10 cm ) , transverse plasma 

wave echoes could be seen but only for Lu) /c ~ 1 . In general, to 

Cc 

reduce collisional effects, it would be more feasible to observe trans- 
verse echo phenomena at higher temperatures or lower densities. 

Next we investigated the effects of temperature anisotropy 
(T x t T z ) in the zero-order particle distribution function. We con- 
cluded that the only interesting effect was if odj ~ 0 or ~ 0 in 
which case the firehose instability may be initiated. A novel use of 
transverse echoes was then suggested and that was to have in a 

stable region and u> 3 in an unstable region so that the instability 
would be excited only at the position of the echo. 

Then we examined the consequences of having propagation just off- 
axis (k x + 0) . We found that if finite Larmor radius corrections were 
neglected, then the usual first-order echo would occur and this would 
be the dominant echo; but in addition, an echo of the opposite type 
would also occur (but at a different location and with a smaller ampli- 
tude) . It was also shown that if finite Larmor radius effects are 
included, then a myriad of associated higher-order echoes appear. 

Lastly we considered transverse plasma wave echoes when = 0 . 
The resultant echo characteristics were summarized in Fig. 6. It was 
concluded that this case was much less interesting than the case for 

B 4 0 and also that for B = 0 the transverse echo would never be 
-o — o 

enhanced by collective oscillations at frequencies and . 
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APPENDIX A 


SOLUTION OF THE LINEARIZED VLASOV EQUATION FOR B q , k x> k z> ALL NON-ZERO 

We consider the linearized Vlasov equation in k, w space as used 
in (2.17) of Part I, i.e.. 


i (k-v- 


Cu) f - SI — ^ = - 

/ 30 




(A.l) 


where 


9 (0) = X 

syr> 


total 


+ eZ - x 




total 


" . ££ 


(A. 2) 


Equation (A.l) has been solved in a particularly convenient form by 

f 33) 

Fried and Hedrick, J and we give here a brief summary of their method. 
We may write (A.l) as 

i - /r x V x COS0) ^ + a = g (0) 
or 


-i §£ f ^ = £« 


_rz 


or 


-2-/ 

^ 1 


/e^j - 


?<*) 


XL. 


(A. 3) 


where 


Co = Co - 4r g V s 


X 

v 


= XlYo 


XL 


_ 

XL 


J-f = A *V/V0 - V"(0 


r (a. 4 ) 
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Then integrating (A. 3) from <j> to <J> + 2 tt gives 




,H(0^Z7r) _ lHW 


{Z3 + 277- 



JZ> 


e 


or using the change of variable <p’ ■* <j>' + 2u and explicitly 
the H's we obtain 



+<\sy*rf 2 ^ ivv'-tXsv/vCa+t') 

J ) £ 


To express in terms of <J> explicitly we let 


E 

~ TOT/9C 


e x Ex. ^-y Ey 


-A- r~ 

* s i 


and take k^ to be in the xy plane (E^ = e z B Q ) • Then (A. 2) is 


^O) - 


_r 


»£ 




^ + £ y^ av- 




+ ^{^ e ( v «jf -v ? ^) 


+£>* x (v i^ _ 

r X > r 3\/ * 2Vy I 


3V x 


i| “ V f ^ ) 


+ ^ *x ( 


v — - v 

* ?v x * v 


Since 

V, = V A c OS- (Z) 
Vy = r//v <zJ 


evaluating 


(A. 5) 


(A. 6) 


(A. 7) 


(A. 8) 
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we find 


2 V. 




SJN <f) 2 

Ml 


2\ u 


2 Cosrf 3 

$ /N $ ^ V x 2# 


but we may ignore the 3/3<j> operators when they occur in (A. 7) 
3f Q /34> = 0 from (2.4) of Part I. Thus 


where 


?,<» - 


f 


077 


E^cosr ^ 


«{M* &•>} 


> ~l~o /c~i 


+ E y w? +■ <3- 

f ‘'/o 4< r 

‘ ^ CO! "* G 


& = y^-y'i 

^ ^ 2V„ V ? 2 Vi 


which vanishes for an isotropic f Q . 


(A. 9) 


because 


(A. 10) 


(A. 11) 
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APPENDIX B 


SYMMETRY PROPERTIES OF THE TRANSVERSE AND LONGITUDINAL 
DIELECTRIC FUNCTIONS 


From the definition of the plasma dispersion function as given in 
(2.45), (2.51) of Part I we can derive the identities 


sun 

i * 

(r) 


(/■*) 


(B. 1) 
(B. 2) 


\ (where * means complex conjugate) . 

We first discuss the symmetry properties of the transverse dielec- 
tric function [as given in (2.44) of Part I] for w = -crwj as occurs 
in (2.70) of Part I, 


GrA**-*'*',, 9 ') = 


/- 




u 




A 


(-<r<x/) 












>4 A 


Az A 


f (B. 3) 


We may write 




&■ .c& , -J- 

> , , ^ t, 


A x-<2_ 4 ^ 




(B. 4) 


and then it follows trivially that 


(- cr£ f -<rtcs,cr) 

+ 


o~ = + / 


cr= -/ 


(B. 5) 
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I 


Also, from (B.4) and (B.2), we have 



= € r _(~ 

+ 

And since 

d _ _ i 

we find 

<3 (- 4 *) 




-l* 


(B. 7) 


(B. 8) 


where we have used (B.5) on the left side of (B.8). 

The longitudinal dielectric function, as introduced in (2.14) of 
Part II, 






may be written 


~ F 




since 


?'(r) =-2-2 rz ± {f) 


(B.9) 


(B. 10) 


CB. 11) 
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Similarly, the transverse dielectric function as used in (2.73) of 
Part I and thereafter 



Let e + (k,<o) represent e T ( k,u> ) or e & (k,u)) and let W represent 
w/a, a)/A, oj-o) /a, or u*(o ./A . Then in this general notation we 

C6 Cl 

have 




F 



CO 

Ar 



(B. 14) 


Now let to be real and consider the roots of 

€ ± (*,«) »o 

in the complex k plane. From (B.14) we have identically 



Taking complex conjugates gives 



or using (B.l) 



(B. 15) 


(B. 16) 


(B. 17) 


(B. 18) 
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which is just 


e ± (-*<*, w) = o 


Now consider (B.17) again and use (B.2) to obtain 




which is just 


e_(^M =o 


Thus if k is a root of e + (k,w) =0 it follows from (B.19) 

is a root of the same equation, and from (B. 21) that k* is 

e_(k,oj) = 0 . Hence if k + iy (k, y real and positive) is a 
+ 

e + (k,w) = 0 , then the other related roots would be located as 
Fig. B-l. 


(B.19) 
(B. 20) 

(B. 21) 

that -k* 
root of 
root of 
shown in 
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k plane 



FIGURE B-l. SYMMETRY PROPERTIES OF ROOTS OF THE TRANSVERSE AND 
LONGITUDINAL DISPERSION RELATIONS FOR REAL to AND 
COMPLEX k. 

Roots of e + (k,to)=0 in the complex k plane for real to 
[showing that e + (k,to) = 0 implies that e + (-k* ,to)=0 , 
e_(k*,io)=0, and (-k,to)=0] . Both the transverse 
and longitudinal dispersion relations obey this 
symmetry. 
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ROOTS OF TIE TRANSVERSE DISPERSION RELATION FOR REAL u AND COMPLEX k 
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I 

1. Summary 

In this appendix we examine the roots of the transverse dispersion 
relation for real u) and complex k . The major results obtained are: 

1. The cold plasma dispersion relation with negligible warm 
plasma damping (C.9), (C.10) is valid for all w except w = u> ce . 

For 0 ) < w L and to >_ , Im k is essentially zero. For 0 £ to £ 

to it is shown that no simple expression for Im k can be derived, 

C6 

and it is concluded that the root comes from the algebraic terms of the 
Z function involved. The basic dispersion properties are given in Fig. 
C-l. 

2. The locus of roots in the complex C plane of the hot plasma 
dispersion relation is given in Fig. C-3. The loci are classified 
[least-damped or cyclotron-damped (for to ~ to cg ) root ; infinite sequence 
of roots; pure imaginary roots] and the corresponding roots in the k 
plane are then investigated. 

3. Numerical results concerning the cyclotron-damped root for u> 
slightly below <o (Figs. C-4a,b,c, C-S) and for to slightly above 

to e (Fig. C-6) are presented. An expression for the value of this root 
at to * to is given by (C.30). It is noted that for to c to £ to 
both the hot plasma cyclotron-damped root and the "cold plasma" pure 
imaginary root contribute to the response of a single transverse excita- 
tion. 

4. The infinite sequence of roots is examined both numerically 

(Figs. C-7 , C-8) and analytically. It is shown that these roots are 

c 2 

caused by the term 2/ni e ^ that occurs in the asymptotic expansion 
of the Z function involved. Expressions are given for the approximate 
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_ r £ 

values of these roots (C.35) and the value of the term 2/rri e”^ at 
these roots (C.32). 

5. Pure imaginary roots are found to occur for oj l < w < . 

For a), < a) < 0 and w £ to <_ ul. the pure imaginary roots are just 
the cold plasma evanescent roots (C.36), but it is shown that these 
roots occur only in one-half of the k plane. For 0 < u < u) cg the 
pure imaginary root must be calculated numerically. 

6. For reference, the locations of all the roots of e (k,w) = 0 

1 + 

are summarized in Figs. C-9, C-10. 

In addition, several investigations involving the transverse dis- 
persion relation are made: 

7. A physical explanation of the cyclotron- damped root results 
is given using Figs. C-ll, C-12. 

8. The cyclotron-damped root may be computed to very good 

f34) 

approximation by using the two-pole approximation for the Z function 
involved, but only for w £ w cg . For u> >_ w ce the cyclotron-damped 
root passes through the region defined by Im C < 0 and |c| ~ 3 
wherein the two-pole approximation should not be used (see Fig. 13). 

9. Effects of anisotropic temperature are examined (T x i T z ) . 

The usual criterion for the firehose instability is obtained. In gen- 
eral it is concluded that the roots for the anisotropic case are essen- 
tially the same as the roots for the isotropic case, except if u> ~ 0. 

10. The transverse dispersion relation in the limit -*■ 0 is 
investigated. Using Fig. C-14 it is shown that Im k of the least- 
damped root is always essentially zero. 

11. Lastly, in Fig. C-15, it is shown that the roots of 
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e T [(w-od /v ),u>] = 0 (v real) in the complex w plane are centered 

i ^ C6 Z Z 

about to = to ce while the roots of e T (k,kv z +w ce ) = 0 (v 2 real) in 
the complex k plane are centered about k = 0 . 


2. Warm Plasma Results 


2.1 Use of the asymptotic expansion of Z(C e )- 


From the symmetry properties discussed in Appendix B we find it is 
sufficient to obtain only the roots of 


= /- 7 fr-* 

+■ (kd-uj 1 


kat. 




<US . tCO-OJ c 

+ n-o (C.1) 


For Maxwellian distributions, the Z functions in (C.l) have the 
asymptotic expansion 

*2. ( 1*) = - -Z- • + ( 'frr cr (/*) e 

^ f 2 ?* vr s + 


-r 


where 


q;(f) * o 
/ 
2 


for 


r > o 

= O 
< O 


(C.2) 


To good approximation we may use 

2 + cn = - f 

for |c| ~ 3. Using (C.3) in (C.l) we obtain 


(C.3) 
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(C.4) 


co 

k z = Ti 






to(uJ-co') co(uJ+a/ ci ) 


■Yi 


jr? to 
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e W, 


f£ -? 
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(C. 5) 


Then assuming j Re k | » | Im k | so 


^ ct (/&.k') Z + "2^ ( &- k) (j», k") 


and equating real and imaginary parts in (C.4), we find 


ff^k - ± f- 


/ - 


< ^°e 
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C o(to-o j ce ) CU(C<s+OU c ^ J 
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(C.6) 


'tfnto 

Jkr*t k ~ ” 


?c 




-/* 2 
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<£Gj)^ C 


CC^ _ 7* 2 

fe „ /y\ fit 


+ cr 6D /f C 


-r 


(C.7) 


Equation (C.6) gives just the cold plasma dispersion relation while 
(C.7) gives the damping associated with a finite temperature. 


2.2 Problem of o (£ ) . 

+ e 

Evaluation of the o's in (C.7) produces some unexpected results. 
From (C.2) we know a + (C g ) and o + (£^) are each positive or zero. 
Thus for co positive it follows from (C.7) that Im k is positive or 
zero. Accordingly in (C.5) , if io > w ce we have Im > 0 , 
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Im £. > 0 so a (£ ) = a (£.) = 2 . However, for 0 < oj < w we 
have Im £ g < 0, Im ^ > 0 so a + (£ g ) = 0, a + (£p = This * ast 
result is surprising since it claims all the damping is due to the ion 
term in (C.7) . 

Thus we are faced with the following dilemma. If we neglect the 

ion term in (C.l) and use (C.3) for Z (£ ) then for 0 < oj < w we 

cannot legitimately find any root, not even the cold plasma root (C.6). 

This occurs because if we assume Im k > 0 we find a (?) * 0 so 

+ e 

Im k = 0 ; if we assume Im k = 0 , then we find cr + (£ e ) =1 so- 
Im k > 0 . 

The solution to this dilemma must be that, for 0 < w < w , Im k 
is caused by the algebraic terms of Z + (£ g ) . Indeed we shall find 
shortly from numerical considerations that the least-damped root, for 
to slightly below w , occurs in the half of the £ plane wherein 

Cc c 

a + (£ e ) = 0. The point is that for 0 < w < and |c| ~ 3 there is 

no simple expression for Im k like (C.7). 

However, we shall find that Im k for the least-damped root is 
negligible except for | £ g | ~ 3 which occurs only when u) is very near 
to ce . Thus the cold plasma dispersion relation (C.6) with Im k = 0 
is a valid approximation for 0 < w < w ce except for w very near to ce . 
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Then (C.6), (C.7) become 


t&L. A 


= 4 - ^ 
- C 




_/ 

2 


cu(cu - u ' ce ') 


A 


A/rr C*/<*Sa £ 



(C.9) 


(C.10) 


2 

We understand that (C.10) is valid only if (Re k) as given by (C.9) 

is positive. In addition (C.10) is to be used only for toj £ 0 or 

oj > u , in which case we shall soon see that | £ I is so large that 
— ce e 

we may simply set Im k = 0 . According to the above discussion, for 

0 < go < u) ce and |£ e l ~ 3 we may also set Im k = 0 . 

We have plotted the cold plasma dispersion relation (C.9) in Fig. 

C-l for w /to =0.4. Note that k goes to zero at to = 0 and at 
pe ce 

the right- and left-hand cutoffs defined by 


co = f 


to 


ce 




CD* + VCcA 


ce 


7>e 


(C. 11) 


In the propagating regions (<o < co^, 0 < u < w cg , to > to^) k is real, 
whereas in the evanescent regions (u^ < to < 0, <o ce < <o < u^) k is 
pure imaginary. 

In Fig. C-2 we have plotted lo gjo^el vs * w ^ or w pe^ w ce = 

0.4, c/a = 1120) where t = (to-u J/ka and k is given by (C.9). 

Note that |? | is extremely large everywhere (|c | » 3) except for 

to = to , in which case (C.9) is not valid and we must solve (C.8) 
ce 

numerically to obtain the desired root. 


186 



- 1.0 


t — r 


t 1 1 i | — n i i i i i i i | rp i — i i i rn i | 

0.14 0 u/cu ce — |.o |.I4 2.0 

tw L <«U ,W R 

FIGURE C-l. THE COLD PLASMA TRANSVERSE DISPERSION RELATION ( ui /u> = .4 ) . 

pe ce 

2 2 

Equation (C.9) is plotted for (Re k) >0 [solid line] and for (Re k) <0 [dashed line]. 
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FIGURE C-2 . 


LOG 10 l s J 

Plotted is 


VS. W ../10 ( uj /ui = .4 , c/a = 1120 ). 

1 ce pe ce ’ 

log.,, jc ! where x, =(ui-co )/([Re k]a) and Re k is given by (C.9). 
i u q e ce 


3. 


Hot Plasma Results 


3.1 Locus of roots in the t, plane. 

We proceed to investigate the hot plasma results obtained by 
solving (C.8) numerically. We first examine the locus of roots of 
(C.8) in the complex plane. Later we shall consider specific 

types of roots and compute their actual values in the complex k plane. 

In terms of C = (w-u) )/ka , the transverse dispersion relation 
(C.8) may be written 






Ccs Cl/ 2 O* 

re J 


tcj 2 

re - 1 



(C . 12) 


2 

For now we assume the E term on the right-hand side of (C.12) may 

6 

be neglected, i.e., we consider 


C W ■ 


(cu-os c ^<^- 


cc/ccs* a? 
re 


(C. 13) 


(7) 

Then in the spirit of Fried and Gould who examined the roots of the 
longitudinal dispersion relation for real k and complex w ) followed 
by Kuehl, Stewart, and Yeh' - 1 (who examined the roots of the longitudi- 
nal dispersion relation for real u> and complex k ) , we note that we 
can find the locus of roots of (C.13) in the complex t plane by 

v 

plotting solutions of 






= o 


(C. 14) 


Then a true solution to (C.13) exists if for a point on the locus given 
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by (C.14) we have 




i£2.6) 


' (co~co c ^ C 2 ' 
oj os'*- a? 

fie 


(C. 15) 


We have calculated the first several branches of the loci of roots 
given by (C.14) and these are shown in Fig. C-3, where for points on 


the solid 

(dashed) lines 

we have 

Re U e 

3 z + U e )] > 0 (<0). 

Noting 




> O 

for 

<Oj < O } Co > <^ ce 

(C 


Coos' 3 - CZ Z 
L fie J 


< o 

for 

O c oj -c co 

ce 


we conclude that in Fig. C-3 the solid lines are valid loci for u < 0 

oru > u while the dashed lines are valid loci for 0 < u < w 

ce ce 

The various loci in Fig. C-3 may be classified according to the 
types of roots they represent: 

(a) The least-damped root in the plane. (For to near io ce we 

shall refer to this root as the cyclotron-damped root.) In Fig. C-3, 

for to slightly below « the cyclotron-damped root is on the dashed 

cc 

loci just above the real axis; as to approaches to , the root 

C6 

approaches c » 0 . For to = to the root is at C = 0 ; for to 
e ce e 

slightly above to £e the root is on the curved solid loci that begins 
at c; e = 0 and progresses downward. 

(b) Infinite sequence of roots. We refer to the infinite number 
of loops of loci that progress downward in the lower-half plane in 
Fig. C-3. For a given value of u we find there is but one root on 
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each loop. And for l w / w ce l ~ 5 say, each root lies very close to the 
-45° line. This will become clearer when we consider specific 
examples (in Figs. C-7, C-8). 

(c) Pure imaginary roots. Note that the whole imaginary C 
axis is a locus because every point on it is a solution to (C.13). 

More will be said about all of these roots when we consider the 
corresponding roots in the complex k plane. 

2 

Figure C-3 is based on Eq. (C.13) which is valid only when the ? e 
term on the right-hand side of (C . 12) is negligible, i.e., when 



Cu( CU - &S ce ) 

S* 


2 

/>e J 





. « z 



or 


2^.2 


r‘ 

e 


« 


(<i7-4/ ce ) Z c 


^ cz.® 


(C.17) 


This inequality holds for Fig. C-3 (for which |c | 55) except pos- 

6 

sibly when oj ~ w cg . But using = (u>-u> ce )/ka , (C.17) is 



(C . 18) 


which we shall find holds when w ~ u) . Thus Fig. C-3 is valid for 

ce 

all values of w . 

If we wanted to compute the loci for |C e | >> 5 we would have to 
include both terms on the right-hand side of (C.12). In that case we 
could not easily calculate the loci because it is not possible to 
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isolate the dependence from the to dependence in (C.12) to obtain 

an equation in the form of (C.14). Instead, the loci would have to be 
constructed by actually calculating all the roots of (C.12) for many 
values of to . Generally though it is not necessary to know all the 
loci accurately for |? | ~ 5. In fact, the least-damped root for 
u) < oj, , to > to (which has |? | >> 5) has already been computed in 
(C.9), (C.10). And the approximate values of the roots in the infinite 
sequences (for |? e | ~ 3) will be obtained in a following section. 

Thus in general. Fig. C-3 should be sufficient for all applications. 

3 . 2 Roots in the complex k plane. 

(i) Cyclotron- damped root (to slightly below to ). 

c© 

As discussed earlier, the cold plasma result (C.9) with Im k = 0 
is a valid representation of the least-damped root of the transverse 
dispersion relation (C.8) for 0 < oo < to provided |e | ~ 3 . For 

Cv c 

ui ~ to , I C I ~ 3 and (C.8) must be solved numerically. For to ~ to 
ce 1 e 1 ce 

we consider the regions to < oo and to > to separately since our 

° — ce — ce 

main interest lies in the former (w < to ) for which cyclotron-damping 

Cc 

measurements would be made. In the latter (to >_ to ) the damping is 
typically so large that it would be difficult to ever measure it exper- 
imentally. In this section we present values of the cyclotron-damped 

root for several cases where oo is slightly below oo . In the fol- 

ce 

lowing section we briefly consider the cyclotron-damped root for to 
slightly above w ce . 

We have solved (C.8) numerically using the following Nyquist 
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complex root finding scheme. We seek solutions in the C e plane of 



'(o/-^ ce ) 3 c 2 

-4- 

co(Us - cc/ <e ~) 

/l 

a? 


Us* 


« 

L re. j 


L /*€ J 

* 


f = O 


(C . 19) 


[as in (C.12) ] . For each desired value of u> we guess a value for the 

root C (such as the cold plasma value or a point on the appropriate 
6 

locus in Fig. C-3) and then compute 


/V ^ 


/ / e£F (je) 

r / F(r) 


2vt 


(C. 20) 


on a closed contour in the C plane that encloses the chosen 5 • If 

we find N = 0 there is no root enclosed by the contour whereas if 
N = 1 there is precisely one root enclosed in which case we compute 


/ / Ig <*F(f e ) 

2W f F(r e ) 


(C.21) 


which gives a value for the root better than our original guess. Using 

this C and choosing a smaller contour we repeat (C.20), (C.21). 

© 

After a few such iterations we obtain the correct value of the root 
. Then using 

K ~ (C.22) 


we compute the value of the corresponding root in the k plane. 

These calculations were performed using the Fried, Hedrick, 
McCune^^ two-pole approximation for the Z function that occurs in 
(C.19). The results so obtained compared favorably (Re k and Im k 
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agreed within about 2%) to some sample results for which Z + (£ e ) was 
computed exactly. Thus use of the two-pole approximation for this case 
was justified. 

In this and the following sections we consider only roots in the 
right half of the complex k plane since the corresponding roots in 

the left half and also the roots of (k,u)) = 0 are then easily found 

using the symmetry properties summarized in Fig. B-l. Now consider the 
first quadrant of the complex k plane wherein (C.19) has only one 

root for 0 < u < a) , and since the Z function in (C.19) has a (£ ) = 

0 there [using the Z function representation (2.51) of Part I], this 
root comes from the algebraic part of the Z function. We have calcu- 
lated this (cyclotron-damped) root for several values of the parameters 

u) /w and c/a with the results displayed in Figs. C-4a,b,c, and 
pe ce 

C-5 . (The parameters of Fig. C-5 are typical values for the UCLA Q 
machine.) 

We note the following features in Figs. C-4a,b,c, and C-5. First 

consider Re k . If w /w is fixed, then for a finite temperature 

pe ce 

as 0 ) approaches oo ce , Re kc/u) ce approaches a finite cutoff value 

instead of becoming infinite as it does for a cold plasma. Moreover, 

as the temperatures increases (c/a decreases) , this cutoff value 

decreases. If c/a is fixed (as in Fig. C-5) then the cutoff value of 

Re kc/u) increases as oj /u) decreases. These results are 
ce pe ce 

explicitly shown by the cutoff formula 


k < 




ce 


V>r C 


2 


V? 4-/ 


which will be obtained in the next section. 
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FIGURE C-5. THE CYCLOTRON -DAMPED ROOT FOR c/a = 1120. 



Since Re k has a finite cutoff value, the phase velocity w/Re k 

does not approach zero (as it does in the cold plasma case) when w 

approaches co . Instead it remains of the order of the speed of 
cc 

light since 


V 

Phase 


co 


CO 

^ce 






. / c 


co*. 




(C . 23) 


where we have noted from Figs. C-4, C-5 that Re kc/w is typically 
of order 10 for w ~ u> ce . Of course the particles that are responsible 
for cyclotron damping are not those that move at v niI , ot . but rather 

r rlAofc 

those that see a Doppler-shifted frequency w-(Re k)v equal to w 

ce 

Those particles have velocity 


co-<o ce (*o c e /) 

V = T = C 

&ikc \ 

/ 



(C. 24) 


which, of course, approaches zero as oj approaches 


to 


ce 


Now consider Im k . Note in each of Figs. C-4a,b,c (for w /w 

E pe ce 

fixed) that for weak damping, increasing the temperature causes the 
damping to increase whereas for strong damping ( u> extremely close to 
u> ), increasing the temperature causes the damping to actually 
decrease. A similar effect occurs in Fig. C-5 (for c/a fixed) where 
for weak damping, increasing w p e / u ce increases the damping whereas for 
strong damping, increasing u p e / w ce actually causes the damping to 
decrease. These effects will be explained physically in Section 4.1. 
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(ii) Cyclotron- damped root (u) slightly above ui ) . 

In the immediate vicinity of oj = w cg the cyclotron-damped root 
has C e ~ 0 • From the power series representation of the Z function, 

% + (n =-2r+ ir 3 - cc.25) 

we find to good approximation that 


‘Z, (f) ~ - 2 t + i V3F 


(C. 26) 


for |c| « 1 . Using (C.26) in (C.12) we find the transverse disper- 
sion relation becomes 


/rc 


rx-ci 




» 

Wi 





Us. 


* a. 


ce 




(C. 27) 


which is correct provided 



us - Us ce 
fccz- 



« / 


(C . 28) 


Since Ikc/oj | is of order 10, say (for oj z w ), (C.28) holds only 

C6 C“ 

for | [(oj/oj ) - 1] — | « 1 • Note that (C.27) is a true quartic 

C C cl 

equation. 

For oj ~ oj , i.e., for |[ (oj/oj ) - 1] — | « 1 the second and 

C6 C6 3. 

fourth terms in (C.27) may be neglected (as can be justifed a posteriori) 
leaving 



CO c 


i)PP ^ce ^ «- 


ce 


(C. 29) 


201 



In this case the principal root of interest is 


kc 


Us 


Cc> ce 

z 

CcJ „ 

ce 

^ _ 


(vs+i; 


(C. 30) 


which gives the value of the Re k cutoff exactly at u> = w ce as 
mentioned earlier. 

For w near w ce , i.e., for | [( w / w ce ) - 1] ~ 1 , we have 

solved the full quartic equation (C.27) analytically and numerically 
evaluated the solutions (which are quite complicated). We found that 
Re k was essentially constant near o> = w while Im k increased 
uniformly as co increased. Thus the cyclotron-damped root is given 
by (C. 30) for u) = io ce and this root is well behaved as u> departs 
from w ce . 

To study the behavior of the cyclotron- damped root for a) > u) ce 
(i.e. for [(a>/u) ) - 1] - ~ 1 ) we must resort back to solving (C.19) 

ce 3. 

numerically. For uj > co the two-polc approximation for the Z func~ 
tion in (C.19) may not be readily used, as will be explained in Section 
4.2. Thus we have used the exact Z function to obtain the results 
presented in Fig. C-6; shown are the hot plasma results [obtained by 
solving (C.19)], together with the cold plasma result (C.9). We have 
considered frequencies in the range w < w ce to u> > to give a 
complete picture of the behavior of the cyclotron-damped root in the 
evanescent region oo ce < w < . (In Fig* C-6, the cold plasma 

results duplicate a small portion of Fig. C-2 while the hot plasma 
results duplicate part of the w p e^ce = curves in Fig. C-5.) 

The dip (at w ~ oa ) in the hot plasma Re k curve in Fig. C-6 

C6 
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FIGURE C-6. COMPARISON OF HOT PLASMA AND COLD PLASMA ROOTS IN THE 

EVANESCENT REGION u < u < u n . 

ce R 

Comparison of the hot plasma cyclotron-damped root of 

(C.8) to the cold plasma root (C.9). (c/a=1120, 

u» /u> =.4) For u>u , Id k of the hot plasma root 

pe ce ce 1 

keeps increasing linearly (even for w>w n ) . 




flattens out as c/a or <*> /w decreases. 

pc Cv 

Note the large difference between the hot plasma root and the cold 

plasma root for « < w < u> n in Fig. C-6. Actually there is no true 

r ce k 

discrepancy because the "cold" plasma pure imaginary root is also a 
valid hot plasma root (as will be discussed shortly). And, in fact, 
both roots are included in the response for a single transverse excita- 
tion. [In Fig. 3 of Part I, only the cyclotron-damped root of 
e T (k,^) = 0 is shown; actually the pure imaginary ("cold" plasma) 
root of e T (k,^) = 0 also produces a residue term for w ce < u> < w R .] 
Since the dominant root is the one with the smallest value of Im k , 
we see in Fig. C-6 that the "cold" plasma root dominates for most of 
the region w < u> < oj r ; only for w ~ w ce is the cyclotron-damped 
root dominant. Thus the sharp transition from pure real k to pure 
imaginary k that occurs at u * a> ce in the cold plasma case is effec- 
tively smoothed out when finite temperature effects are included. 

(iii) Infinite sequence of roots. 

In the above discussion we examined the cyclotron -damped root of 

(C. 19) . For a) > w this root is on the first curved solid locus in 
v ' — ce 

the lower-half £ e plane in Fig. C-3. However, note that this branch 
is just the first of an infinite number of similar branches that pro- 
gress lower and lower in the L plane. And it is important to 

C 

realize that the lower the branch is in the c plane, the less 

V 

damped the corresponding roots in the k plane are, because k 'u 1/ C e 
[see (C.22)]. Accordingly the most-damped root in the plane is 

the least-damped and possibly most significant root in the k plane. 
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We proceed to examine these roots. 

We present first, in Figs. C-7, C-8, the results of two numerical 
examples which show the infinite sequence of roots in the plane 

and the corresponding infinite sequence of roots in the k plane. 

For these examples we chose values of u) in the interval w < oa < w 

C6 K 

although the resultant distributions of roots in the and k planes 

remain qualitatively the same for all values of w . Note that the 
roots in Figs. C-7a and C-8a are specific points on the loci of roots 
given earlier in Fig. C-3. Also note from these figures that as |c | -*■ 
00 , |k| -*■ 0 . Thus as just explained, the most-damped root in the £ 

plane is the least-damped root in the k plane. 

For all values of w the roots in the k plane approach the 
origin along a 45° line as |C e | approaches infinity. There results 
an infinite number of roots in the k plane whose real and imaginary 
parts approach zero. This is a rather disturbing result since it seems 
these roots might dominate the response of a steady-state transverse 
excitation. Fortunately, however, when the appropriate infinite sum 
of residues is computed (as is done in Part I) we find that the response 
is not significantly affected by the infinite sequence of roots except 
at positions very close to the place of excitation. 

Approximate values of i nfinite-sequence roots. 

We now obtain an expression for the approximate value of the roots 
in the infinite sequence in the c plane. Considering only |t | ~ 3 
we use the asymptotic expansion of the Z function [as given by (C.3)] 
in the transverse dispersion relation (C.12). We restrict our attention 
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FIGURE C-7. INFINITE SEQUENCE OF ROOTS FOR w/w = I. 001 . 

ce 

Shown in (a) are the roots in the plane and in (b) 
the corresponding roots in the kc/w plane (c/a=1120, 

C6 

. Not shown in (b) is the root 14.86+i9.56 

pe ce 

that corresponds to the root 0. 133-i0.0856 in (a). 
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(C . 31) 


to Im £ <0 so o (^ ) = 2 and therefore 

e + e 


2 (4) * - T + e r * 

+ 4 

2 

We shall (1) show that the term 2/Hi exp[-r ] causes the infinite 
sequence of roots, (2) obtain an expression for the approximate value 
of these roots, and (3) obtain an expression for the term 2/fri exp[-C g ] 
evaluated at these roots (the latter expression for use in Part I, in 
calculating the infinite sum of residue terms associated with the 
infinite sequence of roots) . 

Using (C.31) in (C.12), the transverse dispersion relation becomes 


21/^i e 




e _ 


£ 




o = 


£ 5 


Oj(as-^ ce ) 
(^-^ce)*c 2 1 


(C . 32) 


We restrict our attention to the fourth quadrant of the plane 

(since any root in the fourth quadrant has the corresponding root 

-C e * in the third quadrant, as follows from Appendix B) . Then 
£ = | Re C | - i |ln> C I and the left-hand side of (C.32) may be 

written as 


-/ 

2V?( € e = 2Y 7f<2 




(C.33) 


O 7 

Note that the factor exp (-[(Re c ) -(Im £ ) ]) equals unity when 
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is on the -45° line in the £ g plane, and that this factor changes 
size very rapidly (for |c | ~ 3 ) when £ departs from the -45° line. 
Thus even though the magnitude of the right-hand side of (C.32) may vary 
considerably, by only slightly departing from the -45° line, the left- 
hand side of (C.32) can match it and a root occurs. (This explains why 
the infinite-sequence roots are always near the -45° line.) Also, at 
a root the phase factor in (C.33), [2(Re £ g ) | (Im £ g ) | + tt/ 2] , must 

equal the corresponding phase factor of the right-hand side of (C.32) 
plus any integral multiple of 2 tt . (This last consideration explains 
the multiplicity of roots.) Thus we have shown qualitatively how the 
term 2/iri exp[-£ g ] produces the infinite sequence of roots. 

To obtain an expression for the approximate value of the roots of 
(C.32), we note that a root occurs whenever the magnitudes (and phases) 
of the right- and left-hand sides of (C.32) are equal. From the above 
discussion we know the roots are near the -45° line. Thus we set 
£ g = pe*^ 4 ^ (p real, > 0) everywhere in (C.32), (C.33), except in 
the factor exp (-[(Re £ g ) - (Ini C g ) ]) . Then, at a root, the phases 

of the two sides of (C.32) must be equal, i.e., we must have 


P 



PHrisa of < 


De 


■ 7r 'in 


+- 


P 


P 3 


+ 2 


/K7T 


(C. 34) 


where n = 0,±1,±2, For |n| > 3 , we have 

tively we may say that the roots in the infinite 

• 7T 


4 = P e 


~ l v - 


2 

p z 2m . Alterna- 
sequence are given by 


P - V ’ 2 «,7T 


= V, S', 


(C . 35) 
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2 

and the value of the term 2v*iri exp[-S ] at such roots is simply 

6 

given by the right-hand side of (C.32) evaluated at the roots. In 
(C. 35) n starts at n = 3 because in using (C.31) we assumed 
|S | ~ 3. It follows from Fig. C-3 that only the first two roots of 
the infinite sequence (which have |? e | ~ 3) need to be calculated 
accurately while the rest of them are given approximately by (C.35). 

[In the above, note that it was not necessary to express Z (E ) 
in terms of Fresnel functions for along the s -45° line and then use 

asymptotic expansions for the Fresnel functions to obtain the approxi- 
mate values of the intercepts of the loci of (C.12) with the -45° line 
(as was done for the longitudinal case in Refs. 7 and 8) since we were 
able to obtain approximate values for the roots and even the value of 

the term 2*^Ti exp[-C ] at these roots by simply considering the 

€ 

asymptotic expansion of the Z function.] 

The infinite sequence of roots in the S e plane, as discussed 
above, occurs only in the lower-half C plane (Im C < 0) . [For 

w C 

2 

Im £ > 0, a (5 ) = 0 and there is no 2/rFi exp(-E ) term to pro- 

e + e c 

duce the infinite sequence of roots.] We now determine where the cor- 
responding roots in the k plane are. Since k = ( w_u ce )/£ e a , 

Im t <0 implies Im k < 0 for u < w and Im k > 0 for w > w . 
e ce ce 

Thus the infinite sequence of roots occurs in the lower-half k plane 
for u) < w ce and in the upper-half k plane for w > o> ce . 

(iv) Pure imaginary roots. 

The hot plasma dispersion relation (C.8) has the pure imaginary 
root(s) 
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(C. 36) 


Ar = :fci 




_c 2 \ CeJ (uy - Co c e ~) 



(for co < co < 0 and co < co < co n ) but only if the root is in the 
v L ce K 

half of the k plane wherein Im ? > 0 . In the half of the k plane 

wherein Im C <0 there may be a pure imaginary root of (C.8) but it 
e 

2 

is not given by (C.35) because the term 2/rri exp[-^ e ] (which goes 

as exp[+|c | ] for C pure imaginary) was not adequately considered 
6 6 

in obtaining (C.36). 

2 2 

For a pure imaginary root, C = -|C | and the hot plasma dis- 
persion relation (C.12) may be written as 




l 


CoujJ' « 2 
„ Pe 

“T 

a/* Jl^l 

pe J 


(C. 37) 


Note that the right-hand side of (C.37) is positive for to < 0, to > w. e 
and negative for 0 < co < co ce . Thus we distinguish the following two 
cases . 

(1) co < 0, co > (0 : From Fig. C-3 we know that for co < 0, 

CG 

CO > co ce any allowed pure imaginary root has C e = +i U e l in which 

case a (c ) =0 . And we note in general that the right-hand side of 
+ e 

(C.37) is very large (except for co extremely close to co £e ) . Thus 
any root of (C.37) will have |? e | ~ 3 in which case the cold plasma 
results are valid; i.e., for to < co^ > to > to^ there is no pure 
imaginary root whereas for co^ < co < 0 and to ce < to < to^ there is one 
pure imaginary root, given by (C.36), that occurs in the half of the 
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k plane wherein Im ? > 0 . 

r e 

(2) 0 < to < 0 ) : From Fig. C-3 we know that for 0 < to < to 

Cw ( 

any allowed pure imaginary root has C = -i|? | in which case 

0 6 

a + (C e ) = 2 . In this case (C.37) may be written as 






CLSCCJ < 2 .^ 

/>e 






Cu (ca-CC^g) 
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Cc/~ 2 

/°e 


✓ 


(C. 38) 


[in which Z represents the principal part integral of Z + (£ g ) as 
given in (2.51) of Part I]. Since exp [-*■ | I ] is a rapidly growing 
function of k | we note that the root of (C.38) will typically have 
|C 0 | - 4 , say. In general, (C.38) must be solved numerically. 

Thus cold plasma theory predicts pure imaginary roots only in the 
two evanescent regions to L < to < 0 and u ce < to < to R . Hot plasma 
theory also predicts these roots occur (but only on the side of the k 
plane wherein Im > 0 ). In addition, hot plasma theory predicts a 
pure imaginary root for 0 < to < to ce . 




3 . 3 Summary of roots in the complex k plane. 

The cold plasma dispersion relation (C.9) with negligible hot 
plasma damping (C.10) is valid in the two propagating regions to < to, 

Li 

and to > w R . In the propagating region 0 < to < to ce except for 
w ~ w ce a root a l so given by (C.9) but with no damping term. For 
w ~ w ce we must solve the full hot plasma dispersion relation (C.8) 
numerically. Some numerical calculations of this root are given in 
Figs. C-4, C-5 for to slightly below u> ce and in Fig. C-6 for to 
slightly above to 

ce 
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In the half of the k plane wherein Im £ < 0 there are an 

infinite number of roots of (C.8) which asymptotically approach the 
origin in the k plane along a 45° line. Some numerical results con- 
cerning these roots are presented in Figs. C-7, C-8, and the approxi- 
mate values of these roots are given by (C.35). 

Also, pure imaginary roots occur for some values of oj . For 
to < a), and to > to n there is no pure imaginary root, whereas for 

to < to < 0 and to < to < to there is a pure imaginary root given by 

(C.36) (that occurs in the half of the k plane wherein Im £ > 0 ) . 

For 0 < to < to ce there is also a pure imaginary root, which can be 

found by solving (C.38) numerically. 

We have summarized these results in Fig. C-9 which shows the roots 
of (k, ) * 0 for all choices of to . We have used the symmetry 
property (from Appendix B) that all roots in the left-half k plane 
are mirror images about the imaginary axis of roots in the right-half 
k plane. 

In Fig. C-10 we have summarized the roots of (k,u>) = 0 and of 

+ 

(k,w) = 0 (shown by crosses) in the upper-half k plane. The 

roots of (k,oj) = 0 are mirror images across the real axis of the 

roots of e.j, (k,oo) = 0 (this symmetry property was also given in 
+ 

Appendix B) . 
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(x) |^<Cl)<0 


Cl)<OJ^ 


0<U)<(jJ 

ce 



FIGURE C-10. ROOTS OF e T+ (k,a>) = 0 [DOTS] AND 
ROOTS OF e T _(k,a>) = 0 [CROSSES] 
IN THE UPPER-HALF k PLANE. 
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Special Considerations 


4.1 Physical explanations of cyclotron-damped root results. 


In this section we want to explain physically the (cyclotron- 
damped) Im k curves of Figs. C-5, C-6. As discussed earlier, cyclo- 
tron damping is caused by those particles whose Doppler-shifted frequency 
is approximately equal to the cyclotron frequency. Roughly, the damping 
particles are those which have v z within 




where v is defined exactly and Av approximately by 

- (/&k) V = ^ ce 


\oij - (^/OCv-kav^ - <*' 


(C. 39) 


in which a is some small number (like 0.001) 


to -co e 

V » — 

AV = ^^ 5-g- 
k^k 


(C .40) 


and the number of damping particles N is roughly 


kV ~ s* -f ( W A V 


(C .41) 


Using a Maxwellian distribution for f (v ) 


V ,)2 


~ 'free e 


(C.42) 


we have 
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(C . 43) 


In each of Figs. 5a,b,c we consider n and B q fixed and examine 

the effects of varying T . Using (C.39) - (C.43) we show in Figs. 

6 

C-ll the particles responsible for the damping. Note that increasing 

T broadens the tails of f (v) and lowers the center peak. For weak 
0 o z 

damping, v is out on the tail of f Q (v z ) and therefore increasing 

T g increases N . For strong damping (a> ~ w cg ) v = 0 and therefore 

increasing T g actually decreases N . This explains the behavior of 

the Im k curves in Figs. 5a,b,c. 

In Fig. C-6 we consider n and T e fixed and examine the effects 

of varying B . In this case increasing u> /u> means decreasing 
o pc cc 

B q (but to maintain the same horizontal position in Fig. C-5 we must 
keep w/(D fixed). In Fig. C-12 we show the particles responsible 
for the damping; note that v and Av change while f Q (v ) remains 
constant. [This may be seen by writing (C.40) as 


Co ~ 


1 c: 

4r 

i ^:) 
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oC ^ c e _ _ 

oC C 


<£a. .fc 1 

f <^0. £cz) f 
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\ ^ce J( 

-c e) 


(C. 44) 


and noting in Fig. C-6 that (Re k c/w pe ) (“pe^ce^ ^ “pe^ce For 
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FIGURE C-ll. PARTICLES RESPONSIBLE FOR ELECTRON CYCLOTRON DAMPING 

WHEN n AND B ARE FIXED AND T IS VARIED, 
e — o e 

(a) Weak damping. 

(b) Strong damping. 
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weak damping, the number of damping particles depends mainly on v 
which is out on the tail of f Q (v ) . Thus decreasing B q decreases 
v thereby increasing N . For strong damping (w ~ w ) v ~ 0 and 
the number of damping particles depends mainly on the velocity spread 
width Av . In this case decreasing B q decreases Av thereby 
decreasing N . This explains the behavior of the Im k curves in 
Fig. C-6. 


4.2 Why the two-pole approximation for Z(£) should not be used 
in calculating the cyclotron- damped root for <d > w . 

(34) 

In this section we explain why the two-pole approximation^ ' for 

Z(C) should not be used in computing the cyclotron-damped root for 

oj > oj .As shown in Fig. C-13, the locus of this root begins at 
— ce 

£ = 0 (for cu = a) ) and progresses downward as o) increases. As will 

e ce 

be explained, the two-pole approximation is in error in the region 

defined by Im £ < 0 and |c| ~ 3. Since the cyclotron-damped root 

passes through this region for u) > u) , the two-pole approximation 

cc 

should not be used in computing the root in this region. 

For Im ^ > 0 the two-pole approximation matches the first two 

terms in the power series [Z + (£) ~ it/rF - 2^] for small £ and the 

first term in the asymptotic series [Z + (s) ~ - I/?] for large £ . 

-C 2 

For Im £ < 0 one is instructed to add on the term 2/ifi e 

Clearly this is correct only for |c| ~ 3 because only the asymptotic 

_ £ ^ 

series (which is valid for |t;| ~ 3 ) has the term i/ir 0 + (t) e 
The power series has no additional term for Im £ < 0 . Thus, as 
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FIGURE C-13. REGION IN TOE 5 PLANE WHEREIN TOE TWO-POLE APPROXIMATION 
FOR Z(C) SHOULD NOT BE USED. 

Shown are regions of the c plane wherein the power 
series and asymptotic series of Z + (c) are useful 
approximations. The two-pole approximation 34 for Z (O 
should not be used for Im ; < 0 and Is |c| s 3 (see 
text). Also shown is the locus of cyclotron-damped 
roots for u 2 u , as given in Fig. C-3. 



presented, the two-pole approximation is in error in the region defined 
by Im 5 < 0 and |c| ~ 3 . 

We are left with the problem of deciding specifically when the 

j-2 

term 2/rri e"^ should be added on. (Presumably this should occur for 
Im ? < 0 and |?| somewhere between 2 and 3.) Furthermore, at the 
place where the term is added on, the two-pole approximation for Z(£) 
will effectively suffer a discontinuity, whereas the true Z function is 
continuous. To avoid this artificial discontinuity, the exact Z 
function should be used in the troublesome region. In summary, for 
Im £ < 0 and 


for 

sTi 

l A 
>-* 


add 

nothing; 

for 

i * Id s 

3 

use 

the exact Z function; and 

for 

3 ~ Id ~ 

3 

add 

r 2 

2/ni e~^ . 


In actual practise, such as in computing the cyclotron-damped root for 
w u) ce it is probably least confusing to use the exact Z function for 
all Im C < 0 . 


4.3 Effects of anisotropic temperature (T x i T ) . 


All of the above work was done for T. = T . For T, i T we 

X z X z 

have from (2.40) of Part I that, neglecting ions. 
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Taking f Q (v z ) to be Maxwellian, we obtain 
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(C .46) 
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where a = 1 - (J x /T^) and Z'(C) = Z(C) . Assuming | ? g | ~ 3 
and | Re k j » j Im k | , and proceeding as in (C.3) - (C.7), we find the 
least-damped root of e^, (k,w) = 0 to have 
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(C.47) 


■aiorj ^ — 


^-r ,C0 t -*V 

2c + l J 


/r- 


or 


ft) 


* 

/>e 


J 


a 


(C . 48) 


If we assume 
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(C . 50) 


then as we found earlier for (C.10), Eq. (C.48) is valid only for 

a) < a) or a) > oj d . For 0 < to < u> (w t u ) we must ignore (C.48) 

L K C6 CC 

and conclude that the root is undamped and given by (C.47). 

In general, the only situation of real concern is if a is such 

that 
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(C. 51) 


a 



which is 
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(C.52) 


For u> = 0 this is 


('-£) 


5 


(C.53) 


which is the usual criterion for obtaining the firehose instability. 
In (C.46) note that the coefficient of Z'(£ e ) 


(18) 


fc of CZ 
leu 



(C. 54) 


is typically very small, in which case the roots of (C.46) are essen- 
tially the same as the roots of the isotropic case (T^ = T^) . Only if 
oj = 0 does (C.54) become significant, in which case the roots of (C.46) 
would differ significantly from those of the isotropic case. 


4.4 The transverse dispersion relation for = 0 . 


In the limit of no external magnetic field, the transverse disper- 
sion relation (C.l) becomes 
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(C . 55) 
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Neglecting the ions and proceeding as in (C.3) - (C.7) we find the 
least-damped root to have 


= ±# /- 
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J^yi A 


'frr Cucc/^£ 
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(C. 56) 


(C.57) 


Equation (C.56) is plotted in Fig. C-14a. Note that both w > 0 and 
a) < 0 have the same dispersion characteristics. From (C.56) we have 


c 



(C . 58) 


and this is displayed in Fig. C-14b. Since c/a is typically very 
large (» 3) we note that |u)/ka| is also very large (except if 
oj ~ 0) ; this justifies our use of the asymptotic expansion (C.3) to 
obtain (C.56) - (C.57) and it also means that the damping as given by 
(C.57) is zero for all practical purposes. 


4.5 Roots of e_. ( Coo-co /v ),«] =0 in the complex u) plane and 

1 ^ C6 2 

roots of Et (k,kv +<jj ) = 0 in the complex k plane (for v 
i ^ z ce “ 

real) . 

Substituting k = (w-u) )/v in (C.8) we find 



CO - cu zt 
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(C . 59) 


is 


225 




FIGURE C-14. THE LEAST-DAMPED ROOT OF THE TRANSVERSE DISPERSION 

RELATION FOR B - 0. 

— o 

(a) Plot of (C.56); the solid line applies for (Re k) 2 >0, 

2 

the dashed line for (Re k) <0. 

(b) Plot of (C.58) « 
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(co-to^ -U> z (cc/-^ C ^) 
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c 2 ^ 


*+(£) ' 


(C.60) 


Note that the argument of the Z function in (C.60) is real (v is real) 
and has no w dependence. Thus, given a value for v z , (C.60) is a 
cubic equation in w which has three exact analytic solutions. In the 
cold plasma limit (a -*■ 0) , and for l v 2 / c l <<: 1 all three roots of 
(C.60) are w = co ce . For a hot plasma, we have calculated the locus 
of roots of (C.60) for a typical case in Fig. C-15b; note that all 
three roots remain very close to uj = w . 

Cv 

Substituting w = kv + u) in (C. 8) we find 

z cc 



(C .61) 


is 


+ # s + (£)] + [4r* * + ©] = ° 


(C.62) 


which is a cubic equation in k . In the cold plasma limit all three 
roots of (C.62) are k = 0 . For the hot plasma parameter values used 
in Fig. C-15b we have plotted the locus of roots of (C.62) in Fig. 
C-15a. Note that the roots of (C.62) may be obtained from the roots of 


(C.60), or vice versa, by use of the change of variable 

/ 

A * w 


(C.63) 


In particular Fig. C-15a may be obtained from Fig. C-15b, or vice versa, 
by use of (C.63) . 
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APPENDIX D 


ROOTS OF THE LONGITUDINAL DISPERSION RELATION FOR REAL 0) AND COMPLEX k 


From the symmetry properties discussed in Appendix B we find it is 
sufficient to obtain only the roots of 


^ ^ < { ta ) ■ ° 

Using the asymptotic expansion of the Z function 




y - Tj* + £ 


and proceeding exactly as in (C.l) - (C.7) we obtain 


k = ± 
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Z VF u> 5 or 
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Cc/ Z 
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<z 3 
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/ ^ -/ 2 


(D.l) 


(D.2) 


(D. 3) 


(D.4) 


where o = a (c ) = a (£.) and t = w/ka , £. = w/kA . None of 

* * v ' 1 C 1 

the problems associated with a (£ ) of the transverse case occur in 
this, the longitudinal case. 

Considering Im £ < 0 so a + (£ e ) = 2 and neglecting the ion 

terms in (D.3) - (D.4) we find 


£ / us .1 e (D.5) 
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(D.6) 


where k , = f2 (w /a) is the electron Debye wavenumber. Equations 

Cl pc 

(D.5) - (D.6) are valid only for |c | ~ 3 since the asymptotic expan- 
sion of the Z function (D.2) was used in their derivation. From (D.5) 
we have 




CO 
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(D.7) 


Setting C = 3 in (D.7) and solving for w we find u> = 1.08 w 

e pe 

Thus only for the narrow range <*) < u ~ 1.08 u are the asymptotic 

pG p6 

formulas (D.5) - (D.6) valid. 

For u) ~ 1.08 a) we must solve (D.l) numerically. Neglecting 
pe 

the ion term, (D.l) is 


£V(£) ' 


pz 


(D.8) 


(8) 

Kuehl, Stewart, and Yeh^ 3 have plotted Re k/k^ and Im k/k^ vs. 

w/u)^ corresponding to the first four least-damped roots of (D.8) in 

(9) 

the C plane. Also, Derfler and Simonen have numerically tabulated 

v 

the first four roots of (D.8). In both cases, for w - 1.4 w the 

pe 

damping of the higher-order roots becomes comparable to, and even less 

( 2 ) 

than, that of the "least-damped" root. However, Gould ^ 3 has evaluated 
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the first few residue terms that occur for a longitudinal excitation 
and concluded that use of just the "least -damped" root is probably a 
valid approximation. 

For computational purposes, it is especially useful to know the 

value of the least-damped root for w near to . Since for to < 

r pe pe — 

to ~ 1.25 a>p e the graphs of Kuehl, Stewart, and Yeh may not be read 
accurately (especially Im k/k^ ) we have also calculated the least- 
damped root of (D.8) for w <_ to _< 1.5 w and this is displayed in 
the top of Fig. D-l. Also shown are the asymptotic results, (D.5) - 
(D.6). In the bottom of Fig. D-l we have plotted |c | vs. to/to 

v DC 

for the true least-damped root of (D.8). Note that as deduced above, 

the asymptotic results (D.5) - (D.6) as shown in Fig. 10 are good 

approximations only for w £ to ~ 1.08 to^ g . 

To obtain the results in Fig. D-l, the true least-damped root of 

(D.8) was calculated using the correct function Z'(£) rather than the 

f341 

two-pole approximation. The latter was tried but found to be very 

poor in this case. More specifically, when the two-pole approximation 

was used in (D.8) for 1.15 oj < u> < 1.5 <*> we found the least-damped 

pe — — pe r 

root had Re k/k^ in error from its true value by 3%-18% and Im k/k^ 
in error from its true value by 15%-40%; and when u went below 
1.15 a) , Im k even went negative. The reason for this poor behavior 
is presumably due to the fact that the two-pole approximation for 
Z'(S) matches only the first term of the asymptotic series for Z'(^) 
accurately whereas the basic asymptotic form of the longitudinal disper- 
sion relation (D.5) - (D.6) requires for its derivation the first two 
terms in the asymptotic series of Z’(£) . 
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Asymptotic 
(D5) / 


Rek 

k d 


True root of 
(D8) 


Imk 

kj 


Asymptotic 
(D6) _ 


True root of 
(D8) 



FIGURE D-l. THE LEAST-DAMPED ROOT OF THE LONGITUDINAL DISPERSION 
RELATION FOR REAL u AND COMPLEX k. 

Top: Comparison of true and asynptotic results. 

Bottom: Plot of |u/ka| for the true root. 




APPENDIX E 


METHOD OF STEEPEST DESCENT AS USED IN EVALUATING 
BRANCH-CUT AND PHASE-MIXING INTEGRALS 


We present here a brief summary of the method of steepest 
(351 

descent as used in both Sections I and II. 

Consider an integral of the form 


r 4- 0 (W) 

I - J yt(w) e c(V4 


(E.l) 


where F(W) is a slowly- varying function of W and <j>(W) is a 
rapidly-varying function of W . Then the integrand of (E.l) has 
saddle points at W q where 

P'W =° (E. 2) 


We assume the initial contour C is along the real axis in the complex 
W plane. Then rewriting (E.l) as 


I € _5Z ^ w) = f /=(w) e 

c 

we note that if we deform the contour 
W = W q along a path such that 


9>(w) - 0 O 0 ) 

W 


(E.3) 


C through the saddle point at 



0 fw) - 0 ( w a ) L 


= o 


(E.4) 


i.e., a path of constant phase, then the integral in (E.3) may be 
approximated as follows. 

We first expand <j>(W) about W = W q 
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(E.5) 


0<V) = 0*(n£)+ -£r {'w- w of + &" / (K') (v-w e ) 3 + •-• 


so 

0 (w) - 0 ( O a* X <t>"( w e ) ( Uf-W 0 ) z 


(E.6) 


We then specify the desired path of constant phase C' by imposing the 
transformation of variable 


(&-K,')* = - i~ 2 


where T is a real variable. Thus (E.3) becomes 


- <p (Vo) 

.ze 



(E. 7) 


(E. 8) 


where the integration contour runs along the real axis in the conplex 

T plane. [Note that the integrand of (E.8) decreases rapidly as | T j 

2 

increases; had we chosen +T in (E.7) the integrand of (E.8) would 
have increased rapidly as |t| increased.] Assuming F(T) is suffi- 
ciently slowly varying, and extending the limits of the T integration 
to T = ± 00 we may write 

— Or~> 


- F ( b>0 V/r 


(E.9) 


since the exponential peaks the integrand of (E.8) about T = 0 . 
Actually the limits at T = ± » are not required; the result of (E.9) 
will be obtained as long as T covers the range -3 5 T < +3 and 
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F(T) is approximately constant over the same range. Combining (E.8), 
(E.9) we obtain finally 


I - 



i-V"(K) 


f(w 0 ) e 


4 0 (V d ) 


(E. 10) 


This result is the first term of an asymptotic expansion for I , more 
terms of which may be obtained (at least in principle) by keeping more 
terms in the Taylor series expansion for <(>(W) in (E.5). 

Returning to the original integral (E.l), we note that the change 
of variable used in (E.7) for the approximate case suggests that we use 
the change of variable 

0O) - 0 (i<\£) = - 7 2 (E. 11) 


(with T real as before) in the exact case (E.l). The path of constant 

phase defined by (E . 11) for -“ < T £ + 00 is the path of steepest 

2 2 

descent. [The path obtained by vising +T in place of -T in (E.l 1) 
is the path of steepest ascent.] To obtain an analytic expression for 
I we must first invert (E.ll) to obtain W as a function of T . 

Then this W(T) is used in F(W) and <^W to give 

I = / f[vj(t)\ (E>12) 

where the integration contour runs along the real axis in the complex 

T plane. Presumably the integrand of (E.12) could be written as a 

-T 2 

power series in T times e and the integration could be performed 
analytically term by term. 

In actual practise the required inversion of (E.ll) may be 
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difficult and one often resorts back to using (E.10) as an approximation 
for I . In this case some comment should be made on the range of 
validity of (E.10). To this end we consider values of W that are near 
W Q and on the path of steepest descent in the range defined by 

(w-w o f < / (E. 13) 

Then (E.10) is a good approximation to I provided: 

(i) Any singularities of 4>(W) occur for values of W far from 
those defined by (E.13). 

(ii) F(W) varies negligibly, i.e., 

L « / 

%.) (E. 14) 

for values of W within the range defined by (E.13). 

(iii) The terms in the Taylor series expansion of <p (W) past those 
terms used in (E.6) must be negligibly small. Using the 
first term dropped to represent the error, we must have 

0 W (K) 

« / 

U"(W 0 ) (E.15) 

which, combined with (E.13), requires 

Vz &'"(%) « 3 (E . 16) 

**• J 

We now apply the above results to some specific cases as used in 
Sections I and II. 
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Branch -cut integral in the k plane. 


We consider Eq. (2.58) of Part I which is of the form 


i = j f(£) e 


0(A) - 



(E.17) 


This form of <J>(k) occurs generally in the branch-cut integral that 
arises in a spatial problem (with a zero-order Maxwellian velocity 
distribution) when the inverse transform of f^(k,oj) [or £j(k,w) , 
etc.] is computed, so the following results should be of general 
interest. The integrand of (E.17) has three saddle points given by the 
three roots of 



2 

2 (CO-& eg) 



(E. 18) 


where k g represents a saddle point. The paths of steepest descent 
defined by (E.ll), i.e., by 

0(/O - 0 (A s ) = - t z 


are given by the roots of 


where 




k * 
£ = 



r 2 • + i - o 

(E.19) 
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and T runs along the real axis from T = - 00 to T = + °°. The 
quantity Z is the characteristic dimensionless parameter associated 
with (E,17). In terms of Z , the saddle points are given by the three 
roots of 

fir* =(^ y ?) 3 = 2^ 2 e 2 (E. 21 ) 

The paths of steepest descent are calculated as follows. We 
choose a value for Z and then select a specific saddle point from the 
three given by (E.21). We then solve the cubic equation (E.19) analyti- 
cally which gives three values of K for each value of T ; two of 
these roots are points on the desired path of steepest descent while 
the third root is extraneous [for T = 0 the three roots of (E.19) are 
k s , K s> and " » the latter being extraneous]. Then by numerical- 

ly evaluating the expressions for the two desired roots for several 
values of T (- 00 < T < +°°) we are able to plot the desired path of 
steepest descent. 

We now display the results of our calculations. In Fig. E-l we 
show the paths of steepest descent for Z - ± 1 . Note that according 
to (E.21) any specific saddle point for Z = +1 changes its phase by 
2 tt/ 3 radians as Z changes to Z = -1 . Thus the resultant distribu- 
tion of saddle points and paths of steepest descent are the same for 
Z = +1 and Z = -1 . To see how the contours actually shift as Z 

goes through zero, we have performed the following calculation. We let 

i 0 

Z = e with 0 <_ 8 £ 2tt so Z takes on values of points on a unit 
circle centered at the origin in the complex Z plane. We then calcu- 
lated the saddle point contours for several significant values of 0 
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with the results displayed in Fig. E-2. The contours for 0=0°, 180° 
are of course the same as the contours for Z = ± 1 given in Fig. E-l. 
Note that as 0 increases from 0° to 180° all three saddle points 
rotate rigidly in the counterclockwise direction. 

Of special significance in Fig. E-2 are the cases 0 = 45° and 
0 = 135° where in each case one of the contours violently changes shape. 
Consider the case 0 = 135° . Then suppose we had been using contour 2 
(which has one endpoint at infinity) to evaluate some integral for 
0° <_ Q < 135°. We note that as 0 crosses 0 = 135° the endpoint at 
infinity of contour 2 is diverted back to the origin and therefore con- 
tour 2 is no longer suitable for evaluating the integral; we must there- 
fore change to contour 1 as 0 crosses 0 = 135°. This behavior is 
known as Stokes phenomenon and the lines in the Z plane defined by 
Z = r e*^ 4 ^ , r e^ 371 ^ 4 ^ (r real and 0 £ r £ °°) are called Stokes 
lines. Note that for 0° <_ 0 <_ 180° contour 1 in Fig. E-2 never 
exhibits Stokes phenomenon and may therefore be used for all values of 

0 within the given interval. 

Thus, in Part I, in evaluating (2.58) we have chosen to use contour 

1 of Fig. E-2. The first term of the asymptotic expansion, (E.10), then 
gives the result (2.61). We now estimate the errors involved in using 
(E.10) as discussed with (E.13) - (E.16). 

First we construct a graph to aid in establishing over what values 

of K , F(K) [as defined by (E.17)] must be essentially constant. We 

2 

scale (E.19) appropriately by dividng it by Z to give 


Ar 
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/r_ + i 


+ c T 






(E. 22) 
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where 




(E. 23) 


We consider only Z real so T is real. Then in Fig. E-3 we have 
plotted "contour 1" of (E.22) for Z > 0 . (For Z < 0 the contour is 
just mirrored about the imaginary axis.) With the aid of Fig. E-l one 
may easily estimate the range of K over which F(K) should be approx- 
imately constant. The dominant portion of the integral comes from 
| T | ~ 2 say [see (E.12)]. Thus we set | T J = 2 . Then given a value 
for Z we calculate | T j from (E.23) and find the two corresponding 
points on the contour in Fig. E-3. The segment of the contour between 
these two points is the range of K over which F(K) should be 
approximately constant. In particular for Z = 1 (cutoffs at | T | = 2) 
a large portion of the contour is needed whereas for large Z , say 
Z = 27 (cutoffs at |t| = 0.66) a smaller portion is needed. 

Upon examining F(k) of (2.58) of Part I with the aid of Fig. 

E-3, one notes that the only poles of F(k) that may be troublesome 
are those in the infinite sequence of poles that come in along the 45° 
line (see Appendix C) in the k plane (and therefore also in the K 
plane). But | (to-tu ce ) /ka | is very large at these poles so the exponen- 
tial in (E.17) should make any pole effects negligible. The slowly 
varying behavior of F(k) is roughly like 1/k 3 (for to/k « c) , 
so using F(k s ) instead of F(k) in (E.9) does of course introduce an 
error, which becomes negligible for large Z . Thus, preferably, we 
should have |z| » 1 . 
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Lastly we find that the requirement (E.16) 

— ; — » / 

|v7^> /// (^ r )| 

by direct calculation from (E.19) reduces to the requirement |z| » 1 . 
This means that the asymptotic expansion method used to derive (E.10) 
breaks down when |z| becomes less than a characteristic length 
|a/(w-u) )| , in which case (2.58) would have to be evaluated numeri- 

cally (along the path of steepest descent for convenience) to obtain an 
accurate result. 

Thus, in Part I, (2.61) is an accurate result for |z| >> 1 and 
it is still a good order- of -magnitude estimate of (2.58) for |z| ~ 1 . 
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K plane 
3 + 


FIGURE E 




. SADDLE-POINT CONTOURS OF (E.I7) FOR Z«±l. 

Shown are the saddle points of (E.17) as given by (E.21) 
and the paths of steepest descent as given by (E.19) 
for 2*±1. This figure holds for Za+1 or 2«-l although 
any specific saddle point rotates through 2ir/3 radians 
when 2 goes through zero, as is shown in the next 
figure. 
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FIGURE E-2. EVOLUTION OF THE SADDLE-POINT CONTOURS OF (E.17) FOR 
Z ■ e l6 AS 0 VARIES FROM 0° TO 180°. 

Saddle points and paths of steepest descent are as given 
by (E.21), (E.19) respectively. Note Stokes phenomenon 
that occurs when 6 * 45° and 0 = 135° . 




FIGURE E-3. SADDLE-POINT CONTOUR OF (E.17) IN THE K PLANE WITH T 
SPECIFIED ALONG THE CONTOUR. 

Shown is the path of steepest descent of <E.17) in the 
K («=K/2 2/ ^) plane for 2 real and with values of T (sT/Z 1 ^ 3 ) 
specified along the contour [see (E.22) ,(E,23) ]. This 
plot is useful in estimating the range of the variable K 
over which F(K) should be approximately constant. 
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Phase-mixing integrals in the v plane. 

£ 

We consider Eq. (2.80) of Part I and Eq. (4.4) of Part II, both of 
which are of the form 


T = f F(v z ) 

- ‘ (Fr*) 2 - 


(E. 24) 


This form of <j>(v ) occurs generally in a spatial problem (with a zero- 
z 

order Maxwellian velocity distribution) when the free-streaming part of 
fjfz.v.t) is integrated over v z to give J_j(z,t) (or Ej , etc.). 
Thus again the following results should be of general interest. 

We note immediately that <j>(v _) in (E . 24) may be obtained from 
<Kk) in (E . 17) by the transformation 


— 


Co- 



(E . 25) 


so all of the graphs for the previous case apply to the new variable 
[(w-u) )/v ] . Of course we want graphs in terms of v directly, so 

we would have to invert the previous graphs point by point according to 
(E.25) to obtain the desired graphs. Instead, we shall write out the 
appropriate equations in terms of v and solve them, although at any 
stage we can check our results by comparing them to the previous case 
using (E.25). 


The integrand of (E.17) has three saddle points given by the three 
roots of 




(^-6c/ ce )? < 


7T 


(E . 26) 
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where v g is a saddle point. The paths of steepest descent defined by 
(E. 11) , i.e. , by 

- 0(V 5 ) = - T* 


are given by the roots of 


Y%V 




= o 


where 


V = 


_v 

C2 


3 3 


2^-^ce) 


C2_ 


(E.27) 


(E. 28) 


and T runs along the real axis from - » to + <» . Once again the 
quantity Z is the characteristic dimensionless parameter. In terms 
of Z the saddle points are given by the three roots of 




3 



(E.29) 


We note for reference that in terms of K = kz as in (E.20), the 
transformation (E.25) is 


fC 


"V" 


(E.30) 


The paths of steepest descent are calculated from (E.27) by the 
method used in the previous case. As earlier, given Z and a specific 
choice of V g two of the roots of (E.27) are points on the desired 
path of steepest descent while the third root is extraneous [for T = 0 
the three roots of (E.27) are V g , V g , -2V s , the latter being 
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extraneous ] . 

In Fig. E-4 we show the saddle-point contours for £ = +1 and 
£ = -1 . Note that the contours for the two cases are not the same 
but are mirror images about the real axis of each other. In Fig. E-5 
we let £ = e^ with 0 £ 6 £ 2tt as earlier. We see that as 0 goes 
from 0° to 180° all three saddle points rotate rigidly in the counter- 
clockwise direction, but only through ir/3 radians. Also we find 
contours 2 and 3 exhibit Stokes phenomenon at 0 = 45° and 0 = 135°, 
respectively. Note that contour 1 does not exhibit any Stokes phenome- 
non for 0° < 0 180° and therefore this is the contour we choose to 

use. 

To aid in establishing over what values of V, F(V) must be 
essentially constant, we construct a graph similar to Fig. E-3. We 
scale (E . 27) appropriately by dividing it by £ to obtain 


V + V 





o 


where 


V 

~ 7 ~ 


~\T 


T 
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'/3 


(E . 31) 


(E. 32) 


As earlier, we consider only 2 real so T is real. Then in Fig. 

E-6 we have plotted "contour 1" of (E.31) for £ < 0 (for £ > 0 the 
contour is just mirrored about the real axis). Also, as earlier, we 
may set | T | = 2 and then given £ we can compute |t| from (E.32) 
to obtain the cutoffs on the contour in Fig. E-6. 
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Thus using contour 1 and the first term of the asymptotic expansion 
(E.10), (2.80) of Part I becomes (2.82). Similarly in Part II, Eq. 

(4.4) becomes (4.8). Roughly the integrand of the first case goes 

3 

'v 1/V while the integrand of the second case goes *v> 1/V . Thus with 

the aid of Fig. E-6 we conclude that F(V) in each case is sufficiently 
constant provided | 2 | » 1 but that even if |z| « 1 the errors 
involved should not change the results greatly. 

Lastly we note that the requirement (E.16) 




[V'rvV) 




1? 



» / 


by direct calculation from (E.24) reduces to the requirement |z| » 

0.12. Thus the asymptotic expansion method used to derive (E.10) 

breaks down in this case only when |z| becomes much smaller than a 

characteristic length |a/(u>-w )| . 

cc 

Thus the free -streaming term in (2.85) of Part I and the echo 
saddle-point contribution (4.8) of Part II are accurate for |z| » 1 
but still good approximations for Jz| a 1 ; these results lose their 
validity only for |z| significantly less than one. 
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V plane 




FIGURE E-4. SADDLE-POINT CONTOURS OF (E.24) FOR 2»±1. 

Shown are the saddle points of (E.24) as given by (E.29) 
and paths of steepest descent as given by (E.27) for 
£■+1 and 2«-l. 
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FIGURE E-5. EVOLUTION OF THE SADDLE-POINT CONTOURS OF (E.24) FOR 
2 » e ifl AS 6 VARIES FROM 0° TO 180°. 

Saddle points and paths of steepest descent are as given 
by (E.29), (E.27) respectively. Note Stokes phenomenon 
that occurs when 8 ■ 45° and 8 * 135°. 






FIGURE E-6. SADDLE-POINT CONTOUR OF (E.24) IN THE V PLANE WITH T 
SPECIFIED ALONG TOE CONTOUR. 

Shown is the path of steepest descent of (E.24) in the 
V (-V/Z 1 ^ 3 ) plane for Z real and with values of T (^T/Z 1 ^ 3 ) 
specified along the contour [see (E.31) , (E.32) ]. This 
plot is useful in estimating the range of the variable V 
over which F^V) should be appro xitoately constant. 
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